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The Scientific Advisory Committee of the International Atomic Energy Agency in June 1959 recommended 
that the Agency undertake the promotion and exchange of information on plasma physics and controlled nuclear 
fusion research. As a specific proposal, the Committee recommended the establishment of an international scientific 
journal in this field. In fesponse to this proposal the Agency now inaugurates “NUCLEAR FUSION.” 

We believe this new journal will contribute to the improvement of cooperation between nations in the difficult 
task of finding ways to extract energy from controlled nuclear fusion. Even if the earnest and concerted efforts 
devoted to this task do not lead immediately to practical energy sources, the endeavor will greatly expand the 
knowledge of plasma -physics besides broadening basic scientific knowledge. The effort will also undoubtedly lead 
to the invention of new and important peacetime devices. 

It is an unusual fact that, although the universe—as it is now understood—consists almost entirely of 
plasma, scientists have made progress only relatively recently in understanding the properties of this “fourth state 
of matter.” As the pages of “NUCLEAR FUSION” show, the properties of plasma are not simple nor are they 
easily investigated. 

All those who read and use “NUCLEAR FUSION” are invited to participate in the journal. Their comments 
and suggestions, as well as their manuscripts, will be welcomed. To those who have helped us in creating the 


journal, I wish to express, on behalf of the International Atomic Energy Agency, my sincere appreciation. 


Le Comité consultatif scientifique de Agence internationale de l’énergie atomique recommanda en juin 1959 
que l’Agence encourage la diffusion et Péchange de renseignements sur la recherche en matiére de physique des 
plasmas et de fusion nucléaire contrélée, notamment par la publication d’une revue scientifique internationale 
sur ces questions. Conformément a cette recommandations, l’Agence presente maintenant le premier numéro de 
sa revue «FUSION NUCLEAIRE». 

Nous sommes convaincus que cette nouvelle publication contribuera a resserret la coopération entre les pays 
qui s’efforcent de résoudre le difficile probleme de la production d’énergie a partir de la fusion nucléaire contrdlée. 
Meme si leur effort intense et concerté n’aboutit pas immédiatement a la découverte de sources d’énergie utilisables, 
I n’en contribuera pas moins 4 étendre considérablement les connaissances en matiére de physique des plasmas 
st a approfondir de plus les notions fondamentales de la science. Cet effort se traduira certainement aussi par 
Pinvention de nouveaux et importantes dispositifs aux utilisations pacifiques. 

Si Punivers — comme on le considére a present — se compose presque entiérement de plasma, il est assez 
ttonnant que les hommes de science n’aient fait que des progres relativement recents dans la comprehension des 
sropriétés de ce qui peut étre considéré comme le «quatriéme état de la matiére». Comme il apparaitra a la lecture de 
(FUSION NUCLEAIRE», ces propriétés ne sont ni simples ni d’une étude facile. 

Tous les lecteurs de «FUSION NUCLEAIRE» sont invités y collaborer, a adresser 4 l’éditeur leurs observa- 
ions, leurs suggestions ainsi que les communications qu’ils voudraient y voir inséré. A tous ceux qui nous ont aidés 


mettre au point cette revue, je tiens 4 adresser, au nom de l’Agence internationale de l’énergie atomique, mes 


emerciements les plus sincéres. 


B wioHe 1959 roqa HayuHo-KOHCYIbTAaTUBHbIM KOMM TeT MexkayHapofHOTo areHTCcTBa 10 aTOMHOM 9HeEprMuu pe- 


cmoco6cTBy!0O tMe OOMeHY MHcbopMalMen 10 MccueOBaHUAM B oOmacTu 
WWOmReHMA KomMuTeT PCKOMECHAOBA 


KoMeHyOBaa ATeHTCTBY HPMHATb MepsI, 


cbu3uKu WIa3MbI U KOHTpPOMMpyeMoro CUHTE3Aa. B KauectBe KOHKpeTHOTrO Tipe 


BbIMyCK MexKAyHapOAHOTO HayuHoro wypHasa B 9Toi ob6macTu. B cOOTBeTCTBUM C 9TMM MpeAIO-KeHUeM ATreHT- 


cTBO B HacToaujee BPeMA TpMcTylaer K MU3AaHMIo WYPHasa «ATePHbIM CUHTe3». 


Mer curaem, 4TO 9TOT HOBbIM wypHam Oyger cofelicTBOBAaTb YJy4UIeHuM1o coTpyqHMuecTBAa Mex*K]y cTpaHaMu 


B PeWIeCHMM TPYAHOM 3aqaum HAaxXOMRMeHUA MyTeM NONy4eHUA 9HeEpruw U3 KOHTpOIIMpyeMoro AyepHoro CMHTe3a. 
Tlaxe ecmm cepbe3Hble M COrmacoBaHHble yCMJIMA, HampaBsIeHHble Ha pelleHve 9TOM Baqauu, cpa3sy MU He NO]- 
BelyT K UpakTM4ecKMM MCTOUHMKAM SHeEprMu, ITA WOMbITKAa, MOMMMO yIJIyOJIGHMA OCHOBHbIX HayUHbIx 3HaHMn, 
3HAYNTENBHO pacluMpuT No3sHaHMA B OOMacTM CbMSUKM Wia3MbIl. QTO ycusmMe HeCOMHeHHO npMBpemeT Takwe K 


HOBLIM M BasKHbIM TeCXHMYeCKUM M300peTEHUAM MMPHOTO XapakTepa. 

HeoOnikHOBCHHBIM CbaKT 3aKIOUAeTCA B TOM, UTO, XOTA KOCMOC — KaK OH MOHMMaeTCA B HacTOANICe BpeMA 
— TOUT NONHOCTHIO COCTOMT M3 TWa3MbI, YUCHbIe JIMUIb CpaBHMTeNbHO HeyaBHO WocTurmM Mporpecca B MOHU- 
M&HUM CBOVICTB STOO «<UeTBeEPTOTO COCTOAHMA MaTepMM»>. Kak l0Ka3bIBaloT ONyOMKOBaHHbIe B WYypHasie «Anep- 


Hb CUHTe3» CTATbM, Uma3Ma OONaqaeT OTHIOAb He MpOCTbIMM CBOMCTBaMM, M HeWerKO TloAAarwoTcn ITU cBOM- 


cTBa MCCIeOBaHMrO. 
Bee unratemm x%ypHana «AyepHbmt cMHTe3» IpMrsiawmiaiorcA YyuacCTBOBaTb B ero M3qaHun. Ux 3sameuaHMa, 
mpeqloxeHMA, a TakxKe cTaThu OysqyT upMBercTBoBaTbcA. OT uMeHu MexmyyHapogHoro areHcTBa MO aTOMHOM 


SHeEpruu A XOTEeT ObI BbIPAaASMTb MOO MUCKPCHHIOLO 6naroqapHOcTb BCeM TeM, KTO MOMOr HaM B co3qaHuu 3STOTO 


wKypHasa. 


En junio de 1959, el Comité Consultivo Cientifico del Organismo Internacional de Energia Até6mica recomend6 
que el OIEA fomentase la difusién y el intercambio de informacién en el terreno de las investigaciones sobre la 
fusidn nuclear controlada y la fisica del plasma. Ademas, el Comité recomend6 concretamente que se iniciase la 
publicacién de una revista cientifica internacional que tratase de estas cuestiones. Como consecuencia de esta 
propuesta, el Organismo lanza hoy el primer numero de “FUSION NUCLEAR”. 

Creemos que esta nueva revista contribuira a acentuar la cooperacion entre las diversas naciones en la dificil labor 
de descubrir procedimientos para la obtencidén de energia partiendo de la fusidn nuclear controlada. Aun cuando los 
mas setios y atmonizados esfuerzos que se dediquen a esta labor no conduzcan inmediatamente a disponer de fuentes 
energéticas, la empresa extendera muy considerablemente el conocimiento de la fisica del plasma ademas de ampliar 
el caudal de conocimientos cientificos basicos. Estos esfuerzos se traduciran sin duda alguna en la invencidn de 
nuevos e importantes dispositivos utilizables en tiempo de paz. 

Es realmente sorprendente el hecho de que, aunque el universo — segun el concepto que hoy se tiene del mismo — 
esta formado casi enteramente por plasma, hasta hace relativamente poco no hayan logrado los hombres de ciencia 
progtesar en el conocimiento de las propiedades de este ““cuarto estado de la materia”. Como se refleja en las paginas 
de “FUSION NUCLEAR’’, la investigacién de las propiedades del plasma no es labor facil ni sencilla. 

Quienquiera que lea y haga uso de “F USION NUCLEAR” queda invitado a colaborar en la revista. Sus observa- 
clones y sugerencias, asi como los originales que envie, seran bien acogidos. En cuanto a quienes nos ayudaron a 
lanzar esta revista, deseo expresarles aqui, en nombre del Organismo Internacional de Energia Atomica, mi sincero 
agradecimiento. 


Senbag 


Director General 
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PLASMA OSCILLATIONS (1* 


I. B. BERNSTEIN ** and S. K. TREHAN *** 


Project MATTERHORN, PRINCETON UNIversity, Princeton, NEw Jurszy, USA 


This article is the first part of a review of the theory of wave phenomena in plasmas. The basic kinetic theory is 
developed taking into account coulomb collisions. From this magnetohydrodynamic equations are derived. It is 
demonstrated that it is legitimate to deal with a closed subset of these equations either in the limit of collision- 
dominated phenomena, or in the limit whete the effective phase velocities of the phenomena of interest ate much 
greater than thermal speeds. In case of collision-dominated plasmas the theory of transport coefficients is discussed. 

These equations are then applied to an extensive treatment of small amplitude wave phenomena in plasmas. 
A discussion of the dissipative effects on hydromagnetic waves is given. Hydromagnetic waves are also considered 
from the Chew, Goldberger and Low theory. Longitudinal and transverse oscillations in current-carrying plasmas 
are also discussed. 

Oscillations of a cylindrical plasma are considered and the phenomenon of ion cyclotron resonance is discussed. 
The possibility of radiation by plasma oscillations by a uniform sphere is exhibited. Some general results on the 
stability of longitudinal electron oscillations in non-uniform plasmas are given. 

A brief treatment of large amplitude electron oscillations is given and the breaking of these oscillations as a 
dissipative mechanism for the organized plasma motion is discussed. 

Part II of this paper, to appear in a future issue of this journal, will be devoted to the discussion of piasma 
oscillations directly from the kinetic theory. 
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I. Introduction 


In this article we shall attempt to present a systematic 
survey of the theoretical knowledge of wave phenomena 
in plasmas. We shall use the term “plasma” to describe 
a fully ionized gas composed of positive and negative 
ions of effectively zero net charge density. For the sake 
of convenience we limit our attention to a plasma com- 
posed of only two species (positive ions and electrons) ; 
the extension of the results to gases consisting of more 
than two species is straightforward and will not be con- 
sidered here. Some of the results are also applicable to 
conducting liquids and gases not fully ionized. 

Astrophysical examples of such systems ate the ion- 
osphere, the solar corona, stellar interiors and the at- 
mosphere of hot stars, the H,, regions of interstellar space, 
and cosmic radio sources. Plasmas are also of fundamental 
interest in the proposed schemes for controlled thermo- 
nuclear fusion, plasma propulsion of space ships, micto- 
wave generators, and gas discharge devices. 

One distinctive feature of plasmas is the special nature 
of collisions, in contrast to the two-body collisions 
familiar in the conventional theory of neutral gases. 
Because of the long range nature of the Coulomb forces 
which govern the encounters of charged particles, co- 
operative phenomena are of major importance in describ- 
ing the dynamics of a plasma. Also these long range col- 
lisions lead to a diffusion in velocity space and are most 
appropriately described by the Fokker-Planck equation. 
The elastic or inelastic two-body encounters of ions of 
electrons with neutrals or perhaps with the core of an 
incompletely ionized atom are usually negligible in 
plasmas and will be considered, effectively, only when 
dealing with dissipative effects in collision-dominated 
plasmas. 

There are two levels of description of plasmas which 
have received major consideration during the past dec- 
ade. The first of these is based on the direct solution of 
the Boltzmann transport equation —usually referred to 
as the kinetic theory description. This method has the 
inttinsic merit of providing a more complete description 
of the dynamics of a plasma. However, this procedure is 
mathematically more difficult to handle. The second 
method is based on using a closed set of moment equa- 
tions to characterize the behaviour of a plasma. This de- 
scription, usually referred to as the hydrodynamic or 
hydromagnetic description, though of limited applica- 
bility (in a sense to be discussed in detail later) has the 
merit of relative mathematical simplicity and is quite 
helpful in the physical understanding of the various 
phenomena involved. This method is applicable either 
when (1) collisions dominate so that one has local thermo- 
dynamic equilibrium and can apply the Chapman-Enskog 
expansion technique, or (2) the phase velocity of the 
wave is much larger than the characteristic thermal 
speed in a collision-free plasma. The latter will be re- 
ferred to as the “Low Temperature Approximation.” 

For plasmas of interest under astrophysical and lab- 
oratory conditions quantum mechanical effects ate un- 
important in considering most of the dynamics of plasmas. 

This work, therefore, divides itself naturally into two 
patts. In part I the kinetic theory description is formulated 
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and its relation to the hydrodynamic description est- 
ablished. The resulting hydromagnetic equations are then 
applied to a variety of wave phenomena. Part I of this 
atticle, to appear in a later issue of this journal, will be 
devoted to the discussion of wave phenomena based 
directly on the kinetic theory. 

The theory described here is concerned mainly with 
the small motions about a steady state of the plasma. The 
general case of large motions is considered only briefly 
in Section V. 

Section HI deals with oscillations in unbounded 
homogeneous plasmas, admittedly a severe idealization. 
However, these treatments do serve to give some in- 
sight into the basic physical processes and the mathe- 
matical methods. The neglect of boundaries would imply 
that the description is valid only for wavelengths much 
smaller than the characteristic length of the system and 
for time intervals short enough that the disturbance has 
not reached the boundaries and been reflected and trans- 
mitted. 

Unfortunately, because of the inherent difficulty of 
establishing in the laboratory plasmas uncomplicated by 
boundary eftects, impurities, etc., there has not been, as 
yet, much experimental verification of the theoretical 
investigations. 


Il. The basic equations 


By and large one is interested only in the knowledge 
of some macroscopic properties of the plasma like the 
mean particle density, mean velocity, etc., rather than 
in the detail afforded by the description of the individual 
motion of the patticles. It is, therefore, appropriate to 
characterize a given class of particles by a distribution 
function f(x,v, 2) such that fd?x d°v tepresents at time ¢ 
the probable number of particles with velocities between 
vandv -- dv, and with positions between x and x + dx. 
It is possible to relate f to the more recondite density in 
the phase space of all the particles, but a discussion of this 
matter is beyond the scope of this article; suffice it to say 
that the introduction of f is compatible with the follow- 
ing notions about particle interactions in plasmas. 

The particles, being charged and possessing no in- 
ternal degrees of freedom at the range of energies of 
interest here, interact primarily through the electro- 
magnetic fields. In particular, when two particles are 
close together it is their Coulomb interaction which 
ptedominates. Consider then such two-body Coulomb 
encounters. They can be characterized, in their center 
of mass system, by the impact parameter and the relative 
velocity. The collisions can then be conveniently divided 
into three classes. 

The first of these corresponds to impact parameters 
less than the mean impact parameter for 90° deflection, 
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where 7, and Z, are the charge numbers of particles in- 
volved, ™, and m, are theit masses, <v?> represents their 
mean square relative velocity, and e is the electronic 
charge. It is usually legitimate to describe such en- 
counters by the familiar two-body collision integrals of 


the kinetic theory of gases when writing the kinetic 
equation for the distribution function. In most cases of 
Interest these collisions are negligible. 

The second class comprises the many-body encounters 
which can be regarded as composed of a succession of 
uncorrelated small angle Coulomb deflections. These 
correspond to elementary Coulomb scattering events 
characterized by impact parameters greater than the dis- 
tance of closest approach b, but less than the electron 
Debye length 

1/3) a <*>] 4 
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where m is the electron mass, <v2) is the electron mean 
Square velocity, and N is the electron number density. 
These encounters ate most appropriately described by 
the Fokker-Planck form for the collision integral in the 
Boltzmann equation. 

When there are many particles in a sphere of radius 
me, (i.c., A, N >> 1), density fluctuations over distances 
ros A, ate suppressed by macroscopic electric for- 
ces. Thus the third class of uncorrelated two-body en- 
counters, those with impact parameters greater than aps 
do not occur per se.These are replaced in the kinetic 
eqation by the correlated effect of many particles through 
the macroscopic electromagnetic fields described by the 
Maxwell equations; this leads to co-operative phenomena. 


1. BOLTZMANN EQuaTION 


The distribution function for each species of particles 
is governed by the Boltzmann equation 


7) 
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where V and V,, represent the gradient operators in co- 
ordinate (x) and velocity (v) space respectively, F de- 
notes the force per unit mass, and (0//O/)con. represents 
ithe time rate of change of f resulting from collisions. It 
| may be remarked here that there is a separate Boltzmann 
jequation for each species of particle. We now discuss 
|the Fokker-Planck form for the collision term. 


1.1 Fokker-Planck Form of the Boltzmann Equation 


| It has been demonstrated by Jeans that when pasticles 
linteract through long range forces, the cumulative effect 
lof weak deflections resulting from the relatively distant 
ncounters is more important than the effect of occasional 
atge deflections. The cumulative effect of weak de- 
ections is most appropriately taken into account by 
titing the collision term in the Fokker-Planck form. 
his reduction is possible only if there exist time inter- 
als A¢ long enough for the particles to suffer a large 
umber of collisions but small enough for the net mean 
quatre change in velocity, <|Av|?>, to be small com- 
ated with the mean square velocity, and to be of 
adet Az. 

Let P(v; Av) denote the probability that in the time 
nterval Az a particle with velocity co-ordinate v under- 
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goes a displacement Av. The distribution function is 
then given by 


L(X, van 
\f (x, v— Av, #¢— AAP (w— Av; Av) BAv. (4) 


We may remark here that in expecting this integral 
equation to be true we are actually supposing that the 
course which a particle will take depends only on the 
instantaneous value of its physical parameters and is 
entirely independent of its past history. In general prob- 
ability theory a stochastic process which has this char- 
acteristic, namely that what happens at a given instant 
of time ¢ depends only on the state of the system at 
time 7, is said to be a Markoff process. 

Since Af and Av are both assumed small, using 
Taylor’s expansion we can write (4) in the form 


Tay. = 
3 af 
FCS Vea) LOsAy) —At~ P(v; Av) 
” | ” | 1 
Av 30 P vs Ay) ailing P(w; Av)| +3 AvAv: (5) 
oy 0) OP O3E es 
seas Pv Av) | nes foal eerenen aay. 
Using the fact that 
| P (wv; Av) d§ Av =1, (6) 


Eq. (5) gives for the time rate of change of f resulting 
from the cumulative effect of small deflections: 


of 0 is ; 
i [Kav 55d Av» , (7) 


where we have retained terms only to lowest significant 
order and 


<Avy=JSf (x, v, 4 P(v; Av) Ava? Av, (8) 
<Av Avy =|\f(«, v,) P(v;4v) AvAvdAv. (9) 


Eq. (7) represents the collision term in the Fokker- 
Planck form. 


Consider the two-body encounters of particles of type 7 
with those of type / and let o,, (0,v) denote the differen- 
tial cross section for scattering of their relative velocity 
vector v = v,; —v, through an angle 0. Then the as- 
sociated momentum transfer cross section is given by 


0; 0) = 2x 504 0,0) (1 —cos 6) sin 649. (10) 


In terms of QO, retaining only the dominant terms in the 
mean square scattering angle, one can write for the con- 
tribution of such two-body encounters to stochastic 
averages, when polarization effects are negligible, 


<Av);; (11) 
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1.2 Fokker-Planck Equation for Coulomb Force 
For the case of Coulomb collisions, the differential 
scattering cross section is given by 
aoe 5 
2,2; 2 (me m;) 1 


‘sin? (0/2) 
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The integral occurring in Q will diverge at the lower 
limit of integration (9 = 0). However, the small angle 
deflections correspond to scatterings with very large 
impact parameters, and the divergence arises from the 
long range nature of the Coulomb force. This divergence 
is eliminated when we take into account the shielding 
that arises from the polarization of charge surrounding 
the scatterer. The polarization screens the scattering 
particle and provides a natural cutoff on the maximum 
impact parameter of the order of the Debye length J,,. 
We thus take for the minimum scattering angle 


2m, mM, v" 


= (m,-+ 13) 
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Using Eqs. (13) and (14) into Eq. (10) we obtain 
hy Ship a aa “| z 1 
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(14) 
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xo M1, M; <0» yes (15) 
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Using the foregoing results in the Fokker-Planck equa- 
tion we obtain after an integration by parts what is called 
the Landau form for the collision term: 


a = View [2a ZP 2; et (72; + ;) 
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It may be noted that, apart from the slowly varying 
logarithmic term, in order of magnitude the collision 
terms are inversely proportional to <v*>,,, the mean 
square relative kinetic energy. Thus in the limit where 
this factor is appropriately large, collisions are negligible. 


1.3 The Moment Equations 


The particle density and the local mean velocity are 
defined as 


N= a) ay (17) 


and 
1 3 
A, fv, Avaev. (18) 
In terms of these moments we can define some of the 
other relevant macroscopic quantities 90 


= > m INV = mass density, (19) 


Uy= => m Nu = velocity of thecenter of mass, (20) 


* We shall adopt the convention that whenever 7 appeats after a 
summation sign, it is the mass of any of the species of the plasma; 
otherwise # denotes the electron mass and M, the proton mass. 
The electron charge is taken to be —e while that of the proton, + e. 
Similarly, we will denote the electron or ion cyclotron frequency 
whereas We and @; will denote the electron and ion cyclotron fre- 
quency respectively. 
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= > e N= charge density, (21) 


\= > e Nu = current density, (22) 


where the summation is to be catried over all the species 
of the system. The random velocity with respect to the 
center of mass is defined to be 


w=v— Up. (23) 
It is convenient to define the other relevant quantities 
with respect to the center of mass velocity. Thus 


= mean velocity of (24) 
each of the species 
with respect to the 
center of mass, 


= x Wik, V,2)a Vv 


7 wwf (x,v,/)d?v  =pressure tensor, (25) 


and 


O_7 | ww w/ (x, v, 4) dv = heat flow tensor, (26) 


where Q is a tensor of tank 3 completely symmetric in 
all its indices. We further define 


P= > P and Q=> Q. (27) 
It is clear that 
S MIN U =0. (28) 
Note that 
i 5 ew 2 
3 Trace P = aw Ks Va? a aN. (29) 


represents the mean kinetic energy relative to the center 
of mass and is referred to as the internal energy, or the 
thermal energy of the gas. This suggests that we can 
define the generalized temperature by the relation 


eee i= : NG: (30) 


2 
where 9 — AT, & is the Boltzmann constant, and T is 
the temperature. We shall always use © to denote the 
temperature in energy units. 

This definition of temperature is natural when colli- 
sions dominate, in which case the various particle dis- 
tribution functions are very closely locally Maxwellian, 
relative to the center of mass velocity, corresponding to 
the same temperature. Occasionally one encounters an 
alternative definition of temperature, where it is defined 
separately for each species of the plasma, as the mean 
kinetic energy of that species relative to its own mean 
mass velocity. The other moments can be similarly de- 
fined. This latter definition is particularly convenient 
when collisions are negligible and will be used in part II 
of this review. The two notions, of course, coincide if 
we consider the motion of a single species of the plasma, 
a picture which is reasonable in investigating high fre- 
quency waves where only theelectron motionis significant. 
In general, however, the two definitions are quite dif- 
ferent and care must be exercised to make sure which is 
being employed in any given context. 


The vector 


1 ; 
7m ww f (x, v, A) dev 


1 
q=;1Q= 
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(31) 


represents the flow of internal energy relative to the 
center of mass and will be called the heat flow vector. 

Since x and v are independent vectors, it is con- 
venient to write the Boltzmann equation (3) as 


ae 1 (2 
Df = apt V-WA+V,-(FA= (Sy (32 
Quite generally, we can write 
e 1 
P=5|E@)+4¥ x BE) Sa me es) 


where E and B denote the electric and magnetic fields 
tespectively, and is the gravitational potential. The 


electromagnetic fields, of course, satisfy Maxwell’s 
equations: 
4 ot 10E 10B 
=e es of a eo, \ —— ee 
ara i cere Vx B=. (34) 
V-E=42¢, and V-B=—0; (35) 


Gaussian units will be used throughout in this work. 
The gravitational potential satisfies the Poisson equation 


V2 == 45G o. (36) 


where G denotes the gravitational constant. 

In order to derive the equations governing the macro- 
scopic quantities, it is convenient to write the Boltzmann 
equation in terms of the random velocity w. This can 
de readily effected by the transformations 


TG, > 7 we), (37) 
0 od, 
oar (OE (38) 
eV Ui: (39) 
V,2Vy- (40) 


Che Boltzmann equation now takes the form 


i= Fy + V-[ (w+ Uns Jrveen=(3),, 


vhere 
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dU, 
dt 


wx B(x,/), (42) 
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nd 


d 0 
pe le ay 


dt 2) 


enotes the time derivative as we follow the motion of 
a element of the fluid. 

Multiplying Eq. (41) by 1, mw, and mww, and as- 
iming that f vanishes sufficiently strongly in velocity 
ace so that all the surface integrals vanish, we obtain 
yf the moment equations 


oN 


+: [Nu+NU,] = (44) 
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dU 
mN - oie —— IP + NUD 
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where (A)” denotes the transpose of the matrix A. We 
now sum Eqs. (44)—(46) over all the species of the 
system and, making use of Eqs. (19)—(22), (27); and 
(28), we obtain 


0) 
5p +V-@U)=0, (47) 


dU- 1 


and 


oP +V-(Q+UP)+P-VU+(P- VU)? 


+ 57(BX P—P x B)—J—eU) E+ 4U xB) 


1 C 


— +50 xB)g—ev)=(2p| (49) 


/ coll. 
where in writing Eq. (49) we have made use of the fact 
that 


DenU=>eNU—U,)=J—ev, 60) 


and we have dropped the subscript on Uy in Eqs. (47)—(49) 
as it is no longer necessary. If we take the trace of Eq.(49) 
and make use of the relation (31), we obtain 


5(229| +¥-(4+5200\+P:vU 


—J—ev).(E+ = x B}=0, (51) 


where the right hand side vanishes since collisions con- 
serve the internal energy of the system. 

We recognize Eqs. (47) and (48) as the well-known 
hydrodynamical equations of continuity and motion. 
Eqs. (47)—(49) are rather useless in their present form; 
they do not constitute a closed set of equations as they 
contain more unknowns than the number of equations. 
One could obtain an equation governing the heat flow 
tensor Q by taking the fourth moment of the Boltzmann 
equation; however, this equation would involve the next 
higher order moment of the distribution function and 
the system of equations is not closed again. 
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THE CASE OF STRONG COLLISIONS 


We now consider the case where collisions dominate. 
Then, to a high degree of approximation, the particle 
distribution functions tend to be locally Maxwellian and 
we may set 


fw) =f Gv) tO Gre. 
where 
i Mm 13/2 
VCR t) aN (x37) moee 
[v — Uy (, 4) ]? 
DO (%, 2) 


and @ is considered small compared to unity. The density 
N, temperature @, and mass velocity Up) are assumed 
to be slowly-varying functions of space and time such 
that their spatial variations over one mean free path, 
and temporal variations over one mean free time, are 
very small. This limiting situation forms the basis of the 
Chapman-Enskog theory of gases and leads to a trun- 
cation of the hierarchy of moment equations. 

The rigorous Chapman-Enskog theory for a plasma 
composed of two of more species, interacting with 
electromagnetic and gravitational fields, is quite com- 
plex. However, the smallness of the electron to ion mass 
tatio and the inability of a plasma to support large 
electric fields—referred to as the condition of quasi- 
neutrality —results in a substantial simplification of the 
Chapman-Enskog expansion technique. 

The condition of quasi-neutrality implies that the 
ratio of the net charge density to the charge density of 
each of the species is very small compared to unity. We 
shall systematically neglect both these ratios compared 
with unity in our discussion of the transport co-effi- 
cients. 

Consider for simplicity a plasma composed of N elec- 
trons and an equal number of protons of masses 7 and M 
and charges —e and +e respectively. Following the 
Chapman-Enskog procedure we can show that, within 
the approximation outlined earlier, the small quantities ® 
ate determined by the integro-differential equations (cf. 
Bas.G@)(16); (52); (93)): 
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The currents in velocity space I’, obtained by an appro- 
ptiate expansion of the Fokker-Planck Eq. (16), are 
given by 


8 


(52): 
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Equations (54)—(57) yield, correct to first order in 
@+ , the following closed system of moment equations: 
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or 
0 
o(2,+U.-V]U,+VP+eV¢ 
—*jJxB—eE=—V-(P—p), (60) 
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where p = (N, + N_)O = 2 NO is the isotropic part: 
of the stress tensor. In the foregoing equations the: 


quantities on the right hand side are determined in terms} 
of ® by the equations: 


J=>¢e\Pwwft OF, (62)) 
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The most comprehensive calculations of transport co- 
efficients for a fully ionized gas are due to Marshall using 
effectively the model outlined here. We shall simply 
state his results. Introducing the auxiliary electric field 
(65; 


ay 


1 il 
De Een Vo <x B 4. INe¥?P> 
and n, the unit vector along the magnetic field, we obtait 


J—ceU,=9,nn-D+o.nx(D xn)+o,n xD 
+x,nn-VO+4 %n x [(VO) Xn] +2x,nx VO, (66) 


where the various co-efficients are given by 
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Here w, denotes the electron cyclotron frequency and 7, 
the effective collision time, is given by 
s ) (5 i @? 

(22)3\0} 2NA@lnA 


where, it may be recalled, that 2/N is the total number 
density. 


(70) 


The heat flow vector is given by 


1 =—4,nn-VO—A.nx [(VO) x n] —A,n x VO 


a10n-D—»n x (D x n)—»,n XD, (71) 
where 
Ava 3.59 A, | 
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It is clear from Eqs. (66) and (71) that, in the presence 
f a temperature gradient, not only will there appear a 
ow of heat q, but also, a flow of electric current. Simi- 
tly, an electric field produces a flow of heat. According 
) the thermodynamics of irreversible processes, these 
put co-efficients are not independent as three of them 

re linked by the Onsager relations 


%=O(, +55), (75) 


here z takes the value 1, 2, or 3. 


In a representation in which the magnetic field is 
ong the z-axis, Eqs. (66) and (71) take the form 


J—eU,=c0-D+x-VO (76) 
id 
nere 
op — 0, 0 
o = | % 02 0 , (78) 
O 0 O71 


d the dyadics x, A, and v also have similar structure. 


PLASMA OSCILLATIONS (1) 


A particularly interesting situation tesults when, in 
the equilibrium state, € = (), and 


A=), Vp=0; WAG) = (0). hive JB =): (79) 
Eq. (78) then yields 
0. 
DS 3 
DD; 7 LD; (80) 
and 
Jz=6,D,, (81) 
where 
o*, + 0. 
= 2s, (82) 
Oy when w,t <1. (83) 
=O) 1203 Awhen @. 0 = 1. (84) 


In presenting the expressions for the components of 
the stress tensor, it is convenient to introduce the co- 
efficient of viscosity, in the absence of a magnetic field, 
©) 


| M \*l2 
where tT, = => (=| T, (85) 
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and the symmetric traceless tensor 
il 1 
A=> VU, + (VU,)"}—FV- Uy. (86) 


Then, in a-co-ordinate system in which B is along the 
z-axis, we have 
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where w, denotes the ion-cyclotron frequency. 


2.1 Weak Magnetic Fields 


In the limiting case of a weak magnetic field On te 
and it is to be observed from the foregoing relations 
that in this limit, to the lowest significant order, 


01 ~ 02, G,= 05 ts — 0", (88) 
De has Az ~0, 
Expressions (76), (77), and (87) now take the simple 


form 
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and 


Q, 
P= ea Ube ah Up + VV Uo)*| : (91) 


On substituting Eqs. (89)—(91) into the equation of 
state (61), we obtain 
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where y = 5/3. From the Maxwell equation ov 
—)B/0¢ we obtain 
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Thus our complete set of hydromagnetic equations 
ane: 
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whete the charge and current densities are determined 


by the Maxwell equations. For an ideal, perfectly con- 
ducting fluid, 01, > ©9, 


il 
E+ —U)xB 0,2 — 0,0 = 0; 


and if the characteristic frequencies of interest are much 
less than the ion cyclotron frequency and the foregoing 
relations reduce to the well-known hydromagnetic 
equations: 


0 
5; + V-(@U)=0, (08) 
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Vp+—J—eU) x B—e Ve, (99) 
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where in writing these equations we have dropped the 
subscript on Ug, as it is no longer necessary. 
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2.2 The A. C. Conductivity 


There is one further important transport co-efficient 
which has not yet been worked out for a plasma. Since 
one is often interested in frequencies of the electro- 
magnetic fields comparable to the effective collision 
frequency, the usual Chapman-Enskog theory needs 
some modification. This can be effected readily by intro- 
ducing terms in 0@ * /0/ on the left hand side of Eqs. (54) 
and (55). Now it is the electrons which primarily de- 
termine the conductivity. For the purpose of an order- 
of-magnitude estimate, we can neglect electron-electron 
collisions. Then, in the absence of temperature and | 
density gradients, magnetic fields, or mass flows, Eq. (55) | 
reduces to 
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where 
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If E varies in time like exp (—7 w/), then so must me | 
and the solution of Eq. (102) is 
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The associated current is then given by 


J=—e\®wwilo(w)P, (105) 


which, on using (O/m)"/: as the unit of velocity, can b 
written in the form 


J=cE, (106) 


where 
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If wt, <1, Eq. (107) yields to the lowest significant 
order 
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whereas if wt, > 1, we obtain 
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3. THE CASE OF WEAK COLLISIONS 


If the plasma is so dilute, or the temperature so high 
that the direct interparticle collisions of the first anc 
second class described earlier are negligible, then the 
particles interact mainly via the macroscopic electro 
magnetic fields. We may thus drop altogether the co 
lision term in the Boltzmann equation (3) in treating the 


dynamics of the plasma. Firstly, we wish to emphasize 
that this description of a plasma via the collisionless 
Boltzmann equation is entirely equivalent to the solution 
of the equations of motion for the particle orbits. This 
can be shown as follows. 


3.1 Equivalence of Collisionless Bolzmann 
Equation and Particle Orbit Theory 
Consider the collisionless Boltzmann equation or 


what is referred to as the kinetic equation in the Russian 
literature: 


of : 
ee yp yy =o. (110) 
The Lagrangian subsidiary equations are: 
dv 
er (111) 
ind 
dx 
Nimes (112) 
et the solutions to these equations be 
V == V (04, Xp: -- Xe, 2) (413) 
ind 
Ki KOs 5 Bayan? - Ng, 2) ; (114) 


vhere the «’s are the integration constants. We further 
uppose that Eqs. (113) and (114) can be solved for the 
’s. Now the general solution of Eq. (110) is an arbitrary 


unction of the constants of motion, ie., 
| 


TEV) = FG, Os Re); (115) 
or substituting Eq. (115) into (110) we obtain 
: of /0 
| dz <P kd UG Ps 1S Vo os] 
=1 
ayes 
fax ; " 
= sa a=? (116) 


nce the ~’s are constants along a particle trajectory. 
hus we have shown that the most general solution of 
q. (110) is an arbitrary function of the integrals of the 
agtangian subsidiary equations (111) and (112). But 
ese are just the equations of motion for the particle 
hits. The equivalence of the particle orbit theory and 
1e collisionless Boltzmann equation was first demon- 
rated by Jeans, and in the astronomical literature is 
ferred to as Jeans’ theorem. When the «’s are not true 
ynstants but adiabatic invariants, the resulting f satis- 
2s the Boltzmann equation to within the accuracy of 
ese invariants. The considerations of the orbit theory 
‘ovide a better insight into the physical factors involved 
a particular problem. However, the Boltzmann 
juation provides a rigorous treatment of complex pro- 
ems which may not be easily accessible to the orbit 


eory. 
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3.2 The Low Temperature Approximation 


Neglecting, then, the right hand side of Eq. (49), the 
equation for the stress tensor now becomes 


) 
wen Ot URE P WV UP y Dy 

é Shs a ; f | 
vale x P—P x B)—(J—eU) (E+ x B 


1 
— E+—UxB (J—eU)=0. (117) 
In order of magnitude |0P/o¢| ~ wv,,> and |V-Q\| 
~ kv ,* where v,, = (O/m)1!2, and w and & denote the 
characteristic frequency and reciprocal length respec- 
tively. Therefore, it is clear that if 
7) 


E> tn» 


Z (118) 


the heat flow tensor can be neglected in Eq. (117). For 
the linearized equations of motion w/k is the phase 
velocity of the wave, while for a non-linear system it 
might be termed the effective phase velocity. Thus if the 
phase velocity of the wave is much larger than the 
thermal speed of the ions or electrons, Kq. (117) te- 
duces to 


SP +V-(UP)4P-VUL@-VUy 


M1 6 


+78 X P—P x B)—(—eU) (E+, U xB) 


—(Bf=uxB\g—eu)=0. (119) 

Thus, by virtue of the approximation (118) we are 
able to obtain a closed set of moment equations. This 
scheme will be referred to as the Low Ti emperature Ap- 
proximation and will yield results correct to first signi- 
ficant order in the parameter of smallness Rv, |/@. 

Clearly, one can decouple the zth moment equation 
from its dependence on the (w + 1)th moment in an 
analogous fashion, and obtain the corresponding higher 
order thermal corrections. In all such calculations care 
must be taken to discard any sputious roots which may 
occur, that is, roots which are incompatible with the 
basic assumption (118). 


3.3 Chew, Goldberger and Low Approximation 


If we now assume that the magnetic field is so strong 
that the electron or ion cyclotron frequency is very large 
compared to the characteristic frequencies of interest, 
we can make a systematic expansion of the Boltzmann 
equation (110) in inverse powers of B. This has been 
carried out by Chew, Goldberger, and Low. An equi- 
valent procedure is to carry out this expansion in the 
equation of motion (119) for the stress tensor. If we 
assume that the Joule dissipation is negligible, i.e., 


[E+ Ux B|g—eU)| < Bx P|, (120) 


then the leading term in Eq. (119) gives 


BXP—PxB. (121) 
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It immediately follows from this equation that 


P=P 1, Pan, (122) 


where n denotes a unit vector along B. Eq. (122) iene 
plies that the pressure tensor is diagonal in a local rec- 
tangular co-ordinate system, one of whose axes points 
along n. In the plane perpendicular to n the pressure 
is a scalar of magnitude P,, and the pressure along B is 
P, which, in general, is not equal to P,. Equation (119) 
now teduces to 


d 
4 (Pr4 (P,P) an]+ d+ 
(Pi —Ppnol VU L[P,15 (Pi, 


P,)nnl-VU 


+ {P, I+ (2,—Pp)nnj:-VU=9, (123) 
where we have now set 
it 

Bo UX B%, (124) 


Eq. (124) implies that the medium has very high elec- 
trical conductivity so that it is unable to support, in its 
rest frame, any significant electric fields. Taking the 
trace of Eq. (123), we obtain 


d ih ibe oe 
Bl Prp 5?) lg a Ut YC 
+ (&.—P)an:VU—0. (125) 
The nn component of Eq. (123) is 
d 
Fe ie oh een ee (126) 
From Eqs. (125) and (126) we readily obtain 
d 
Sime oF ee aa DY UO) ell? (127) 
Using the relation (124) in the Maxwell equation 
cV <x E = —0B/0é, we obtain 
0B 
iio= Vie(Ur B): (128) 
It is convenient to write Eq. (128) as 
dB 0 
2 (2+ U-v)B=B-¥)U—B(V- 0). (129) 
From Eq. (129) and the equation of continuity 
d 
2 +4+0V-U=0, (130) 
it can be shown readily that 
Vue 
narVU=37,\ | (131) 


where B denotes the magnitude of B. On using this 
result in Eqs. (126) and (127), we obtain the two 
adiabatic relations 


d | PB 
5 (—-} 0 (132) 
and 
d | P, 
ii\pB)=° (133) 
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Starting from the Boltzmann equation (110), Chew, 
Goldberger, and Low have derived the same adiabatic 
relations on the assumption that the heat transport along 
the lines of force is negligible. | 

Equation (133) implies that, if one moves along with 
an element of the fluid, the average transverse kinetic 
energy per particle is proportional to B. This is com- 
patible with the well known constancy of the magnetic 
moment of a charged particle in the adiabatic approxima- 
tion. It can be shown that Eq. (132) implies the con- 
stancy of the longitudinal invariant of Chew, Gold- 
berger, and Low. Thus for a low density plasma in the: 
presence of a strong magnetic field, in the framework 
of the low temperature approximation of on the as-: 
sumption of negligible heat transport along the lines: 
of force and when the conductivity is very high, we: 
obtain the following complete system of equations: 


00 , 
sp eo, (135)) 
dU 0 
OT o(240-V|U 
1 
Vera @ eU) xX B—eVe, (136) 
oB 
a Voc (UeGn): (137)) 
P= P.T+-(P,—P) 04, (138) 
d | P, B® 
ee |=0, (139) 
and 
i Pe 
5 53 =o: (140) 


where the charge and cutrent densities are determine 
by the Maxwell equations. 


Ill. Wave Phenomena in Plasmas 


1. GENERAL REMARKS 


We shall now discuss the behaviour of plasmas foo 
states initially close to a steady state. This implies the 
existence of a time-independent, self-consistent solutiot 
of the appropriate equations of motion, jointly with 
the Maxwell equations, subject to proper boundary con 
ditions. The notion of closeness implies the existence o 
some dimensionless parameter of smallness, in terms 0 
which one can expand the time dependent quantities o 
the equations governing the system. In most situatio 
the appropriate parameter is the ratio of the macr 
scopic velocity to the phase velocity. The lowest ord 
system of the expanded equations is just that governing 
the steady state. To next order one obtains a set of line 
equations, whose coefficients do not in general depene 
on time, which describe the small motions about th 
steady state. 

If there exists no small initial perturbation who 
amplitude grows monotonically or sinusoidally in time 
then the system is said to be stable; otherwise, unstabli 
Ab at the absence of any dissipative mechanisms the sma. 
motions of a stable system may tend to persist i 


| mechanism they may tend to die out, the system te- 
turning to the steady state. 

If the system is unstable, after a certain interval of 

time the amplitude of the motion will become so large 
that the linearized equations will cease to be valid. The 
subsequent behaviour of the system will be governed by 
non-linear equations, whose analysis, in general, is dif- 
| ficult to carry out. 
_ There are two plausible situations into which these 
instabilities might possibly carry the system. The first 
| of these might be termed a turbulent State, characterized 
»by randomly fluctuating electromagnetic fields and 
| motions, which could manifest itself by an enhancement 
) of diffusion and resistivity over and above that resulting 
jfrom collisions. The second situation is associated with 
jinstabilities whose exponentiation time is short, can 
affect significantly the heavy particles and may lead to 
coherent macroscopic motion of the plasma, resulting 
in a radically different final state. The former of these 
has been conjectured as a cause of the loss of contain- 
ment in the stellarator type of controlled fusion devices. 
j}An example of the latter has been observed in pinch 
effect experiments where the discharge has destroyed 
itself upon achievement of a critical current. 

Since the coefficients of the linearized equations do 
not depend explicitly on time, one can seek a complete 


set of characteristic solutions whose time dependence is 
of the form 


| 
| 
| 
tte whereas in the presence of a dissipative 


wre 


exp (—/@/). 


(1) 


his leads to a characteristic value problem for w which, 
in general, may be complex. If there exists an w for 


hich 


Imw>0, (2) 
hen the system is unstable. If the coefficients also 
appen to be independent of one or more of the space 
atiables, one can seek solutions of the form exp (¢kx). 
he solution of the characteristic value problem, sub- 
ect to the appropriate boundary conditions, then leads 
© a relation between w and & which may be written 
mplicitly in the form 


F(w,k)=0. (3) 


telation of the form (3) is usually teferred to as the 
aspersion relation. The solution of Eq. (3) for w in terms 
f & would immediately disclose whether the system is 
table or unstable. 

Sometimes it is convenient to characterize the family 
f steady states by a certain parameter, say 8. As one 
hanges / continuously, the system may pass from 
ability to instability for some determinate value of #, 
y PB =f; 1.e., w passes from the real axis into the 
mplex w plane (with Imw > 0) at 6 = f,. The state, 
aracterized by 6,, which marks the boundary between 
e stable and the unstable regions is referred to as the 
arginal state. eT es 
In the case of static equilibria for non-dissipative 
stems governed by the hydromagnetic equations, it 
s been shown by Bernstein, Frieman, Kruskal, and 
ulsrud that the linearized equations of motion are self- 


> 
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adjoint; consequently the characteristic val 


ues ~" ate 
real. Therefore, 


OF 30) (4a) 
corresponds to stability and 
ate << (()) : (4b) 


to instability. However, for more general types of 
equilibria, the characteri 


stic values w may be complex. 
In case the imaginary 


part of w is positive, the system 
will be unstable; following Eddington, this type of in- 
stability is referred to as overstability. 

In addition to questions of stability which are usually 
posed for isolated systems, there is the important class 
of driven problems. In these situations a system in a 
steady state is excited by an external agency, perhaps at 
a fixed frequency. When the amplitude of excitation is 
small, the problem can again be treated by the same set 
of linear equations as those employed for the analysis of 
stability, but with appropriate modification of the 
boundary conditions. In the cases where a Fourier 
analysis in space is in order, one is again led to a dis- 
persion relation, but now it is to be viewed as determin- 
ing the wave number & (in general complex) as a func- 
tion of the prescribed frequency w, assumed teal. The 
fact that & can be complex permits waves to grow or die 
out in space, giving rise to the phenomena of amplifi- 
cation and evanescence. Problems of reflection, trans- 
mission, scattering, and amplification of waves by plas- 
mas fall into this category. It must be emphasized that 
when any particular problem is properly posed with due 
tegard to boundary conditions, there will be no con- 


fusion about the interpretation of the resultant dispersion 
telation. 


2. UNBOUNDED PLASMAS IN THE ABSENCE 
OF EXTERNAL FIELDS 


2.1 Longitudinal Electron Plasma Oscillation 


We first consider the simplest case of electron oscil- 
lations in a uniform plasma, neglecting the thermal 
motion of the electrons. The ions will be assumed to 
be at rest and to form the uniform neutralizing back- 
ground for the electron charge in equilibrium. Let NV 
denote the electron density (which is equal to the ion 
density) in equilibrium. In the perturbed state let the 
density be denoted by N +, where n is considered 
small so that all terms quadratic or higher in ” can be 
neglected. The equations of motion now ate (cf. Egs. 


(11-44) and (II-45)): 


S24 NV v=, (5) 
and 
Ov 
ee ee 6 
aa = eH, (6) 


where —e and m denote the charge and mass of the 
electron respectively. For longitudinal oscillations the 
electric field is determined simply by the Poisson equation 


V-E=4me=—42zen. (7) 


is 
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Eliminating v and E from Eqs. (5)—(7) we obtain for 
the equation of motion 


07 7 4 op Ne? n ; 
ro —= — a SSS ESS SS (ye WH 3 (8) 


where 


UL 


oe N e2\"ls 
Ory == | = 


It follows from Eq. (8) that the density fluctuations 
oscillate sinusoidally with the characteristic frequency 
w@,- The frequency @, is referred to as the plasma 
frequency. Formula (9) was first obtained by Tonks and 
Langmuir in 1929. 

The absence of space co-ordinates in Bq. (8) shows 
that these waves ate non-dispersive —their group veloc- 
ity vanishes and there is no tendency for a wave 
packet of this type to propagate through the plasma. 
‘As a result we can specify the phases of the electron 
displacement in such a way as to obtain a traveling wave 
a wave, however, which moves continuously through 
a fixed region without ever progressing beyond. 


2.2 Transverse Electron Plasma Oscillations 


We shall now consider the transverse (or electro- 
magnetic) plasma oscillations, neglecting the thermal 
motions of the electrons. The equations governing the 
problem are the linearized equations of motion (cf. 
Eqs. (1-44) and (II-45)) and the Maxwell equations 


0 
ST tNV-v=0, (10) 
mo = cB, (11) 
4a Ne 10E 
Vb == Soa (2) 
and 
10B 
From Eqs. (10)—(13) we readily obtain 
er 
Sie Or Ee NG E. (14) 


Taking E to be of the form expi(—w/+k-r), 
Eq. (14) leads to 


(o?—o,2—24)k x E=0, (15) 
(o? — w,2) (k-E) =0. (16) 


For transverse oscillations k x E does not vanish and 
Eq. (15) gives the dispersion relation 


and 


=O. a 0k: 


(17) 
while for longitudinal oscillations k - E does not vanish 
and Eq. (16) gives the dispersion relation already dis- 
cussed. In contrast to the longitudinal oscillations, the 
transverse oscillations propagate in the plasma with the 
group velocity. 


LY 


xe) ck 
= Fl + @e/@Ryh 


(18) 


ME eas 


i. 


-2% x 8980 N'sec?. (9) 


2.3 Effect of Thermal Motions on Electron 
Plasma Oscillations 


We shall now treat electron plasma oscillations, taking 
into account the thermal motions of the electrons using 
the low temperature approximation developed in 
Section II. In the equilibrium state, we shall assume the 
plasma pressure to be a scalar; however, we do allow 
for anisotropy of the material stress tensor during the 
perturbations. The perturbed state is characterized by 


v, N+, pl+Pp (19) 


where p = NO denotes the electron pressure in equili- 
brium, and p,, the perturbation in the stress tensor. The 
linearized equations of motion now are (ciaqs. (4) 
to (II-)46): 


ay ¢=0% (20) 
Ov 
LI sre) pi—NeE, (21) 
and 
m) 
2 p= —pW-I—plVvt Wy. 22) 


Equations (20)—(22) must be solved in conjunction 
with the linearized Maxwell equations 


4a Ne 10E 
and 
10B 
Via ae (24) 


Introducing &, the Lagrangian displacement of an ele- 
ment of the electron fluid from its equilibrium position, 
to a first approximation, we can write 


0§ 
i (25) 
Eqs. (20) and (22) can then be readily integrated to give: 
n=—NV-& (26), 
and 
Pi Ve ea (27) 


where the constants of integration have, by definition,. 
been set equal to zero. Assuming the space time behaviout' 
of all quantities to be of the form exp7(—ow/+k-r), 
the foregoing equations yield 


) 0 ., e 
and 
ekxkxE=+42N@0?§—a’?E, (293 


From Eqs. (28) and (29) we obtain for longitudinas 
oscillations the dispersion relation 


: ) 
or =a, +37 KR, (30! 


while for transverse oscillations we obtain 


(2k) (oF) + oto? =0. GL 


Eq. (31) leads to a quadratic in wm? whose roots ate 
given by 


2 z 2 2 22 ke Cc 1 ‘ Re — 
fp ARERR Ofm | *h ; 
a] i [w,? ao c2 p2 aE (R2 O|m))?{ < (32) 


In our approximation @? “> £2 O/m. Retaining terms to 
the lowest order in &? O/m, we obtain for the two trans- 
verse modes of oscillation (distinguished by the sub- 
scripts 1 and 2): 


oO; =o," ct Re ar or, 


(33) 


and 


eo Z| + (@z7/e* k2) 
The root @,? corresponds to the transverse oscillation; 
however, the root @,” 1s incompatible with our basic 
approximation that w?> &? @/m and therefore must be 
discarded. 
It may be remarked here that if, instead of using the 
equation for the stress tensor, we had used the classical 
adiabatic relation 


(34) 


d 
epi G lee Al aly (35) 
we would have obtained for longitudinal oscillations the 
dispersion relation 
2 een nee 

wt = wey 7 (36) 
instead of the relation (30). However, in the absence of 
collisions, it is necessary to allow for anisotropy in the 
pressure during the oscillations. 

The effect of thermal motions on electron plasma 
oscillations was first considered by Thomson and later, 
by Bailey. They obtained the dispersion formula 

0 uy 
wo? =o, + a (37) 
for longitudinal oscillations. It can be seen readily that 
the appearance of the factor O4?/m instead of 3 OK?/m 
ot y O?/m results from neglecting the change in tem- 
perature of the gas resulting from the density fluctua- 
tions. This amounts to considering the oscillations to 
be isothermal and not adiabatic. This is incorrect for the 
high frequency electron plasma oscillations; for a con- 
sistent approximation it is necessary to use the complete 
set of equations (20)—(22) to treat the problem. 

In part II of this survey this problem will be discussed 
from the standpoint of kinetic theory, and it will be 
seen that expressions (30) and (33) determine correctly 
to first order in © the real part of m. The higher order 
thermal corrections can be recovered by going to higher 
order moment equations. However, an essential con- 
sequence of the treatment from the kinetic equation in 
the absence of collisions, which is not recovered on 
using the truncated set of moment equations devoid of 
collision terms, is the phenomena of Landau damping. 
That is, the kinetic equation yields for isotropic equilibria 
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a negative imaginary part effectively proportional to the 
number of particles moving with the phase velocity of 
the wave. This number, of course, tends to be ex- 
ponentially small since, by assumption, the thermal 
spread is much less than the phase velocity in question. 
It must be emphasized here that the longitudinal and 
transverse oscillations are strictly uncoupled only in the 
case of a non-relativistic plasma and in the absence of 
any external magnetic fields, temperature or density 
gradients. The presence of an external magnetic field or 
inhomogeneities in plasma density and/or temperature 
result in the coupling of longitudinal and transverse 
modes, and the behaviour of the plasma, in general, is 
quite complex. This will be discussed in detail later. 


2.4 Ton Oscillations 


So far we have considered only electron plasma oscil- 
lations on the assumption that these oscillations are too 
rapid for the heavy ions to follow, which implies that the 
ions, therefore, may be considered at rest. Another class 
of oscillations, which is possible in a collision-free 
plasma, is the so-called ion oscillations. These are so 
slow that electrons see them as quasi static and con- 
sequently are distributed according to the Boltzmann 
distribution. During the oscillations it is appropriate to 
treat the ions and electrons as having different tempera- 
tures. However, we shall restrict ourselves to the case 
of longitudinal oscillations only. 

The linearized equations of motion now are (cf. 


Ras. (11-44)—(II-46)): 


Ie ay ee (38) 
Ov 
MNS =—V-p,+ Nek, (39) 
) 
7 Pi=—pV-vI—p[Vv+( vy], (40) 
V-E=47e(;—42,), (41) 
and 
ye N (exp On = 1| (42) 


where #, and 7, denote the perturbation in the ion and 
electron equilibrium densities (assumed equal) respec- 
tively, O, is the electron temperature; p; — NO, where 
Q, is the ion temperature, p, is the perturbation in the 
ion material stress tensor, and M is the ion mass. 
denotes the electrostatic potential defined as E = — Vp. 
Following the normal mode procedure outlined in § 2.3 
of this section, we obtain for the dispersion relation 


if OF 
2 eee ee oeeeee 3 
oY = Ta (eA aor +3k (43) 
where w,,;, the ion plasma frequency, is given by 
‘4 ot IN e\ "Is 
On; = i (44) 


and we recall that A,, the Debye shielding distance, is 
given by 


o ) (45) 


tho peri 


15 


I. B. BERNSTEIN &S.K. TREHAN 


Physically the reciprocal of A, acts like a screening radius 
for electrical forces in ionized media. Thus the potential 
in the neighborhood of a sphere charged with q units ot 
electricity is effectively given by 
22 & An? 46 
p=! exp (Aor) (46) 
in an ionized medium neat thermal equilibrium, com- 
pared to g/r ina dielectric. 


liven ie., the wavelength of the disturbance 
is large compared to the Debye length, formula (43) 
reduces to 


O, a L2 
i ee Ore 


C0” = Wy” B® Ay® + 3 A? (47) 
It is clear that if the ion temperature is very small 
compared to the electron temperature, the effect of the 
ion thermal motions can be neglected. We shall show 
in part II, using the kinetic description, that if QO, is not 
small compared with @,, the ion oscillations are strongly 
damped. It may also be observed from Eq. (44) that the 
ion plasma frequency is much smaller than the electron 
plasma frequency. 
The low frequency ion oscillations differ from the sound 
waves in the mechanism responsible for the organized 
motion, Ion oscillations are produced by long range 
Coulomb forces whereas sound waves are produced by 
short range collisions. 


3. UNBOUNDED COLD PLASMAS IN THE PRESENCE OF AN 
EXTERNAL MAGNETIC FIELD 


We shall now consider oscillations in a uniform 
plasma in the presence of a constant external magnetic 
field B,. We shall assume the plasma to be at zero tem- 
perature; this amounts to neglecting the thermal motions 
of the ions and electrons. The effect of thermal motions 
will be taken into account quite rigorously via the 
Boltzmann equation in part H. We shall further assume 
that each constituent of the plasma has in the equili- 
brium state a constant velocity u. The motion of the 
constituents of the plasma will, in general, give rise to 
a curtent distribution J. We shall, however, assume that 
the magnetic field which results from J is negligibly 
small compared with the external magnetic field. 

The linearized equations of motion for each of the 
species of the plasma are*: 


ov é Z 1 
sea ae 2 oes, AB Ae pero Bal ts) 


On 

ay Na a) (49) 
where the perturbations in the electric (E) and magnetic 
(B) fields are given by the Maxwell equations 


Amt 10E 
eS é Ue - On” 


* The notation used here (u for the velocity in the steady state 
and v for the perturbation in the velocity) is different from that 
used in section II. This should, however, cause no confusion. 


(50) 


16 


10B 
VES (1) 


DE 
and 


J=>e(Nv+2u). (52) 


The summation Eq. (52) is to be carried over both the 


_ species of the plasma. The cyclotron frequency is defined 


to be 
e By 


Mo 


= (53) 

We shall now assume that all quantities have space 

and time dependence of the form exp/ (—wt +k-r). 
Eqs. (48) and (49) then yield 
N 


= > {k- V) 


fae (54) 


and 


1 
_—i@—k-u)v=—|E+74 xB)—w,xv- 65) 


The appearance of the factor w — k-u instead of w is 
typical of a current-carrying plasma; the term k - u ac- 
counts just for the Doppler shift in the frequency w due 
to the motion of the electrons or the ions. From Eqs. (51) 
and (55) we obtain 


—i(w—k-u) (w,° V) = = (ee:E) |. 


~~ w.:[ux(k XE] 66) 
and 
—i(w—k-u) (w, X V) = = (we <x E)+ 
ow, x fu X (k X EI] — uw (ev) + oe. (57) 


MO 
Using the results (56) and (57) into Eq. (55), we obtain 
g i 


1 O2—(@—k -u)? 


i(—k-w E+ 54 (koe B)| 


+o.xEB+ hu, x {ux kx B} 


Dem \(w.-E) + © wo [u xX (k x BI} (58) 


From Eqs. (52) and (54) we obtain for the current density 
due to oscillations 


J=>eN(vta es). 


Equations (50), (51), and (59) lead to 
(w? — c? k®) E + c?k (k- E) 


(59) 


+4nioyeN(vtuc= ya) =0, (60) 


where v is given by Eq. (58). On substituting for v in» 
accordance with Eq. (58) into Eq. (60) we obtain, in 
general, an equation of the form 


RB = Q, (61) 


where R is a 3 x 3 Hermitian matrix. In order that 
Eq. (61) have a non-trivial solution, we must demand 
that the secular equation 


07 (62) 


be satisfied. This results in the required dispersion re- 
lation. We shall now discuss the implications of Eq. (62) 
in several special cases. 


3.1 Oscillations in a Static Plasma 


We shall first discuss the case of oscillations in a static 
plasma, ie., where the mean fluid velocity vanishes 
everywhere. Setting u = 0 in Eq. (58) and substituting 
the resulting expression into Eq. (60), we obtain an 
equation of the form (61). Taking k to be along the 
x-axis of a Cartesian system of co-ordinates and assuming 
that the magnetic field lies in the xz plane making an 
angle @ with the z-axis, the elements of the matrix R are 
given by 


= SA Oe ee 
2 +4 ih Wy" "Cos? 8 
il Dies @,2 — w 
5 WO We cos 6 


oo — fy = > io, Se 
12 21 PO — ow 


2 
Oo" 
9 Cc 

By=+Ry=— > o 2 


—, sinf cos 
O,2 — 


wn” 
2 Or! 
R — 2 — c2 h2— a5) 2 m,2 : : 
NEM Pen > 2 — oe 
3 WO ¢ sin f) 
Rog = — Ry = 7 . OF mee Nae. 


, (63) 


OW,” — w 
a — 9 > 2 
NX i w,2 sin? 6 
5 9 9 Cc 
R33 = 0” LEOne To Oe gy 8 2 
Lach pa CO —— 


where w,” — 42 Ne?/m. 


3.11 No Magnetic Field Case 


As the simplest illustrative example of Eq. (62), con- 
sider the case of electron oscillations. Setting w, = 0 in 
Eq. (63), we obtain Eq. (64), see below on this page, 
where w, now denotes the electron plasma frequency. 
We thus obtain for longitudinal oscillations (k along E) 


"= 0,7, (65) 


and for transverse oscillations (k perpendicular to E) 
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3.12 Propagation along the Magnetic Field 


For propagation along the magnetic field, 0 — 0, 
and Eq. (61) takes the particularly simple form (cf. Rqs. 
(61) and (63)) of Eg. (67), see bottom of this page. 


For longitudinal oscillations we thus obtain the relation 


=n 


I 


Da 2 LNs 2 
Dy = Wye [1 +- aI OO mans 


OS (68) 
where w,, denotes the electron plasma frequency. The 


transverse oscillations are circularly polarized and are 
described by the equation 


lw’ == c C= > One ms 2, IE — OF (69) 
where 
ye == E =e 1B > (70) 


the plus or minus sign corresponding to the right or left 
handed circularly polarized waves respectively. The 
oscillations are, therefore, governed by the dispersion 
relation 


aw? — ¢2 kK? _q, 2 


Oa ON, ae) 
for the right handed circularly polarized waves. Here 
the electron and the ion plasma frequencies are denoted 
by w,, and w,,; respectively, while the electron and ion 
cyclotron frequencies are denoted by w, and w, respec- 
tively. A relation similar to (71) applies to the left 
handed circularly polarized waves. 


a. High Frequency Oscillations 


If m > w,, the dispersion relation (71) reduces to 


7) 
2 a 
ow? —¢2 k? —_q,, ko Oy — OF 
In the limit when w > w,, Eq. (72) yields, correct to 
the lowest significant order in w,/@, 


(72) 


1 


We 
~ warner 


OO = 0? £2 Ew, 1 (73) 


Co" = Og? OP RE, (66) It is clear that a necessary condition for the validity of 
5 5 y, © 
which are in accordance with Eqs. (8) and (17). this relation is that w,2 < w,,” + c7k?. 
| wo — wy” — 6? k? 0 0 malig 
0 wo” — ow, —c? kK? 0 EY '=0 (64) 
0) 0 o—o," || £, 
Oy” ; ow 
oo ee4 > 2 woe tis oe — ot 0 F 
a” a 
‘amie » Oy” a 5 Da » Oy” 02 —o 0 y{=9. (67) 
0 0 oo = oy Zz 
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b. Oscillations neat the Electron Cyclotron Frequency 


Another case of interest is the one when the frequency 


of oscillations is close to the electron cyclotron frequency. 
In this case we find that Eq. (72) leads to 


oe (74) 


w= —w,|1—--; Se Rill c 
t De 2 (c2 bo One?) 


It is clear that the necessary condition for this approxi- 
mation to be valid is that 


Ong << | @,2 —C? &? |. 75) 
We can write Eq. (72) in the form (for @ > ,) 
j 1+ (eo) 
¢ ul ie (@,/@) =, (@ye"/ a)’ 
whete 17 = w/& is the phase velocity of the wave. The 
behaviour of the phase velocity as a function of the fre- 


quency of oscillation is illustrated in Figs. 1—3 for 
some values of the field strength. 


(76) 
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Fig.1 The behavior of Z = (V/c)?, where V is the phase velocity 
of the wave, as a function of x = w/w, for electron oscillations for 
propagation along the magnetic field for the case when w, = 0.1M,. 
The negative values of Z correspond to imaginary values of &, 
which correspond to absorption of waves in a bounded plasma. 


Fig.2 The behavior of Z = (V/c)? as a function of x = w/w, 
for electron oscillations for propagation along the magnetic field 
for the case when w, = Wy. 
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Fig. 3 The behavior of Z = (1 ‘/c)? as a function of x = «|, for 
electron oscillations for propagation along the magnetic field for 
the case when w = 10 @,. 


c. Oscillations much below the Ion 
Cyclotron Frequency 
Now consider the low frequency oscillations such that 


the condition wm <a, is satisfied. Then to the lowest 
significant order we can write Eq. (71) as 


(63) @ 
(1 +2) Og. == Ue 


O,; W, 


a) 
ot —c2®k2 + w,7 x 


(77) 


Using the fact that w,,?/w, = @,,?/@, and Oy, (O07 = 
c24a,N M/B,2, Eq. (77) leads to 


eae 
2 
oT ae ey (78) 
where 
Bo 
A= GaN My" (7) 


denotes the Alfvén speed. Eq. (78) is the dispersion 
relation of Astrém for the extra-ordinary hydromagnetic 
wave. It must be emphasized here that the hydro- 
magnetic wave, in principle, is just a special case, in the 
appropriate frequency region, of the well-known trans- 
verse electromagnetic waves. 

In the limit when c/A > 1, Eq. (78) reduces to 


ties 


i (80) 


that is, the phase velocity of the wave is just the Alfvén 
speed. In the other limit when c/A < 1, we get 


wt = 6222, 


(80a) 


which corresponds to the usual electromagnetic modes. 


d. Oscillations near the Ion Cyclotron Frequency 


Let us now consider the case when wm ~@,. This 


implies that w <w,, and the dispersion relation (71) 
reduces to 
Gots gh Re a2 ee 


- »* W—O; 


(81) 


This leads to the frequency of oscillation 


Wi” 
aad, h | o,? — (c? asia alk 


The condition for the validity of this approximation is 
that 


(81a) 


OF > w,2+ ck, 


3.13 Propagation Transverse to the Magnetic Field 


For propagation perpendicular to the magnetic field, 
6 = 2/2, and we obtain the secular Eq. (82) see below. 


One of the modes of oscillation corresponds to the usual 
transverse oscillation with electric field along B, and 
satisfies the dispersion relation 


wo = R41 YD wy? ~ c2@h2 1 Oe (83) 


The other two modes of oscillation are determined as 
toots of the equation (84), see bottom of this page. 


Consider the case ee aw < ck. Then to the lowest 
significant order, the 22 element of this determinant 
must vanish. We thus obtain 


2 2 
ee ee Dog 6 
' oF — a ' w,2— w 


(85) 


In the limit when w, < w <w, < w,, Eq. (85) reduces to 


2 2 
Orn) 
pr _ pe = QO. 86 
ot * w2 (86) 
The frequency of oscillation is then determined by 
oy 2 ay 2s 
We" Wy;"\"l2 : 
— | at) = (@, @,)" (87) 
pe 


This frequency of oscillation is referred to as the /ower 


hybrid frequency. 

Now consider the limiting case where the conditions 
wo? > w,? and w2>w,,7 are satisfied. Eq. (85) then 
leads to 


Qu (w,? a Dye”) /2. (88) 


This frequency of oscillation is referred to as the upper 
hybrid frequency. 

For the high frequency oscillations, w > w,, w,; 
we may neglect the motion of the ions, and Eq. (84) 
leads to 
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The behaviour of the phase velocity as a function of the 
frequency corresponding to the roots of this equation 


are illustrated in Figs. 4—6 for some values of the 
field strength. 
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Fig.4 The behavior of Z = (V/c)?, where V denotes the phase 
velocity of the wave, as a function of x» = = w/o, for electron 


oscillations for propagation transverse to the magnetic field for the 
case when w, = 0.1, 


Fig.5 The behavior of Z = (V/c)?, where V denotes the phase 
velocity of the wave as a function of x = w/w, for electron oscil- 
lations for propagation transverse to the magnetic field for the case 
when @, = @,, 


ol a 3.14 Weak Magnetic Field 
amass 2 2 LE ey , eak, Magn 
[uP 4-0, is | | + Wy ae eae i] | 
ao? w,? We now consider the effect of a weak magnetic field 
Vpsnigy 8 coe tgl (89) on the high frequency electron oscillations described 
é@ 
; wo ay W,* — c2k? 0 0 
5 oO,” ; 2 DW, 
2 ¢2 5; —f Pw 
0) uo ke ae we 22 —— we | ==) (82) 
| 2 
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Fig.6 The behavior of Z = (V/c)’, where V denotes the phase 
velocity of the wave as a function of x = w/w, for electron oscil- 
lations for propagation transverse to the magnetic field for the 
case when w, = 10 a,. 


in § 2 for arbitrary direction of the magnetic field. The 
weak field approximation, specifically, implies that 
w,< @. Retaining terms up to second order in @,/@, 
the dispersion equation now becomes (cf. Eqs. (61) 
and (63)) Eq. (90) see bottom of this page, where we 
have put 


Wy? = W,2 + 62k; (91) 
@, is the frequency of oscillation of the transverse wave 
in the absence of a magnetic field. Since we have assumed 
w,/@ to be small, to a fitst approximation, for the 
transverse oscillations, wm can be replaced by w, in all 
elements of the matrix R except Ay, and Ry. Then it is 
clear that R,, is the leading term in R, being of order 
c2k?2 and Eq. (90) gives 


CRE, =0. (92) 


Since & +0, we must require that /, vanish. Thus to 
first approximation in @,/w, the longitudinal and 
transverse modes of oscillation remain uncoupled. 
Retaining terms up to first order in w,/m, Eq. (90) 
reduces to 


ie 
w? — W,2 —1Wy¢" a COS On I as 
0 
ae =; (93) 
Wp —cos8 w*—w, ds) 
One an 0 y 


In order that this equation have a non-trivial solution, 
we must demand that the determinant of the co-efficients 
vanish; this leads to the dispersion relation 


wo 
a? =) + w,,2 — cos 0. 
aed ay 


(94) 


These modes of oscillation are circularly polarized with 


B,=+iB,,. It can be verified readily from Eqs. (90) 
and (94) that to first order in w,/w, E, vanishes. 

To carry through the calculations consistent to the 
second order in w,/m, we now substitute @ in accord- 
ance with Eq. (94) in all elements of the matrix R except 
R,, and Ry. in Eq. (90). Retaining terms up to second 


order in w,/m, we obtain 


mel Rs . 
— Oy” +. Gace Gy? (1 + cos* 0) + 


Dns Os 2 { One e : Open wo, 2 \ 
ee 4 cos? i aa Bae cos 6 +- Oot Wee cos? 8 { 
Wg\* f ees Ce 
4. (2 ie —cos 6)? + sin? @ 2 8 2 (1 + cos? 6) 
; a," | Is 
+ sin? 2a 1 (95) 


In particular, for 6 = 0, we obtain 


re: 9 Ve 2 Oe , L Ope” 
= Wry, oe Dive oa! = 2 we |, (96) 
and for 6 = 7/2, Eq. (95) gives 

i] Oe (ipa: Z 


We now find that the longitudinal and transverse 
waves ate coupled; their polarizations are given by 


1S . ss @ re @ 
BaF [1 Soe en os] (98) 
and 
Piet ak a et. 
EB =—7 2 ke Dye geet 6 (99) 


7] 0 


Let us now consider longitudinal oscillations which 
ate characterized by the plasma frequency w = w,, in the 
absence of a magnetic field. To a first approximation we 
can put wm = ,, in all elements of the matrix R except 
A, and obtain the secular equation (100) below. 
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—71Wy.@, CosO (100) 


ee) 
w”? —o,,2—w,? sin? 6 


The electrostatic oscillations are characterized by the 
condition that ck => wm. Thus the leading elements of the 
matrix are K,, and &,,, and to compute the corrections 
to the root thus obtained, it can be readily verified that 
it is sufficient to demand that the co-factor of R, 
This yields the dispersion relation 


, Vanish. 


/ 2 2 
OO, OO, 


w= w,,2 ji }. ls a soi sin? | (101) 
= pe 


The correction terms are small as long as w,2/w,,2 and 
@,”/c2k* are both much less than unity. 

‘Tt may be remarked here that for transverse oscillations 
the correction terms to m2, to lowest significant order, 
ate of order w,/m, while for longitudinal oscillations 
the correction terms, to the lowest significant order, are 
of order (w,/@)?. 


3.15 Strong Magnetic Field 


We shall now consider the case of a strong magnetic 
field which implies that m,/m > 1. Retaining terms up 
to second order in w,/m@, Eqs. (61) and (63) lead to (102), 
see bottom of page. In writing Eq. (102) we have used 


Tee oe + 
a © Fr M 
47 NM c? oe 
ae = se (103) 
and 
igen ot. (104) 
O, o, ; 
One of the roots of Eq. (102) is clearly given by 
2B 
i (105) 
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this is the dispersion formula of Astrém for the extra- 
otdinary hydromagnetic wave. The characteristic fre- 
quencies of the other two modes are determined by the 
algebraic equation (assuming ¢/A > 1) 

cA 2 


we (c? 2 + w,,”) + 02k 


ao rs @,," cos" 0 =). (106) 


The roots of this equation are given by 


2 


o = - {2 + ow,” 
= [c? al Dae 


In hydromagnetic waves ,, < ck, and Eq. (107) 
reduces to 


| i= 


this is the ordinary hydromagnetic wave of Alfvén. 


)2 4 c2 2 w,,? cos? OJ". (107) 


(108) 


PLASMA OSCILLATIONS @) 


4. INSTABILITY IN A PLASMA CARRYINe CURRENT 


We now consider the case of a uniform plasma com- 
posed of an equal number of electrons and protons in 
the absence of any external magnetic field. The ions will 
be assumed to be at rest, and the electrons, moving with 
a constant velocity 7. We shall neglect the magnetic field 
which results from the electron current, a procedure 
which is plausible if the associated magnetic forces are 
much less than the electric forces considered. For the 
problem at hand, Eqs. (58) and (60) reduce to 


(w® — ¢2 k2) E + 2k (k- E) 
+427 SeN |v Lu poe [= Om (109) 
fu Yo leo ‘ 
and 
Ze@ 1 
Nie raven ear \E+3 a xX (k x B)]. (110) 


We first consider electrostatic oscillations only and 
take k to be in the direction of the electron stream. 


Eqs. (109) and (110) then lead to the dispersion relation 
Ones 


2 
One 
=o § (w—ku)’ 


(111) 


where m,; and w,, denote the ion and electron plasma 
frequencies respectively. It is now convenient to intro- 
duce the dimensionless parameters 


a) ku Mm 
eres I~ Wn,’ and are (an) 
The dispersion relation (111) then takes the form 
i € 
TF (90) 9) = =p | ae ile (113) 


A schematic plot of F (x, y) as a function of x for 
y = 1is shown Fig. 7 for a hydrogen plasma. If the line 
F = 1 lies above the minimum of the curve, then it is 
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Fig. 7; The plot of F (x,y) given by Eq. (113) as a function of x 
Hoye jy) —= al “for a plasma composed of electrons. and protons. 
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clear that there exist four real roots of Eq. (113) and 
that all modes of oscillation are oscillatory. However, 
if the line / — 1 lies below the minimum of the curve, 
then there are only two real roots and there must exist 
a pait of complex conjugate roots, one of which has 
a positive imaginary part and is, therefore, unstable. 
The minimum of the curve F' (x, ) is determined by the 
condition that, for a given y, F (x, 9) be stationary, i.e., 


te) : — 
e) =2 =| ce) 
LOG io €> 


Fig.8 The behavior of y, given by Eq. (116) as a function of e, 
the ratio of the negative to the positive ion mass. 


) 
This leads to 
e's 6 
i eon 2D ife< ilies (115) 


The value of y for which this minimum is just tangent 
to the line /’ = 1 can be readily obtained by substituting 
this value of x into Eq. (113) and solving the resulting 
equation for y. Denoting these minimum values by x, 
and y,, we have 


Je = [ele + (1 — esp] (116) 


The behaviour of y, as a function of € is illustrated in 
Fig. 8. If y, >_y, Eq. (113) admits four teal roots. Thus 
a cold streaming plasma is sfab/e to all perturbations 
having wavelengths smaller than the critical wavelength 


A, given by 

u Dee \ee 
A= 2n— ae ! ze 
TG 


(117) 


ie 
a l1—ze), ife<1. (118) 


For values of y close to y,, by Taylotr’s expansion, 
Eq. (113) can be written as (cf. Eq. (114)) 


i 
(J Je) 1s (avs) = 2 (% — %e)” 1a es) =0. (119) 


This leads to 
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it is clear that this expression is valid so long as 


3 
Ve Sy ‘In. (121) 
Eq. (120) implies that 
dx 
SS SS CO J>y,. DD, 
Im ae Gol AS Voy), (122) 


For values of y close to the origin, and on the assumption 
that x/y <1, an approximation which is legitimate for 
the unstable roots (cf. Eq. (116)), we obtain 


(123) 


Equations (122) and (123) give the behaviour of x near x, 
and the origin as a function of y. For the continuous 
range of values of y between 0 and_y,, Eq. (113) has been 
solved numerically, and the results are illustrated in 
Fig. 9 where the real and imaginary parts of x are plotted 
as a function of ). 

The growth rate of maximum instability can be 
readily computed from Eq. (113) by observing that at 


the maximum 


xetiehy. 


dx 
a= = 0). 
dy 


From Eq. (111) we obtain 


(124) 


(125) 


Fig.9 The behavior of the real and imaginary part of the fre- 
quency of oscillation (x = w/w,) as a function of the wave number 
of the disturbance (&k = w,,y/u) for the unstable modes of oscillation. 


Since Im dx/dy vanishes at the point in question, it 


follows that the denominator occurring in Kq. (125) 
must be real; we take 


La ot Re at yee ye 


& |S 


(126) 


where « is a real positive number of order €/s. This 
choice of phase will be justified a posteriori. We use an 
asterisk to denote the value of any quantity cortespond- 


ng to the maximum growth rate of instability. Substi- 
uting Eq. (126) into Eq. (113) and separating real and 


maginary parts, we obtain 
it 
ye = | DS a a (127) 
ind 
11i+2.a , 
es: eae (128) 
l'o the lowest significant order, Eq. (128) gives 
WY} te 
— l= (129) 


Substituting this value of ~ into Eq. (127), we obtain 
o the lowest significant order 


(130) 


From Eqs. (126), (128), (129), and (130) we find that 


o E Ms if! E€ Is 3 3 "I 
Rete) last) Pelz) +. 


The maximum growth rate of instability, therefore cor- 
responds to (for a hydrogen plasma) 

* aU 1/ 

Im ~ = (5| ee (132) 

ne 
Therefore the smallest e-folding time is 18 w,,-1. For 
a plasma with a particle density of 10™ per cc, this time 
is of the order 10~! seconds. These instabilities grow 
so rapidly that it seems plausible to view them as provid- 
ing a mechanism for converting the organized energy, 
associated with the relative motion of ions and electrons, 
into the fluctuation energy of the oscillations. This 
provides an additional mechanism for the effective 
dissipation of the energy of organized plasma motion 
over and above collisions, and may consequently lead 
to enhanced resistivity of the plasma, compared to the 
corresponding value for a stable plasma. The behaviour 
of the real and imaginary parts of x* as a function of ¢ 
is illustrated in Fig. 10. 
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Fig, 10 ‘The dependence of the maximum growth rate of insta- 
bility as a function of ¢ is shown by the solid curve. The cotrespond- 
ing real part of the frequency is shown by the broken curve. 
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It may be emphasized here that these results are valid 
only so long as the phase velocities of the waves ate 
much less than ion and electron thermal speeds. When 
this problem is examined from the standpoint of kinetic 
theory in part II, it will be seen that the thermal correc- 
tions transform certain stable regions in the (#, «/A) 
plane into unstable regions, albeit with an exponentially 
small instability rate. Moreover, the region of the 
(v4, w/k) plane, corresponding roughly to the relative 
velocity u comparable to the ion thermal speed, is 
rendered stable. 

The case of two interpenetrating beams of electrons 
(or ions) can be deduced also from the analysis given 
here by setting € = 1 and by going over to a reference 
frame in which one of the beams is at rest. In this case 
one finds that the reciprocal of the smallest e-folding time 
is .866 w,, in contrast to (1/18) w,, as obtained earlier. 


5. PLASMA CARRYING CURRENT IN A MAGNETIC FIELD 


We now consider the case of a plasma, in which the 
electrons are moving with a constant velocity u relative 
to the ions, in the presence of an external magnetic field. 
We shall assume that the magnetic field which results 
from the steady state motion of the charges is negligible 
compared to the external constant magnetic field By. 
This is plausible if the fractional change in the magnetic 
field induced by these external currents is small over 
the wavelength of the perturbation in question. The 
equations for the problem at hand are (cf. Eqs. (58) and 


(60)): 
(w? — c222) E + ck (k-E) 


lee 
+4210 > eNv4 BoE a = (()- (133) 
and 
et 1 
~ mw O, a (wo —k : u)2 {o. (w,: E)— (wo —k-u)?E 


+ cae vy [ee (we) —k (kup (W-E) 


—i(mo~—k-u)w, X E—/w, x k(u-E)}. (134) 


We shall first restrict ourselves to the case when 
B,, u, and k are along the z-axis. Then w,  k vanishes 
and Eq. (134) reduces to 


5 Se ee | 


v~ mo oO =o keu) | 


— (w — ku)? (E 


— (@— kat fee 


ku)? 


— i (w—ku) ow, X Ej, (135) 


where n denotes a unit vector along the y-axis. From 
Eqs. (133) and (135) we obtain for longitudinal oscil- 


lations 
2 =>o pe a ae (136) 
? (w—ku)?| * 
and for transverse oscillations 
w—ku 
jor ae Dee wea | es = 0 (137) 
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where 


eae (138) 


Lif, 


represents the electric vector for the right and left 
handed circularly polarized waves respectively. Lhe 
upper sign in the denominator in Bq. (37) corresponds 
to £,, while the lower sign corresponds to Ee. 


When E. does not vanish, we obtain for longitudinal 


oscillations the dispersion relation 


iL 
ae 139 
> (wa — ku)? : i) 


It is to be observed from this expression that the lon- 
gitudinal oscillations ate unaffected by the presence of an 
external magnetic field (cf. §4). Consequently, our 
earlier remarks are applicable in this case as well. For 
propagation along any arbitrary direction the present 
discussion of longitudinal oscillations is also valid if one 
teplaces & by &, the component of k along the magnetic 
field. 

We shall now discuss the implications of the dispersion 
formula for the right handed circularly polarized wave 
only; a similar discussion applies to the left hand 
circularly polarized waves also and will not be given 
here. We have for the dispersion relation (cf. Eq. (137)): 


w—ku 


x 

2 2 22 2 

(2 == ¢2 pk? — a) =(0), 
SS > w—ku—o, 


Denoting the stream velocity for the ions and electrons 
by w,; and #, respectively, Eq. (140) takes the form 


(140) 
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where, it may be recalled, mw, and w, denote the ion and 
electron cyclotron frequencies respectively. Introducing 
the non-dimensional parameters 


(141) 
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we can write Eq. (141) as 
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We shall now discuss the implications of the relation (143). 


Let us first consider the case when there is no stream. 
Eq. (143) then reduces to 


ap ee a Ge" 
ie y EX—X Sieg 
(144) 


The analysis of Eq. (144) is most readily effected by 
means of Fig. 11, where F’ (x) is plotted schematically 
as a function of x. Observe that for large x, F(x) >—x2, 
F (0) = «?, and there are simple poles at x — — z and 


F(x) =1+e—x 
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x = €x, Thus the curve F’ (x) consists of three branches 
as indicated in Fig. 11 and there are clearly four real 
roots; consequently, all modes of oscillation are 
oscillatory and are s/ab/e. 


Fig. 11 A schematic plot of F (x) given by Eq. (144) as a function 
of x. The singularities occur at the electron and ion cyclotron 
frequencies. The curve F (x) intersects the x-axis at four points 
corresponding to four real roots. 


5.1 Stream Instability (Transverse Oscillations ) 


We shall now consider the transverse oscillations for 
the case where we have two beams of fast electrons 
traversing the plasma. Suppose the fraction 1 — f of the 
electrons is moving with velocity — u,, and the fraction B 
is moving with velocity u,. The net current due to the 
electron streams is 


J=eN|[—(1—f) u, + Bug]. 


In order that the net current due to the electron streams 
vanish, we must demand that 


(1 — f) 4 = Bug. (146) 


The corresponding results for the problem of one beam. 
can be obtained immediately from these results by 
setting #7, and 6 equal to zero. For the beam of rela- 
tivistic electrons, the only modification in the dispersion. 
relation is to replace the electron mass m by 


(145) 
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where ”m, is the rest mass of the electrons. The dispersion: 
relation now becomes 
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where w,, now denotes the electron plasma frequency’ 
defined with respect to the electron rest mass, and, it will| 
be recalled, that k and u are both along By. Introducing ; 


the dimensionless variable given in Eq. (142), we can 
write Eq. (148) in the form 


) 3 fe 
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Firstly, we observe that if we have only one beam of 
electrons, Eq. (149) reduces to 
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It turns out that Eq. (150) has real roots, corresponding 
to oscillatory motions, for all values of the parameters 
of physical interest. 

We now consider the case of two beams of electrons. 
The first point which is to be observed is that the effect 
of adding another beam is to introduce a new pole as 
shown schematically in Fig. 12. As w, increases from its 
zero value, the additional root moves to the left in Pior IZ. 


(150) 


Fig.12 A schematic plot of F (x) given by Eq. (150) as a function 
of x. Note the additional branch of the curve introduced by the 
ptesence of interpenetrating beams of fast electrons. The curve F (x) 
intersects the x-axis five times corresponding to all real roots. 


When 4, is so large that the gap introduced in the curve 
of Fig. 12 straddles the x-axis, two real roots are lost 
and a pait of complex conjugate roots are acquired. 
One of these must have a positive imaginary part and 
cotresponds to instability. This is illustrated in Fig. 13. 

When instability does occur it has a simple physical 
interpretation; namely, the condition that the gap 
straddles the x-axis is that 


‘ x + tp + (z/72) ~ 0. 


This can be interpreted as saying that, in its rest frame, 
the electron beam sees the circularly polarized wave as 
having a frequency equal to the electron cyclotron 
frequency, and the in-phase rotating component of the 
electric field is resonant with the particle motion. Since 
instability occurs only when the gap straddles the x-axis, 
it is clear that for this to occur the new pole must lie 
between the electron and ion cyclotron frequencies, 
.e., between — x and €z. 


(151) 
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Fig. 13. A schematic plot of F (x) given by Eq. (150) as a function 
of x for the case where the interpenetrating beams have velocities 
corresponding to unstable modes of oscillation. Observe that F (x) 
now intersects the x-axis only three times, corresponding to three 
real roots and two complex roots. 


3.11 Oscillations much below the Ion Cyclotron Frequency 


We now consider the consequences on the hydro- 
magnetic waves given by Eq. (78) of the introduction 
of beams of fast electrons such that the small gap straddles 
the real axis near the unperturbed root (denoted by 
subscript 1): 
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where in writing Eq. (155) we have used the fact that 
(l= P)e—= P52 (156) 


Making use of Eq. (152), we can write Eq. (155) as 
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For the case under consideration x/ex < 1 and, further, 
since we ate interested in the case when x ~ x,, we 
may set 


So eo (158) 
Eq. (151) then reduces to (assuming x/% < 1): 
il a ) 
P(x) =— Ria 
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Eq. (159) then leads to the following quadratic equation 
TO Dg 
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The roots of this equation are given by 
J x g 
x= 2 fe 5 + Je 4 = 
1 g x 4 ry} 
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It is clear from this expression that the maximum 
instability rate occurs when the term in the braces 
vanishes; this is indeed the approximate criterion given 
in Eq. (151). The spread in velocity of the electron beam 
about that given in Eq. (151) for which there is instability 
is given by (cf. Eq. (163)): 
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5.12 Oscillations near the Ion Cyclotron Frequency 


Let us now consider the case when the physical 
patameters are such that in the absence of the beam the 
interesting root is given by Eq. (81) — namely, near 
the ion Crete frequency. For the case under con- 
sideration, x < zy, and we shall assume that Dp is ot 
order x. Denoting the root in question by the subscript 2, 
we have (cf. Eq. (81)) 
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Equation (149) can then be written in the form 
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The roots of this equation are given by 
x=-4 Al EX Ne X_ | | 
2 : ‘a P a 
1 2) Gay) Ata 


The maximum growth rate of -instability occurs when 


x ! dt = 6 
2 Pi 
which is again similar to the approximate criterion. 
stated in Eq. (151). The spread in velocity about that 
given by Eq. (151) for which instability does occur: 
is given by 
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The instabilities due to the presence of streams of: 
electrons (or ions) may be responsible for microwave: 


noise in the laboratory plasmas and the radio emission: 
from astronomical objects. 


(172) 


6. SuMMarRy oF Sections HI2 —III5 
We shall present here a summary of the dispersion: 
relations for the various modes of oscillation of a plasma.. 
a) Static Plasma in the Absence of a Magnetic Field 
(i) Longitudinal electron oscillations (w? > &? O/m)) 


— 
OW? = Wy? + 3 = ke (1739) 


(ii) Transverse electron oscillations (wm? > k? O/m)) 


(6) oO, 2 


pe 
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(iii) Ion oscillations (O; < O,) 


w= w,,2 + c2 k? + ke 
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b) Static Plasma in a Magnetic Field, k || B, 
(iv) High frequency oscillations w > w,,;, 
We? K Wye + 2K? 


(175) 
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(vy) Oscillations near electron cyclotron frequency 
(w,? > w,? + c2&?) 
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(vi) Oscillations much below the ion cyclotron 
frequency 
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(vii) Oscillations near the ion cyclotron frequency 
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>) Static Plasma in a Magnetic Field, k | B 
(viii) Lower hybrid frequency (w < ck), 
Da —O <= Wn <= Oe 
@ = (w, @,)'!2 (180) 
(ix) Upper hybrid frequency, o >, 0 > o,; 
@ = (w,? + w,,)'l2 (181) 


1) Stream Instabilities 
1. Longitudinal oscillations, in the presence of a 
stream of electrons with constant velocity u, are 
unstable to all perturbations (along u) having 
wavelengths greater than 


E || m7 ae 
Es -(1+59) - 


(182) 


Smallest e-folding time = 18 w,,-1. 
This instability is unaffected by the presence of 
a magnetic field along u. 

2. Transverse oscillations, in the presence of two 
streams of relativistic electrons, are unstable near 
the Alfvén frequency and the ion cyclotron 
frequency, k and u being assumed to be along B. 


7. HYDROMAGNETIC WAVES FROM CHEW, GOLDBERGER, 
AND Low THEORY 


We have shown earlier that, for a low density plasma 
n a strong magnetic field, the kinetic equations for the 
ylasma reduce to equations for a single fluid if we 
reglect the heat flow along the lines of force and/or 
ssume the low temperature approximation. We shall 
10ow discuss the propagation of hydromagnetic waves 
ising this formalism. 

Consider, then, a spatially homogeneous plasma in 
tatic equilibrium in a uniform magnetic field, which we 
hall assume to be along the z-axis of a Cartesian system 
yf co-ordinates, with the pressure tensor given by 


P=p,1+(pj—P,) 00, 


vhere n denotes a unit vector along the field lines. In the 
erturbed state, let the various quantities be denoted by 


= 6 0, Pi af 6p» P| ar 6 pi), Bee OB, (184) 
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The equations governing the perturbed quantities are 
obtained by linearizing Eqs. (I-135)—(II-140). These are 


Ov , 1 : 
0s, SN re (V x 6B) xB, (185) 
oF 5 6 
aye 0 = — 0 YF OAPD (186) 
0 os 
5708 == Vx B)), (187) 
, OB O00 
OP, =p, ae nel ge (188) 
and 
. SB 00 
5 2 a ze 3 
dp Ping + 38RD (189) 


where 05 — B- oB/B. Introducing the Lagrangian 
displacement §,v = 08/04, Eqs. (186) and (187) can 
be readily integrated to give 


6g=—oeV-& 
OB=-V x6 xB), (191) 


where the constant of integration, by definition, has 
been set equal to zero. It follows readily from Eq. (191) 
that 


(190 
and 


eee (192) 


where we use a subscript | to denote the component of 
a quantity in the plane transverse to B. Thus 


Vi= Von. (193) 


The component of a quantity along B will be denoted 
by the subscript ||. We shall assume that all the perturbed 
quantities behave like exp 7 (— wt + k-r). Using the 
foregoing results, Eqs. (188) and (189) reduce to 


q 0p, = —p, (2V,-&, + 7k) §))) (194) 
an 
OV pj = — pi (V1: + 374) &i) . (195) 
From the form of P, it follows that 
V-6P=V,0p,+2(p| —p,) 4 On 
+n (p\—p,) V-6n + 7k 6p))]. (196) 


The change 6n in the unit vector along the lines of force 
is obtained readily from Eq. (191): 


én=n-VéE,+[(n-V)nn]-&. (197) 
Using the relation (197) into (196), we obtain 
V-6P= Vp, + (p:—P))) RiP e, 
+ #& (i — 21) Vr + ai] (198) 
From Eq. (191) it follows readily that 
(V x 6B) x B= B(V,V,-&,— 7) - (199) 


Using the foregoing results the equation of motion 
(185) splits up into the two equations 
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and 


(ow? — 3p Ry?) €)) = Ap, (Ky §1)- (201) 
On substituting for 4) in accordance with Eq. (201) 
into (200) and forming the scalar product with k,, we 
obtain the dispersion relation 


: B?\ é 2 
ow + & |p Pi ead 2k,2( p+ go 
4 2 ke? 
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If 6 denotes the angle which the propagation vector 
makes with the y-axis, 
ki =kcosb, je, = asin’. (203) 


After some simplifications, Eq. (202) reduces to the 
following quadratic equation for w*: 


2 f Ba 
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The roots of this equation ate given by 


2 | B2 
ee | p, + 2p cos?6 + p, sin?) 


9 


C25 
ee (|B as he 
+ Do \ax + p,(1-+ sin? 0) —4p\cos 0) (205) 
Ia 
+ 4p,2sin®0 cos? 


For propagation along the lines of force, 6 = 0 and 
Eq. (205) leads to the following two modes of oscillation 
(distinguished by the subscripts 1 and 2): 


k2 | B2 
oP = maa Pie (206) 
and 
3k 
Wo? = 0 P\l- (207) 


It is clear that the mode of oscillation corresponding to 
Eq. (206) is unstable if 
2 
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That is, if the distribution function for the ions and/or 
electrons is strongly peaked along the lines of force, 
the hydromagnetic waves propagating along the field 
lines are unstable. 

For propagation transverse to the field lines, 0 = a/2 
and we obtain 


pee See? 
ot |e +2): (209) 
The corresponding phase velocity of the wave is given by 
(60) 2P "Ye 
Vp ib = =| ) (210) 
where P = p, + (67/82) denotes the total transverse 


pressure. 
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For any arbitrary direction of propagation there exist 
unstable modes of oscillation for all 9’s less than 6, as 
determined by the equation 


= + p, (1+ sin? 6,) = sin? @, + pjcos?6,. (211) 


As an example, if we take p, — B?/8a and p, = B?/2x, 
we find from Eq. (211) that 0, = 26° 45’. 

It will be shown in part II of this survey that in the 
presence of a magnetic field, whenever the distribution 
function 18 anisotropic’ G.e., Pj) << P, of Py P\\), we 
have instabilities, although their e-folding times tend to 
be rather small. 


8. HyDROMAGNETIC WAVES: GENERAL THEORY 


In our previous simplified discussion of hydro- 
magnetic waves we neglected the effects of the finite 
temperature, viscosity, electrical conductivity, thermal 
conductivity, and other transport co-efficients of the 
plasma. These shall now be taken into account in order 
to exhibit the coupling of the hydromagnetic and sound 
waves, and to discuss the damping of these waves due 
to dissipative effects. However, for the sake of sim- 
plicity, we shall assume that the ion gyration frequency 
is much smaller than the effective collision frequency 
so that all the transport co-efficients may be taken to be 
isotropic. Eqs. (II-94)—(I-97) then form the basis of 
out discussion. — 

In the equilibrium state U=0O and B= B,é,, 
where B, is constant and p and @ are uniform. For small 
departures from the equilibrium state let the various 
quantities be denoted by 


Vi, Papi» 04 Cis and B, + B,. (212) 


The equations governing the small departures are 
obtained by linearizing the equations of motion 
(II-94)— (11-97). These are: 
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where s = (yp/@)s denotes the sound speed and D, is 
obtained from the linearized form of Eq. (I-89). Thus 
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The temperature and the pressure are connected by the 
relation 0 — Mp/o; the linearized form of this gives 


M | Q 
O1= 0 (>. aly al (2183 


Introducing &, the Lagrangian displacement of an 
lement of the fluid, and assuming that all quantities 
ehave like exp 7(— wt +k-r), after some straight- 
orward reductions Eqs. (213), (215), and (216) lead to 
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Pi= PQ), (221) 
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t must be emphasized here that in writing the foregoing 
xpressions we have assumed the various transport 
o-efficients to be isotropic. However, this may not be 
fue, in general, for oscillating electromagnetic fields. 


Using the foregoing results in the equation of motion 
214), we obtain 
, ul I il u 
9 5 i 9 2 ° i 
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where n is a unit vector along B, and A is the usual 
Alfvén speed. Assuming that k lies in the (xz) plane 
making an angle 9 with the z-axis, Eq. (224) leads to the 
‘ollowing equations for the different components of &: 


(224) 
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8.1 Absence of dissipation 

Let us first consider the limiting case of an ideal 
ydromagnetic fluid for which all the transport co-effi- 
ients vanish. In this limiting situation « = 1, p= 1, 
nd uw = 0. Eq. (225) then leads to 


w= +Ak cos 6; (228) 


ese modes of oscillation remain unaffected by the finite 
lasma pressure. For Eqs. (226) and (227) to have a 
on-trivial solution, we must have (setting « =—1, 
= 1, i= 0): 

@w2 —A? k2 — 52 2 sin? 9 — 5? &2 sin cos 6 
s2 k2 sin cos 6 w2—s? k2 cos? 6 


=0. (229) 
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This leads to the following quadratic for @?: 
aot — @? (A? + 52) R2 1+ A? 52 £4 cos? = 0. 


The roots of this equation are given by 


ii : {a +4 52) f2 
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In the case of a strong magnetic field, s/ A < 1, and 


Eq. (231) leads to the following roots (distinguished by 
the subscripts 1 and 2): 
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and 
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In the case of a weak magnetic field, A/s <1, and 
Rq. (231) leads to 


@,? = R* (52 + A? sin? 0) (234) 


and 
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The principal result which stems from including the 
finite plasma pressure is that the hydromagnetic waves 
remain unaffected only if the motion of the material 
(i.e., §) is normal to the plane containing the magnetic 
field and the propagation vector; otherwise the disturb- 
ance breaks up into coupled hydromagnetic waves and 
sound waves. The chief application of this result is to 
waves ofiginating in the surface layers of a sunspot 
where A may be comparable to s. These remarks may 
be of some significance also in laboratory plasma 
diagnostics. 

That the @? are real in the case of ideal hydromagnetics 
also follows from the self adjointness of the equations 
of motion which may, in general, be written in the form 


—wg§ =F ©), (236) 


where F is a linear operator acting on §. The self adjoint- 
ness of F immediately leads to the conclusion that w? 
is real. Consequently, the criterion for stability is that 
there be no negative characteristic value wm*. Also, it can 
be shown easily that a necessary and sufficient condition 
for this circumstance is that the potential energy 

100, 
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always be positive. This latter condition is a statement 
of the energy principle governing the stability of ideal 
hydromagnetic system. This criterion, while derived on 
the assumption of collision-dominated plasma, can be 
demonstrated to give a lower bound on stability as 
compared with the conclusions of the theories based on 
the collisionless Boltzmann equation in the limit of a 
strong magnetic field. 


(236) 


8.2 Exffects of dissipation 


We now consider the effects of dissipation on the 
various modes of oscillation discussed in the preceding 
section. Eq. (225) gives 
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w — A? 2 cos? 0 (1-4 
If c2h2/4 700, and {iw/o A? are both much less than unity, 
to the lowest significant order Eq. (237) yields 
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We now consider Eqs. (226) and (227). In order that 
these admit a non-trivial solution, Eq. (239) see bottom 
of this page. 

On the assumption that ¢?k/420,0, MA?/Ow, and 
[A, — (1,%,/0,)] [k2.M/wo] are all much less than unity, 
we obtain in the strong field case (s/A < 1): 
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Under the same approximations for the weak field 
case we obtain 
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It is clear from Eqs. (240) and (241) that the damping 
of the hydromagnetic waves is 


REN fl ae 2 se 
lee agro, oe. oes \ 


yey 2 i all 


0, / @ y 
(244) 


Since s/A <1, we may write this approximately as 
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This suggests that only the resistivity and the viscosity 
of the medium are significant for the dissipation of 
hydromagnetic waves. In order that these waves be 
realized in any physical situation, it is essential that the 
time required for their amplitudes to be reduced by 
a factor of ¢ 1 of their initial amplitudes be much larger 
than one period of oscillation of the wave. Thus we 
must have 
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Since & = 22/A, this relation implies that the wavelength 
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of the disturbance must be greater than the critical 
wavelength A, given by 
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IV. Oscillations of Bounded and 
Non-Uniform Piasmas 


1. OSCILLATIONS OF A CYLINDRICAL PLASMA 


We shall now investigate the oscillations of a cylin- 
drical plasma immersed in a uniform external magnetic 
field which will be assumed to be along the z-axis. 
In the equilibrium state the plasma is at rest and the 
electron charge density is neutralized by the ion charge 
density. We shall assume that the plasma is so dilute 
that direct interparticle collisions are negligible and that 
particles interact through long range electromagnetic 
interactions only. For the sake of simplicity we may, 
therefore, neglect the plasma pressure. For small depar- 
tures from the state of static equilibrium the equations 
of motion for the electrons and ions, assumed singly 
charged, are 


Ov, il 
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together with Maxwell equations 
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where 
J=eN(v,—v)- (4) 


We have neglected the displacement current as we are 
interested only in low frequency oscillations satisfying 
the condition wm < ck. From the definitions of the mass 
density and the mean mass velocity, we obtain 
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Ley i 
== |B x Bo (5) 
and 


Ov il 
Os = 1x Bo: (6) 
neglecting terms of order w/M. We now need a relation: 
connecting v,, v, and J. From the definitions of J and v 
we readily obtain 
1 

~ Nel 
From Eqs. (5) and (7) we get 
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As we are interested in frequencies of oscillation much 
below the electron cyclotron frequency, we may drop 
the term on the left hand side of Kq. (8) and obtain for 
Ohm’s law: 


1 ; 1 
Be Yo Beane) X B= 0. (9) 
From Eqs. (3) and (9) we then obtain 


oB 2 1 
37 = VX (¥ x Bo— 3y,] X By). (10) 
Equations (6) and (10) form the basis of the charac- 
teristic value problem which must be solved subject to 
appropriate boundary conditions. We now assume that 
the #, 0, and x dependence of all quantities has the form 
exp 7 (— wt + m0 + kz). Equations (6) and (10) then 
lead to 


—3B=V xf x B) Xn] Xn+ ix (V x B) x a} 
(11) 


where A denotes the Alfvén velocity and n is a unit 
vector along By. On carrying out the normal mode 
analysis, we obtain the following equations for the 
various components of B: 


ee ea pine 8 a Se 
oS Aik| 22 eB. +k ga — 2)" |, 
(12) 
mk. | aw ow dB, 

L 0222 \.R = gece e 
B= A : ta Ss — 2A) BR kOS | (13) 
and 

Mee pd Boel re IP Nap 
aA et dE Del Ce 
where 
1 cw? \2 
Se A 4Q2 
4=(#—Se] —A@, (15) 
1 1 
= Ak — (w? | A?) RQ?’ (16) 


and (2 =w/w;. The solution of Eq. (14) which is 


regular on the axis of the cylinder is given by 
a aaa Bb, m (zr) > (17) 


where B, is an arbitrary multiplicative constant. Having 
determined B,, B,, and By, are readily obtained from 
Eqs. (12) and (13). This determines the perturbation 
in the magnetic field inside the plasma. 


1.1 The Perturbation in the Vacuum Field 

‘The vacuum magnetic field satisfies the equations 
Baw V0; (18) 

The solution of Eq. (18) which is bounded at infinity is 
given by 


and 


fe = BK, (kr), (19) 
where B, is another constant, K,,, (x) is the Bessel func- 


tion of purely imaginary argument (the notation is that 
of Watson), and we have suppressed the factor 


exp 7 (m0 + kz). 


PLASMA OSCILLATIONS (1 


1.2 [ke Boundary Conditions 


The boundary conditions which must be satisfied at 
the perturbed boundary of the plasma are 


AB 11) == 0), 20 
and ey 

(es 2 

Ale |= 0, (21) 


where A(X) denotes the jump experienced by the 
quantity at the interface. The linearized forms of 
these equations are 


BP=Be (=P (22) 
and 
Dae =e Be (r= 1 (23) 


where have now used the superscripts p and v to indicate 
the quantities referring to the plasma and vacuum 
respectively. Using these boundary conditions, and after 
some straightforward reductions, we obtain the dis- 
persion relation 


9 . we w 5 Fm (48) 
3 ) | L2 - 2k) 4 Caner 
RA | mw Q ri (é Aa Dk | ak m ERI 
= KR Ge (k k) 4 


For a specified & A, and, 2 Has. (15) and G6) can 
be solved to give wm in terms of vy. When the w thus 
obtained is substituted into Eq. (24), it admits an infinite 
number of discrete roots for y for prescribed values 
of k, A, and @; denote these roots by y,. The disper- 
sion relation then takes the form 


cw | Ww | | cw il 


DE Sez tyre ee A) 


1.21 Low frequency oscillations 


Consider first the case when w/w; < 1. Then Eq. (25) 
yields the roots 
On = RA (26) 
and 


(27) 


which correspoad to the ordinary and the extra-ordinary 
hydromagnetic waves. 


aie AP op 


1. 22 Oscillations near the ion cyclotron frequency 


We now consider the case when w/KA <« 1. Eq. (25) 
then yields 


9 2 Oo? wo? 
or = 0," |) 1 ~ R2A® AB (RB 2) | 


These oscillations are referred to as ion cyclotron oscilla- 
tiONs. 

In ion cyclotron oscillations the ion velocities are 
circularly polarized and the sense of rotation is the same 
as that for a free ion moving in a magnetic field. The 
electric field is elliptically polarized. One may decompose 
this electric field into two circularly polarized compo- 
nents, one of which will have the same sense of rotation 
as the ion motion in the wave. However, the ion velocity 
is 90° out of phase with that component of the electric 


(28) 
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SPECTRA FOR ION CYCLOTRON WAVES 
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growing spiral orbits and there is 
strong power absorption by the ions 
from the electric field. There are 
no waves at this very low density 
limit. 

For moderately high densities 
(Fig. 14) there are ion cyclotron 
waves. The different vertical lines 
indicate the different modes of os- 


t 
Bo 


B> 


WAVES } 


DAMPED WAVES | 


WEAK POWER ABSORPTION BY 
IONS FROM ELECTRIC FIELD 


cillation, and for the high density 
condition all of these modes occur 
at magnetic fields considerably stron- 
ger than the cyclotron field By. Once 
the ion cyclotron wave has been set 


NO WAVES 


STRONG POWER 
ABSORPTION BY 
IONS FROM THE 
ELECTRIC FIELD 


a HIGHER DENSITIES LOWER DENSITIES oe ee do a cages ae 
| LONGER AXIAL WAVE-LENGTHS SHORTER AXIAL WAVE-LENGTHS Coa ae eee a ae 
; there is no damping of the wave. 

The second picture in Pig. 14 
2uf = oo n= radial mode number shows the spectrum for ion cyclo- 


Fig. 14 


field which rotates in the same sense as the ion. Con- 
sequently, there is no absorption of energy by the ions 
from the waves. 

However, in the foregoing discussion we have not 
taken into account the thermal motion of the ions, 
and/or we have neglected also the collisions (short range) 
between the ions and electrons or the like particles. 
In either case damping of the ion cyclotron waves will 
result; in the former case we have Landau damping and 
in the latter case we have the usual collisional damping. 
When damping occurs there appears a component of the 
elliptically polarized electric field which is in phase with 
the ion motions. This leads to absorption of energy by 
the ions from the electric field of the wave. Conse- 
quently, this phenomenon is of fundamental interest 
in connection with the heating of the plasma. 

In the experimental situation one has the plasma in 
a magnetic field and one superposes on this a transverse 
electric field which oscillates at the ion cyclotron 
frequency. As the ion motion is in phase with the 
electric field, energy is fed into ion motion, and the 
ions will move in a spiral of increasing radius. 

At higher densities, say 101° ions/cm, this single 
particle picture is no longer valid. Ion currents induce 
an electric field which has a decelerating effect on the 
ions, and heating a moderately dense plasma at the ion 
cyclotron frequency becomes rather inefficient. However, 
if we keep the frequency of the oscillating electric field 
fixed, but go to slightly stronger magnetic field, we find 
the existence of natural oscillations of the plasma. These 
natural oscillations, which we call ‘‘ion cyclotron waves,” 
are, in fact, the short-wave, low-density limit for the 
extra-ordinary hydromagnetic waves described by 
Alfvén and Astrém. 

Fig. 14 shows the spectra for ion cyclotron waves. 
The resonances are indicated as lines and are plotted 
against magnetic field strength. We first look at the 
picture denoted by “no waves.’ Here we have low 
densities and, as a consequence, single particle accelera- 
tions of the ions in the electric field. The ions move in 


a2 


tron waves for somewhat lower den- 
sities. The various modes now occur at magnetic field 
values which are somewhat closer to By. There is still no 
damping of the waves. In the third picture the density is 
quite low and the modes now occur at magnet ic fields very 
close to B,. This is an intermediate case between the 
ion cyclotron wave resonance and the single particle 
resonance. For this condition we have waves, but the 
ions begin to move in growing spiral orbits and pick 
up energy from the electric field. In this case we have 
ion cyclotron waves which are slightly damped. 

The transition from the undamped waves to the slightly 
damped waves is particularly interesting because it 
provides a method for thermalizing ion cyclotron 
waves, that is, for transforming the wave energy quickly 
into heat. It turns out that it is equally valid to talk of 
wavelengths instead of densities, and it is the short 
wavelengths which are damped. To illustrate a plasma 
heating scheme, we make an analogy to ocean waves 
running up on a beach. A cross section of an ocean 
beach is shown in Fig. 15. The waves are moving in 
toward the shore and, because the water is getting more 
and more shallow, the wave length of the ocean waves 
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becomes shorter and shorter. Finally, in the shallow 
water the waves are unable to propagate and they 
“break;”? the wave energy is transformed into heat. 
Similarly, one might use an induction coil to generate 
rather long wavelength ion cyclotron waves and allow 
these waves to propagate along the magnetic lines of 
force through a region where the magnetic field tapers 
off and becomes somewhat weaker. The wavelength 
becomes shorter and shorter and, finally, when the 
wavelength is sufficiently short the waves damp out and 
the wave energy goes into heating the ions. It woud be 
in this “magnetic beach” region that the real heating of 
the gas would take place. 


2. RADIATION BY PLASMA OSCILLATION 


We now wish to exhibit the possibility of radiation 
by plasma oscillations. Firstly, it may be remarked here 
that this problem is mathematically quite analogous 
to the scattering of electromagnetic radiation by a 
dielectric. Consider, then, a simple equilibrium config- 
uration in which the ions are held fixed with a uniform 
density NV and the electrons, whose thermal motions 
will be neglected, are distributed likewise inside a 
sphere of radius RK. From the linearized equation of 


motion 
2 
als eE 


Mh y= == 
07 


(29) 


and the linearized Maxwell equations, we readily obtain 
eV XV x E=(w’—o,”) E, (30) 


where we have assumed the time dependence to be of 
the form exp (—/mf), and mw, denotes the electron 
plasma frequency. From Eq. (30) it follows that 


V -[(@? —@,”) E] = 0. (30a) 
Inside the sphere w, is constant; consequently, if 
w? +,” Eq. (30a) leads to 


Vie Bi2= 0 (31) 


for both inside and outside the plasma. 

Because of the symmetry of the problem, it is con- 
venient to seek a solution of Eq. (30) in terms of vector 
spherical harmonics: 


A= 6,1," (059). (32) 
Ws 
B= W+ Ds v Yi" (8, 9), (33) 
and 
Caan VX" OO], G64 


e+] 
where Y™, is the normalized spherical harmonic. A, B, 


and C constitute a complete set of basic orthogonal 
vectors on the surface of a sphere; namely, 


P47;9,9):-Q¢4%6,p4d2=0, 


where P and Q are any two distinct vectors A, B, and C. 
The normalization condition 1s 


S[P 4m”; 0,9) -P (CAZES 9,9) dQ= O11, Orn my" 


(35) 


(36) 
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The spherical harmonics satisfy the equation 


Sa, ] : - 
VY" 0,9) =—sM/+1¥"(0.9). G7) 


We now write 

F G H 

ES Sie BG: (38) 

where /', G, and H/ are functions of r alone. On substitut- 
ing Eq. (38) into Eq. (31) we obtain 

af ap 

be e+ 1)| G. (39) 
On substituting Eq. (38) into Eq. (30) and making use 
of the relation (35), we obtain the two independent 
equations 


ad 

a EG w,2) ve 1)] #0 (40) 

} ¢ } 

and 

dG if 1 1 

ie + eye [ales — allt 0 Fo. 
(41) 


The problem, then, is to solve Eqs. (39)—(41) subject 
to appropriate boundary conditions. 


2.1 The boundary conditions 


From Eq. (30a) it follows upon integration across 
the interface that we must have 


(wo? — w,") E,-n=@*E,-n, (42) 


where the subscripts 1 and 2 refer to the plasma and 
vacuum sides respectively and n denotes the unit outward 
normal to the surface. Since only A has a radial com- 
ponent, it follows from Eq. (42) that 


(w? — w,”) PF, = wo F, (43) 
It follows from cV x E = — 0B/0¢ that the tangen- 


tial components of E are continuous at the interface. 
This requires that 


Ae Pp == IR. 


GiSG, atr="h, (44) 
and 
baa, atnaok. (45) 
Integrating Eq. (40) across the interface, we get 
Gils Olas - 
fpr oy aM IP = IR (46) 
The conditions (45) and (46) may be combined to read 
d 
Ce = In Hy. (47) 


dr Ear 


2.2 Solution of the problem 
The solution of Eq. (40), which is regular at the center 
of the sphere, is given by 
Fj hier youd, (ey) (48) 
where 


(49) 
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and H, is an arbitrary multiplicative constant. It can be 
seen readily from Eqs. (39) and (41) that F satisfies the 
same equation as H; consequently, a similar solution 
holds for F. G is then determined by Eq. (39). 


2.3 Tle case of a perfectly conducting boundary 

An interesting situation which may approximate 
cavity resonators is the one where the plasma is enclosed 
by perfectly reflecting walls. The boundary conditions 
to be imposed in this case are that E, and E, vanish 
at r = R. This requires that 


Tivtl, (Ar B) = 9. 


Let the roots of Eq. (50) be denoted by /;,,,. From Eqs. 
(49) and (50) we than obtain the dispersion relation 


(50) 


2 


Ce 
wo =o," + RJ” (51) 


I,m. 
It may be remarked here that as K > co, mw > @,, the 
natural plasma frequency. 


2.4 The case of free boundary 


In the case of a free boundary, the inside solution 
must be matched to the one outside with due regard 
to the boundary conditions. The equations governing 
the vacuum field can be obtained readily from Eqs. (40) 
and (41) by setting w,? =0. Thus H/ satisfies the 
equation 


gd /7+1 

cat |e] zH=0, oe) 
where &, = w/c. The solution of Eq. (52) is given by 
H = rh [Hy HS), (her) + Ho Ais), er), (53) 


where H,and H,” are the two kinds of Hankel func- 
tions whose asymptotic forms are given by (Ayr > 1) 


. 2 4 tlx—Gl/4) 2 +1) 
pig PS eA / (54) 
and 
2 Nai FIR) (29 =F 1). 
Gees a ; (55) 


In Eq. (53) Hy and Hy, ate constants of integration. 
We have assumed the time dependence of the wave to 
be of the form exp —/w#. Therefore, in order that at 
large distances from the sphere we have only an outgoing 
wave, we must set (),=—0. Thus the permissible 
solution is 


feb ila beberle (56) 


Applying the boundary condition (47), we then obtain 


d 


d 
Fp yyy, rt) = gla SY), (her). 


(57) 


We now festrict ourselves to the special case where 
k,Rk > land &,k > 1. We may then use the asymptotic 
expansion 


Jel) ~{zz) 08 |[e—G erty] (58) 


on 


and Eg. (54) for H,'?. Eq. (57) then leads to 


— ky tan [Ay RZ + 1)| = Aho. (59) 
Using the relation 
ate cans 3; in =e, (60) 
we obtain from Eq. (59) 
= (fo 5 (204/42) 
1, OF (oO OF)% (61) 
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The notion of a radiative solution implies that, at a 
fixed point in space, the amplitude of the perturbation 
should die out asymptotically in time on account of a 
finite amount of energy available for radiation in any 
admissible initial perturbation of the plasma. This 
implies that has a negative imaginary part. We con- 
sider oscillations close to the plasma frequency and set 


wo=W,+6—7y, (62) 


where y is positive, and 6 and y are both assumed small. 
On substituting Eq. (62) into Eq. (61) and retaining 
terms to the lowest significant order, we obtain 


0 ae rR 1) 
vy 2 Ale c 


(63) 


We now consider the physically interesting situation 
where the time taken by the wave to traverse the radius 
of the sphere is much larger than the period of a plasma 
oscillation. That is, 


R aT 
peepee: (65) 
where 
R il 
T= z= and Tp — oe (66) 
p 


From Eqs. (63) and (64) we obtain to the lowest signifi- 


cant order 
T,\3 /\)2 
v=o (2) fa-+1 +3) 
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This establishes the possibility of radiation by plasma 
oscillations. However, it is not possible to determine 
the intensity of the radiation emitted within the frame- 
work of the linearized theory presented here and in the 
absence of any knowledge of the excitation mechanism 
for the oscillations. 

In the preceding calculations it was tacitly assumed 
that w? + w,”. There are, however, situations for which 


(67) 


and 


(68) 


wo” —q@,”. It will be recalled from § 2 of Section TI 
that within any region of a cold plasma of uniform 
density one can construct arbitrary wave packets of lon- 
gitudinal waves characterized by 


V x E=0 and V-E+0, (69) 


which are non-dispersive and which do not couple to 
the radiation field characterized by 


We =. (70) 
The longitudinal oscillations, of course, can exist also 
in the interior of the sphere considered here, but clearly 
they do not radiate. 

Those radiating modes involving the functions F 
and G, discussed previously, are associated, however, 
with a fluctuating surface density and do radiate. 

The sharp transition between plasma and vacuum 
assumed here is merely a convenient mathematical 
device for representing a rapid, but continuous, tran- 
sition in density. Clearly the radiation associated with 
the fluctuating surface charge can be viewed as arising 
from a coupling between longitudinal and transverse 
fields induced by the density gradient. In fact, one can 
verify readily that for the case where the density is a 
function of position, the decomposition into solenoidal 
and irrotational waves is no longer possible. For the 
general case, however, the calculations become much 
more difficult and the radiation associated with charge 
separation has, as yet, been explored only incompletely. 


3. OSCILLATIONS IN NON-UNIFORM PLASMAS 


We now consider the longitudinal electron oscillations 
in non-uniform plasmas in which temperature and/or 
density gradients are maintained by some agency. The 
model we take for the steady state is one in which the 
electron velocity distribution is assumend to be iso- 
tropic against a fixed background of positive ions, 
resulting in an electric field which balances the pressure 
forces in the equilibrium state. The force balance con- 
dition 1s, then, 


Vp=—WNeE, (71) 
where p (= N@) and E denote the pressure and the 
electric field in equilibrium respectively. It is also 
assumed that there is no magnetic field. If we assume O 
to be constant, Eq. (71) can be integrated to give 

N= Nesp (72) 
where @ is the electrostatic potential given by E = — VP 
and J, is a constant. Eq. (72) is the well known density 
law for an isothermal atmosphere. On the other hand, 
if we assume NV to be constant, Eq. (71) gives the energy 
integral 


In general, it follows from Eq. (71) that 
p=f(N), (74) 


where f is an arbitrary function of its argument. 
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The equations governing the small departures from 
the state of equilibrium are (cf. Eqs. (I-47), (II-48), 
and (II-119)) 


= He Vin 20; 


(75) 
BOW) | 
m IN Ww ————— V Py — nek = NeE,, (76) 
and 
0 
sp Pit V-(vp)i+pVv+p(Vv)7=0, (77) 


where # denotes the perturbation in the electron density 
and p,, the perturbation in the stress tensor. The per- 
turbation in the electric field E, for the longitudinal 
oscillations under consideration is given by the Poisson 
equation 


V-E, = —42en. (78) 


Introducing §€, the Lagrangian displacement of an 
element of the fluid, Eqs. (75), (77), and (78) can be 
integrated to give 


CN VG) (79) 
; PizetaeeN (Ep) Lap (V Ei pV.E) 2, 180) 
E, = 47 NE. (81) 


Using the foregoing results in Eq. (76) we obtain 
— 


mN°? == pV 5 (pe) +5¥p (§-V>) 


== af ValE-VN)-+V. [pVé+ > (V ENE aa Viasaes 
(82) 


We now multiply Eq. (82) by &* and integrate over 
the volume of the plasma. Assuming that § vanishes at 
the boundary, Eq. (82) leads to 

2 


ot fmargiede=\] arcl¥ (28) le Vell 


+|(G) BVI] +2700 Bom VE + mNeng?| Bae 
(83) 


Now all the terms except the second on the right hand 
side of Eq. (83) are positive definite. Thus the system 
will be unstable only if the second term is the dominant 


one on the right hand side of Eq. (83). 


3.1 Density gradients 


Let us first specialize to the case where there exist 
density gradients in a plasma at uniform temperature. 
Then f = NO, where @ is a constant and the second 
and the third terms on the right hand side of Eq. (83) 
cancel each other. Consequently, w? is always positive 
definite and the system is s/ab/e. 


3.2 Temperature gradients 


We now consider the case when @ = @ (x) and N 
is constant. Then the third term on the right hand side 
of Eq. (83) vanishes. For instability to occur the second 
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term must be capable of having a large negative value 
compared to the fifth term. In order of magnitude, if / 
denotes the characteristic scale of variation of 9, then 
for instability we must have 


NO . 
m Nw? < aa (84) 
Eq. (84) can be written in the form 
ee ear ce Ve) 
i MW =a Co) 


Thus, for instability to occur we must have significant 
temperature gradients over lengths smaller than the 
Debye length. However, at such small lengths the 
macroscopic theory is inapplicable. 

We may thus conclude that electron oscillations in 
non-uniform plasmas with arbitrary density or tem- 
perature gradients are macroscopically stable. 


V. Non-Linear Oscillations 


We now consider large amplitude electron oscillations 
in a uniform plasma and, for the sake of simplicity, 
restrict ourselves to one-dimensional oscillations only. 
Furthermore, we neglect the thermal motions of the 
electrons. Let x, denote the equilibrium position of an 
electron and X(x,) denote its displacement (not neces- 
sarily small) from the equilibrium position. The displaced 
position of the electron is then given by 


x= xy + X (x) (1) 


In moving the distance X (x,) the electrons in the x, 
plane pass over an amount of positive charge eX per 
unit area, where NV denotes the equilibrium density of 
the electrons. If the ordering of the electrons in the 
x-direction is not changed, then all electrons which 
were originally on the positive side of a given electron 
remain on the positive side, while all those which were 
on the negative side remain on the negative side. We 
shall assume that the system is of infinite extent and that 
no charge enters or leaves the system at infinity. The 
equation of motion of an electron is given by 


a2 X 
iy = — CHE, (2) 


where, by Gauss’s theorem, the electric field is given by 
AI ee. (3) 


From Eqs. (2) and (3) we obtain for the equation of 
motion 


aX 


ee gels (4) 


where w,,” = 4.Ne?/m. The general solution of Eq. (4) 
is given by 


X = X4(x9) sinw,t + X4(x9) cos w,f, (5) 


where X, (x) and X, (x9) are some arbitrary functions 
of x). Thus each electron executes simple harmonic 
motion about its equilibrium position, independent of its 
amplitude and the motion of the neighboring electrons 
so long as the ordering of the electrons is maintained. 
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The condition that the ordering is maintained can be 
readily derived in the following manner. Let x) + eG 
denote the displaced position of an electron originally 
at x, while x, +Ax,+ X +A4X, denotes the 
displaced position of an electron originally at x + PAF 
In order that the ordering be maintained, it is clear that 


‘Ax, + AX > 0; we thus obtain 


0x 
5 Ss (6) 
x0 


as the required condition. 
An interesting special solution to Eq. (4) arises when 


X,=0and X,= Asin kx, (7) 


where A is a constant and & denotes the wave number 
of the disturbance in the x-direction. The solution of 
Eq. (4) is then given by 


i= Ko Bid A sin kx COS Wt (8) 
and 
E=4nNeA sin kx cos yf. (9) 
At ¢ = 0 Eqs. (8) and (9) reduce to 
x =x) + Asin kx (10) 
and 
B=4nNeAsinkx- (11) 


We may solve Eq. (10) for x, as a function of x (this 
is most simply done numerically) and thus obtain F as 
a function of x. The ratio £/2.., 1s plotted in Pie. 45 
as a function of x for various values of A. For small 
values of A the curves are essentially sine waves. For 
large values of A one gets a multiple valued curve which 
is, of course, meaningless. What has happened is that 
the ordering of the electrons has been destroyed. The 
crest of the wave falls into the trough and we have a 
wave breaking process somewhat similar to the breaking 
of water-waves, and we may expect the waves to be dis- 
sipated by this breaking mechanism converting the wave 
energy into the thermal energy of the electrons. 


Fig. 16 The behavior of E/E yay as a function of x for various 
values of A at time ¢ = 0, 


A convenient mechanical model which exhibits the 
breaking process is a series of elastic pendulums con- 
strained to move along a straight line. For small dis- 
placements the pendulums swing independently of each 


other. However, if the amplitude of the displacement 
becomes large, the bobs of the pendulums strike one 
another, destroying their periodic motion. 

The breaking of large amplitude waves has been 
investigated by Dawson for the one dimensional case 
under various initial conditions. Firstly, consider the 
case when the initial displacement vanishes and the 
initial velocity is given by 


== Va sink xy. (12) 
From Eqs. (6), (8), and (12) it follows that the critical 
breaking amplitude is 


Ve= oI k. (ars) 
In his numerical calculations Dawson investigated the 
case when the amplitude of the wave was taken to be 
7% larger than the critical amplitude given by Eq. (13). 

Secondly, we consider the case when, in addition to 


the initial velocity given by Eq. (12), there is super- 


imposed on the particles a certain amount of random 
motion. The random motion was taken to be 


Vin = Va exp (—V2/V,2) (14) 


whre I’,, is a constant which is proportional to the 
thermal energy. The effect of thermal motions on the 
breaking of waves is exhibited in Fig. 17 for the case 
where the thermal energy is 4°% of the wave energy 
The lower curve in Fig. 17 shows the wave form for the 
wave at a time two plasma oscillations later. It is clear 
from this curve that the amplitude of the wave has 
reduced considerably and that the wave energy is being 
converted into the thermal energy of the electrons. 
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Fig.17 The behavior of E as a function of x when the amplitude 
of the wave is 7% larger than the critical amplitude at which the 
wave breaks, for different initial conditions. The crosses show 
the wave form for a cold wave just as it is starting to break. The 
squares show the wave form for a wave with initial random energy 
equal to 4% of the initial wave energy, also at the start of breaking. 
The circles show the wave form for a warm wave two plasma 
periods later. 


Fig. 18 shows the velocities of 270 pendulums after 
three oscillations for an initially cold wave. It is apparent 
that most of the particle velocities lie on a sine wave. 
However, there exist a number of particles with very 
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Fig. 18 The velocities of the particles after three oscillations. 


high velocities, both positive and negative. These are 
particles which managed to get in just the right phase 
with the wave so as to ride it, and they account for the 
major loss in the energy of the wave. 

In order to get a quantitative idea of the energy loss 
of the wave, we Fourier-analyze the displacement and the 
velocity in terms of the equilibrium position, i-e., 


N 
J . AIX, 
A (505) Se sia (15) 
n=1 
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where L denotes the linear dimension of the system. The 
Fourier co-efficients A, and B,, are given by 
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Fig. 19 The ratio of the energy in the first harmonic to its initial 
value for a number of different cases. 
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For the Fourier analysis it is convenient to use the 
equilibrium position rather than the actual physical 
position since for such a choice there is no coupling 
between the different modes without breaking. The 
results of the Fourier analysis are shown in Figs. 19 
and 20. Fig. 19 shows the ratio of the energy in the first 
harmonic to its initial value (for a number of different 
cases) as a function of time. The cold wave loses very 
little energy during the first oscillation and a half, after 
which there is an almost continuous decrease in energy. 
As is clear from the diagram, the presence of a small 
amount of random motion causes the wave energy to 
drop to about 10% of its initial value over a period 
of one oscillation. 

The absolute values of the amplitudes of the various 
Fourier modes for the case of forty-five pendulums 
without any random motion is shown in Fig. 20 after 
2.25 oscillations. As is clear from the figure, there is no 
tendency to feed energy from the fundamental into any 
one of the higher modes. The curve is flat (with only 
small fluctuations superimposed on it) all the way out; 
only the fundamental has appreciably more energy than 
the rest of the modes. The energy of the wave seems 
to go into individual particles rather than into Fourier 
modes. 

One may expect the system to exhibit a Poincaré 
cycle and return to its initial state after a sufficiently 
long interval of time. Simple probabilistic arguments 
indicate that this interval of time will indeed be very 
long even for a small number of charge sheets. For 
example, such arguments indicate that a system of five 
sheets should spend approximately 1/400 of the time 
in a state for which the energy of the fundamental is 
more than 95% of the total energy. A system of nine 
sheets would spend only one part in 1.6 x 10° of its 
time in such a state. 
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Fig. 20 The absolute values of the amplitudes fot various Fourier 
modes for the case of forty-five particles per half wave length after 
2.25 oscillations. 


The recurrence time was investigated by Dawson 
following systems of five sheets and nine sheets for 
about sixty oscillations. The systems were started in the 
fundamental mode with amplitudes larger than the 
breaking amplitude. It was found that on the average 
there were five crossings per sheet for the five-sheet 
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case, and ten crossings per sheet for the nine-sheet case, 
in one-half an oscillation time. The state of the system 
was sampled for 120 values of the time. It was found that 
the five-sheet system returned to a state in which the 
fundamental had 95°%, of the energy once during this 
time, while the nine-sheet system returned to a state in 


’ which the fundamental had half the energy once during 


this time. Fig. 21 shows the sum of 120 time-averaged 
velocity distribution functions for the nine-particle 
system. The smooth curve in Fig. 21 shows what would 
be expected if the motion was roughly ergodic. These 
results seem compatible with the probabilistic estimates. 

The foregoing discussion indicates that the breaking 
of waves can be aneffective mechanism for the dissipation 
of wave energy into the thermal energy of electrons and 
that the presence of a small amount of random motion 
can greatly enhance this mechanism. 


100 


f°) 
Se ee oR eo] © th eS 
VELOCITY 


Fig. 21 The velocity spectrum of the particles. 
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RAYLEIGH-TAYLOR INSTABILITIES OF A MAGNETICALLY 


ACCELERATED PLASMA 


T.S. GREEN AND G. B. F. NIBLEtTtT 


A. W. R.E., ALDERMASTON, BERKSHIRE, ENGLAND. 


High speed photography has been used to study the compression of a deuterium plasma by an axial magnetic 
field generated by a wide single-turn coil. In this experiment the inductance of the coil is small compared with 
that of the external condenser bank so that the plasma implodes at a constant acceleration of 5 & 10° cms/sec?. 
Flutes which develop on the outer surface of the plasma are interpreted as magnetohydrodynamic Rayleigh-Taylor 
instabilities produced by the inward acceleration of the interface separating magnetic field and plasma. The observed 
growth rates are in agreement with simple theory. It is shown that the ratio of the diffusivities of magnetic field 
and momentum determines whether the plasma resistivity or its viscosity is the dominant mechanism for damping 


the short wavelength instabilities. At the temperature and densities of this experiment the resistivity is primarily 


responsible for the damping and the wavelength of “maximum instability’ 


observed value. 


Introduction 


The rapid compression of relatively dense plasmas by 
means of axial magnetic fields has recently received con- 
siderable attention [1]—[9] principally because this is a 
possible method of heating a plasma to thermonuclear 
temperatures. This paper describes a high speed photo- 
graphic study of such a compression process in which 
the plasma is pinched to the centre of the discharge tube 
by a magnetic field generated by a wide single-turn coil. 
A stationary plasma confined by a uniform axial magnetic 
field is expected to be neutrally stable against hydro- 
magnetic disturbances [10] but in the absence of crossed 
magnetic fields there is no restoring force to prevent the 
development of interchange instabilities of the flute 
type [11]. In the present experiments the plasma col- 
lapses towards the centre of the tube with approximately 


QUARTZ 


> 


calculated on this basis agrees with the 


constant acceleration and flutes which subsequently 
develop on the outer surface of the plasma are inter- 
preted as Rayleigh-Taylor instabilities produced by the 
inward acceleration of the interface between the magnetic 
field and the plasma. These instabilities are analysed on 
the basis of Taylors’ theory and the damping by resi- 
stivity and viscosity is discussed. It is shown that the 
observed growth time and wavelength are in accord 
with simple theory. 


Experimental Arrangement 


A diagram of the experimental system is shown in 
fig. 1. A copper coil 21 cms wide surrounds a quartz 
tube 8.0 cms internal diameter and is connected by a 
single spark gap switch to a bank of capacitors sup- 
plying 3.5 kilojoules at 24 kV. The total circuit in- 
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Schematic diagram of the experimental arrangement. 
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luctance is 0.14 “H of which the coil contributes about 
).02 4H; thus about one-seventh of the voltage appears 
nitially across the coil and since the inductance of the 
xternal circuit is so much greater than the coil induct- 
nce, variations of the latter during the discharge have 
ittle effect on the voltage and current waveform. A 
nagnetic probe at the centre of the coil measures the 
ield on the axis of the tube. 

Discharges in deuterium at initial pressures in the 
ange 50—200 microns have been photographed axially 
s shown in fig. 1 using two high speed rotating mirror 
ameras. The streak camera records motion across a 
liameter and has a writing speed of up to 20 mms per 
nicrosecond. The other camera is fitted with 30 separate 
enses and gives a series of photographs of the whole 
ube cross-section at a framing rate of 8 xX 108 pictures 
er second. The images are stabilised by using a Miller 
yptical system and are 8 mms in diameter. 


Experimental Results 

Fig. 2 shows a streak photograph of a discharge in 
leuterium at an initial pressure of 80 microns. The first 
alf-cycle in which the gas is ionised is not shown on 
he photograph; the gas current starts to grow at the 
eginning of the second half-cycle of the external cur- 
ent and this time is labelled ¢ = 0. The outer edge of 
he imploding luminous sheath is interpreted as the 
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boundary between magnetic field and plasma and has a 
constant inward acceleration of 5 < 1012 cms/sec? lasting 
for about a microsecond. The inner boundary of the 
sheath is taken to be a shock wave which reaches a 
terminal velocity of 5.5 x 10® cms/sec. After the shock 
has reached the centre of the tube the plasma column is 
further compressed by the still rising magnetic field and 
then oscillates for about three periods before it expands 
to the wall when the external field drops. As the field 
changes sign and increases again a new implosion begins 
and the compression process is repeated. The magnetic 
probe on the tube axis becomes luminous when the 
plasma reaches the centre at about ¢ = 1.5 usec. One of 
the flute instabilities can be seen as a luminous streak 
separating from the lower half of the plasma at about 
this time and reaches the wall at ¢ = 2.5 msec. 

Fig. 3 shows framing camera pictures of a discharge 
which is similar to that of fig. 2 except that the initial 
pressure is greater and there is no magnetic probe pre- 
sent. The upper row of photographs (up to 2.40 Secs) 
show successive stages of one discharge whilst the lower 
row consists of photographs selected from two other 
camera runs. The discharges were repeatable in their 
characteristics so that fig. 3 shows faithfully the progress 
of a typical discharge during the first half-cycle and the 
beginning of the second. 

The cylindrical implosion shows up clearly and the 
plasma is seen to be concentrated in a cylindrical annulus 


1-20 


2:40 


2-65 2-89 313 3-38 3-62 


5:86 


5:20 


413 4.40 467 4-94 


ig.3 Framing camera photographs of a discharge in deuterium at an initial pressure of 150 microns taken using a high speed 
meta which viewed the discharge axially. The time in microseconds from the beginning of the first implosion is shown under 
ch photograph. The exposure time of each frame is 0.12 microseconds. 
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with an inner core of trapped axial magnetic field. The 
framing camera has a depth of field of about 10 cm and 
was focussed on the mid plane of the coil. For the first 
two microseconds the plasma is in sharp focus but be- 
comes increasingly blurred towards the end of the half- 
cycle; this is because the plasma is being ejected from 


the ends of the coil. During the new implosion at the 


beginning of the next half-cycle the plasma is again in 
focus. 

Irregularities in the plasma cylinder develop during 
the implosion phase and these subsequently become 
well-defined on the outer surface of the plasma and are 
believed to be instabilities of Rayleigh-Taylor type. The 
irregular inner surface of the quartz discharge tube prov- 
ides nonuniform starting conditions and these perturb- 
ations grow in response to the inertial forces produced by 
the inward acceleration of the interface between plasma and 
magnetic field. The growth time of the instability is a 
fraction of a mictosecond and the typical wavelength, 
or wavelength of ‘“‘maximum instability,” 1s about 
1.5 cms. The shape of the instability—sharp flutes of 
dense fluid separated by round-ended troughs of the 
lighter fluid —is characteristic of the later stages of 
Rayleigh-Taylor instabilities. 

The temperature of the plasma in fig. 3 has been 
estimated in a number of different ways. The shock 
velocity gives a temperature of 3 * 107 °K but this is a 
lower limit since at these temperatures joule heating is 
comparable with shock heating. Moreover estimates of 
the temperature from shock velocity are imprecise in 
this range where the temperature is a sensitive function 
of the shock velocity. The temperature has also been 
estimated from the measured damping of the radial 
oscillations [12] assuming that this damping is purely 
resistive; this method gives an estimated temperature 
of 6 X 104°K. Pressure-balance calculations yield a tem- 
perature of 7 X 104°K and a particle density of 1017 cm ™*. 
Thus none of these estimates is precise but they indicate 
that the temperature is of order 6 * 104° K and the particle 
density or order 101 cai? 


Growth Rate 


The instability of an accelerated fluid interface was 
first studied by Rayleigh [13] and the theory subsequently 
developed by Taylor [14]. Experimentally, observations 
on liquid-liquid interfaces have been made by Lewis[15] 
and Allred and Blount [16] among others. If @, and @, 
are the densities of the lighter and heavier fluids respec- 
tively and a is an acceleration applied normally to the 
interface, then an initial disturbance of amplitude A, 
and wave number & (A =2a/A where A is the wave- 
length) is amplified in the form 


A= Avent 


= a a1) i 

i= |\—|==——— lar 

lr + Or @ 
The interface is unstable only when the acceleration is 
directed from the lighter to the heavier fluid and equa- 


tion 2 shows that the growth tate increases indefinitely 
with increasing wave number (decreasing wavelength). 


(1) 


where 
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Kruskal and Schwarzschild [17] were the first to dis- 
cuss Rayleigh-Taylor instabilities of an accelerated inter- 
face between an electrically conducting fluid and a 
magnetic field and they described in detail the way in 
which the inertial forces produce charge separation and 
particle drift which drive the instability. They showed 
that though the magnetic field produces the plasma 
motion the gross nature of the instability and its growth 
rate are the same as for ordinary fluids. 

We now apply Taylor’s theory to the observed in- 
stabilities though this comparison can only be ap- 
proximate since Taylor’s theory is a linear theory for an 
incompressible fluid of infinite extent and none of these 
idealizations apply in our case. Since we can regard the 
magnetic field as a fluid of zero density equation 2 be- 
comes 


(3) 


In our experiment Ais 1.5 cm anda is5 < 101° cms/sec?. 
Thus the characteristic growth time is 


1 A \4 
-=(5 |? = 2 x 10-*secs 
22a} 


— (4) 
which is of the order of the observed value. 

Conditions in this experiment are particularly favour- 
able for the growth of Rayleigh-Taylor instabilities be- 
cause of the constant acceleration for a period of one 
microsecond, and since the instability grows like e”’ and 
n is proportional to a's, the shorter the time taken to 
reach a final velocity i.e. the greater the impulsive nature 
of the acceleration, the less the growth of the disturbance. 
In the limit of a pure step-function in velocity there is 
no growth at all. 


Damping Mechanism 


Taylor’s theory ignores effects such as surface tension | 
and viscosity which will eliminate or damp out the: 
However, subsequent investi-. 
gators [18]—[22] have considered these effects and shown | 
that surface tension and viscosity act as damping proc- 
esses which produce a wave number &,, of “maximum 
instability.” The effect of viscosity, for example, in-- 
creases with decreasing wavelength because the shearing ; 
forces increase, and the wavelength of maximum in-- 
stability is essentially that wavelength at which the rate ¥ 
of doing work by inertial forces is equal to the rate of: 


shorter wavelengths. 


energy dissipation by viscous forces. 


Ina plasma the instabilities will be damped by the resis- - 


tive diffusion of magnetic field into the plasma as welt 
as by viscosity, and the inertial, viscous and resistive 
contributions to the growth rate # can be written: 


Inertial term: (a R)’l2 
Viscous term: vk 
Resistive term: R14 06 


where y is the kinematic viscosity and o is the electrical. 
conductivity of the plasma. The inertial term represents: 
the rate at which the instabilities grow in response to: 
inertial forces whilst the other terms are the reciprocals: 


of the characteristic times fot momentum and magnetic: 


field to diffuse across a wavelength and are obtained by; 
observing the similarity of the diffusion equations for: 


momentum and magnetic field, a similarity pointed out 
by Batchelor [23]. The relative importance of viscosity 
and resistivity as damping processes is given by the 
dimensionless parameter P where 


[ees 4nov = R,,/R, (5) 


P is the ratio of the diffusivities of momentum and 
magnetic field and is also the ratio of hydromagnetic to 
hydrodynamic Reynolds number. 

The conductivity and viscosity of a fully ionised plasma 
are given by Spitzer [24] and Marshall [25] and we may 
write for deuterium, 


1-16-2% (AT,)*! 


tela m,'l 22 ln A °) 
lis ee 
am, 8whamshelnA @ 


where 7’, and 7; are the electron and ion temperatures, 
m, and m, the electron and ion masses, # is the ion density, 
2 is the electronic charge, ¢ the velocity of light, & is 
Boltzraann’s constant, /“ the dynamic viscosity and In A 
is the logarithm of a collision integral tabulated by 
Spitzer which is approximately seven for our conditions. 
These equations show that the electrical conductivity 
depends on the electron temperature and the viscosity 
on the ion temperature though at our temperatures and 
densities relaxation times are so short (much less than a 
microsecond) that 7, and 7, are effectively the same. 
Moreover collision mean free paths are short compared 
with the Larmor radii so that it is permissible to use a 
simple scalar conductivity and viscosity. 
Inserting numerical values for o and y gives 


SLU 
— = 
n 


P=42z06 i (8) 
showing that P is very sensitive to temperature. The 
Jamping effects of viscosity and resistivity will be ap- 
Sroximately equal when P = 1 and for deuterium at a 
particle density of 10!7 cm? this corresponds to a tem- 
erature of 2 X 105°K. At temperatures above 2 X 105°K 
he instabilities are damped primarily by viscosity and 
it lower temperatures by resistivity. Thus at our esti- 
nated temperature of 6 < 104 °K it is resistivity that 
Jetermines the wavelength of greatest instability. 

Equating the inertial and resistive terms in the growth 
ate gives 


Am = (7/2)8 a * 08 (9) 


howing that /,,, is insensitive to a. Since a is proportional 
o.7,/2, 4,, is inversely proportional to the electron 
emperature. Putting 1,, = 1.5 cms, a=5 x 10! cms/sec? 
nd In A — 7 gives T, = 5.8 X 104 °K, in agreement 
vith the estimated value. 

At higher temperatures where viscosity provides the 
lamping we should expect 


An = Qna* v3 (10) 


o that the wavelength is again insensitive to the ac- 
eleration but since » varies as 7,’/, A,, should vary as 
"ls and should rise rapidly with increasing temperature. 


RAYLEIGH-TAYLOR INSTABILITIES 


Conclusions 


Analysis of high speed photographs of the instabilities 
which develop during the compression of a plasma by a 
rapidly rising uniform axial magnetic field Suggests that 
these are produced by inertial forces resulting from the 
inward acceleration of the boundary between magnetic 
field and plasma. The estimated temperature and density 
of the deuterium plasma is6 X 104°K and 10!7cm-3 and 
under these conditions the instabilities are damped by 
electrical resistivity. The agreement of the observed 
growth rate and the observed wavelength with simple 
theory tends to support this conclusion. In more power- 
ful discharges the instabilities should be damped by 
viscosity and this increases rapidly with temperature. 
The growth of the instabilities would also be limited by 
a more impulsive acceleration. 
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UN CRITERE NECESSAIRE DE SAB Pea HYDROMAGNETIQUE 
POUR UN PLASMA EN SYMETRIE DE REVOLUTION 
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GROUPE DE RECHERCHES DE L’ ASSOCIATION Euratom-C.E.A. FONTENAY-AUX-ROSEs, SEINE, 


FRANCE 


Pour ce travail, on établit une condition nécessaire de stabilité pour un plasma en symétrie de révolution. Elle 


correspond a des déplacements localisés au voisinag 


e des surfaces a pression constante sur lesquelles les lignes de 


force magnétique sont fermées. Ce critére se réduit A celui de Suydam en symétrie cylindrique. 


Introduction 
En géométrie cylindrique B. R. Suydam [1] a donné un 
critére de stabilité qui est en fait un critére vis-A-vis des 
déplacements limités au voisinage d’une surface 0 = cons- 
tante, surface qui dans ce cas est un cylindre. 
En effet nous exprimons le déplacement — sous la 
forme 
E = [E, (r), Eo(r), & (r),] eike time 


On peut dans le principe variationnel donné pat Bern- 
stein et autres [2] intégrer sur 9 et x et minimiser algé- 
briquement sur &, et €, et obtenir 


: . (Qitgme Os - - 
OV em = | 2rd eee Ss ht icme.) t 
f=kB,+ (mr) B, 
ou ) g=kB, — (mr) Bo 
h= (2/r) J, Bo 


& et m étant convenablement choisis, f(r) s’annulera au 
yvoisinage d’un cylindre quelconque dans le plasma; soit 
r =a ce cylindre. Supposons un déplacement localisé 4 
3on voisinage. Nous pouvons alors remplacer r, g(r) et 
b(r) par a, g(a), A(a) et f(r) par (r —a) f’ (a) et Ecrire 


=n posant 7—7r—a 


pe ' , a 
80 = 20 | eee HORA 


—n7 


a perturbation étant limitée entre les 2 cylindres de 
fayOn 7 — 7,7 + 7 OU 7 est trés petit. 

OW peut encore s’écrire en posant ¢ = 7x et en inté- 
yrant par partie le terme en é &’: 


4 +1 

QD 12, ce 

6v= ot ee (x2 £2 4 AE?) dix 
if 
1VEC 2 
i : p 
A = apa — 28h —b (+ Be) = aap 
ft = Bo|(rB,) 


Pour obtenir la derniére forme de A ona tenu compte 
le la relation d’équilibre du plasma: 


p+ B,B’, + BoB's + (Ber) = 0 


Or on peut démontrer (voir annexe) que 


. i 
J (x8 E24 AB) dx >(A44),) | dx 
=I — il 


et qu’il existe des déplacements qui permettent d’appro- 
cher d’aussi prés que l’on veut Pégalité. 
Ainsi pour avoir 6 W > 0, il faut 


IS TT 
pat ele] 20 


c’est l’expression de Suydam. Elle fournit ainsi un critére 
nécessaire et suffisant pour les déplacements localisés 
dans le cas de la géométrie cylindrique. 

En géométrie de révolution la symétrie du probléme 
€étant moins forte, une seule intégration dans l’expres- 
sion de 6W peut étre faite facilement. Pour continuer 
nous ferons alors ’hypothése des déplacements localisés 
aux voisinages de surface a égale pression de type patti- 
culier correspondant au choix de Suydam en symétrie 
cylindrique. 

Nous obtiendrons finalement un critére qui généralise 
le critére de Suydam, évidemment nécessaire car cor- 
respondant a des déplacements possibles, mais probable- 
ment pas suffisant en général pour la stabilité vis-a-vis de 
tous les déplacements suffisamment localisés autour 
d'une surface a pression constante. 


Les équations fondamentales 
Léquilibre dun plasma a lapproximation hydro- 
magnétique est déterminé par les équations suivantes: 


divB=0 
B Xx rot B= — gradp 


En symétrie de révolution, il est commode d’écrire: 


Ben wendy (1) 
ie r 
qui est la solution générale dans ce cas de divB = 0 
r = distance d’un point a l’axe de révolution 
oe = coordonnée le long de Vaxe de révolution _ 
eg = vecteur unitaire perpendiculaire au plan méri- 


dien passant pat ce point 
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Pp — pression du plasma (définie 4 4 a pres) 
T et y = sont des fonctions quelconques. 

La deuxiéme équation d’équilibre fournit les renseigne- 
ments suivants: 7’ et p ne sont fonctions que de , elle- 


méme définie par Péquation suivante: 


apie t 2 5 
a oes 2) 
avec 
0 Lom OF 
pe are ror | ox? 


Dans ce travail, nous étudions la stabilité des équilibres 
régis pat ces équations. 

oO E . , . r , 

Le point de départ sera le principe d’énergie développé 
part I. B. Bernstein et autres [2] que nous prendrons sous 
la forme: 


8, =3| detlO+ (n-8) (i x mp 
+ yp (div 6) —2 (n-B)°(j Xn) (B-V) (3) 
ou OC = rot (& x B) 


n — est le vecteur unité normal 4 la surface p = constant 
€ — est un déplacement atbitraire, satisfaisant a la con- 
dition n- § — 0 4 la limite du plasma. 
La stabilité pour les déplacements qui ne déforment 
pas la surface du plasma sera obtenue si OW, teste tou- 
jours positif. 


Le critére de stabilité 


Nous emploierons le systeéme de coordonnées ortho- 
gonal suivant, introduit dans [2], adapte a Péquilibre 
dont on recherche la stabilité: Les surfaces de coordon- 
neées seront: 

wy = constantes, 9 = constantes, 7 — constantes 

Les surfaces 7 = constantes sont choisies normales 
aux surfaces y et 0 de telle sorte que (y, 0, %) forme un 
systéme orthogonal direct. 

L’élément de volume s’écrit: 


dt = Jdwd0dy 
Dans ces coordonnées, lopérateur gradient s’écrit: 


oh pbs Cone: 
dy | r 00 | JB * dy 


V=7Be, 


B: longueur de la projection de B sur la plan méridien. 
De l’équation d’équilibre et de Péquation (1) on tire: 


J=rotB=jes-+- BI e, 


avec 
ne Oey (J B) 4 
J J oy (4) 
et 
> oie ae 
i pe (5) 
wo Gate ponies et ls 
dp dy 
Un déplacement & est déterminé par ses composantes 
(Ey Eo Ey) 
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En posant 


r Bey =X, ijraV¥,6,/B=ZaY—(T/Z—U 


et Popérateur 


Onn, eo 
TEN ve Yr 
. L’équation (3) s’éctrit: 
Rae caksc 1 . 
Va aA Deods peas 
ei 00 Bee) lh ab Ox | J el 
ale Tie Gi at oy tree 
a ae b2| 2/KX? 6 
rip y Fam ary) ‘a y (6) 
avec 
_ TT d(logr) jad(log BJ) 
ta oes 
re Ow r oy 
TT’ doflog (B/N _,,d(log B/) 
4 po oy Pp ow ia (7) 


Remarquons que K < 0 est une condition suffisante 
de stabilité. 


Cette condition, dans le cas de symétrie cylindrique 
B=BC) f= 27=—y (Aw ecet 7 T= 0 
condition bien connue dans ce cas. 
Posons: 


X => X "cos mb + X,™ sin m6 
U = > (U,"/m) sin m0 — (U,™/m) cos mO (8) 


m 


Z= XZ," cos mO+ Z," sin mO 


m 


On obtient pour chaque # (m~ +0) 


(5 Wyn = Z\awa A (s+ whe X,) 
' [2 ee x,) | t 


ox — re 
P11 dU, 9 J Xe)? 
oy =) 


J We Ox 
fou pod #s)* 
J\mdyi* ~ oy =| 


OL 
mili 2 pleasent | 
JB LU, + 


; 2 

+ os Xs 

OX. j= 2 

+75/U, + S24 x] 
yp | 0 QZ. 2 
$2P 2X) +I + Gt wie Z| 


oy 
pe OL O25 as 
PS UX) +I Ut ¥ — mis 2, 
~2](X P+ X8)(9 (ogy —4 5, (logs) 


L'indice w a été supprimé sur X|X, UU, Z,Z,. pout 
simplifier. 


Posons Ag ae 
Oo = Chg ai iU,™ 


Lm = 2" + 125" 
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ee Aue et 

L’indice m de X UZ sera supprime chaque fois qu’il 
: ies 

my aura pas de risque de confusion. On peut écrite 


2 aa LOZ ihe 2 
(OW 2) m = A dwdy be B2| ox ee pe Xx 
LON OHTANE la OK, 7 oP 
See a | ee | 
eee © As a Oh oe eA 
Play) Pye 4 v ints | 
3 ae eee 2 
= _— or sae SS o ly | 
27 |s,5 (log r) lee (log By) | (9) 


Dans la suite nous utiliserons le résultat suivant relatif 
aux solutions périodiques d’une é€quation différentielle 
du type suivant: 


OW : 
vi + 24(x)Z + b(y)=0 


On suppose que a(x) et b(x) sont des fonctions de 
période yo. Si 
x Ko 
$a(x)dy= \ a(x) dx +2 Ka (Kentier) 
x 
on a une solution périodique, de période ,, et une seule 


1 


1 x Xo 24 
2 (x)= ay § a(x) ae \ b(x) expi | a (u) dudx (10) 


x 


En effet, on peut vérifier directement que cest une 
solution. Elle est périodique. En effet: 


%+2 x0 x 
1 


1 —expz $ a(x) dx \ b (x) exp | a(u)dudx 


L+ Xo L+ho 


7 Xo) = 


et faisons le changement de variable 


4=Yo+” 
Z(x% + xo) 
1 LT ho f 


—— j ——/, 
1—expi$a(x) dx Jou alae a(v)dvd. (x) 

x x 
Enfin il n’y a en pas d’autre car on obtient la solution 
générale en ajoutant a cette fonction périodique parti- 
culiére une solution générale de Péquation homogéne 
(0z/d7) +2a(z) z=), Cest-a-dire 

2 Exp | a (x) dx 
yal 

solutions qui ne sont jamais périodiques. 
. Revenons a (6 WF) ,, et minimisons d’abord sur Z. 
En posant 


Le 
=P 2 xX) +JU+ Soins Z| 


Soit/= f(y) exp im eae a| 


f doit étre périodique en yx (période Xo). Ceci ne peut 
se produire que si 


ieee, Sie 
§ Lin ay=2ka 
$ indique Pintégrale prise sur une période. 


Sauf pour un ensemble de points particuliers (de mesute 

nulle) cette relation n’est pas vérifiée. Donec C0) = 0) 

La minimisation impose donc que f = 0. 

Encore faut-il trouver une solution Z convenable Cest- 

a-dire périodique en x déterminée par f = 0. Ory —0 

est de la forme étudiée ci-dessus 

eee Z+JU+4 o Geo) —10) 
eee oy 

avec 


a(y)=—wm/JT/r 
On Fae 
by) =JU +s I% 
fonction périodique de période y, et 
$ a(x) dx = —m $ (J Tr?) + 2k a en général 


La solution cherchée qui donne le minimum pat rap- 
port a la variable 7 est donc 


il 
~ Texp(—in§ (ITIP) dy 
Ueaigta x 
\ (72 +9 /¥) ex [| \ (7 1e%) du| ds) 


INTERPRETATION DE 9 Cir) ax 


La ligne magnétique tracée sur la surface y = cons- 
tante et qui passe en un point M de cette surface a 

our €quation: 

‘ pe sae or 
LG == oF — 7 

Soit: 

Gea Tire) ay 

00 = $ JT Ir) dy (12) 
Ainsi la quantité $ (/T/r?) dy n’est autre chose que 
Pangle 0 pris entre deux point Met M’, M’ étant le 
point ot la ligne magnétique recoupe pour la premiére 
fois le paralléle passant par (/. Cet angle est indépendant 
du point M. ; 

Les cas particuliers o1 6, =2 ka signifient que la 
ligne magnétique se referme. En général 0, +22 et 
la ligne magnétique ne se referme pas. 

Posons alors 


iq - +U 
e= a eae x (13) 
a oe 7 (4 | x 
Entre les 7 ¢ & il existe la relation: 
- o LiméE =0 (14) 
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Moyennant cette relation, on peut tirer univoquement 
X et UA partir 7 € & périodiques en x en utilisant une 
fois encore le théoréme (10) appliqué a Péquation: 


Ox er 
ox r 
On peut alors écrire: 
ey |e|? p fle 
On Fm == 3 andy | re DB 5 = } Po) xl 


+78 |n+ 4X] —2s\x}x} (15) 
avec: 
=#1 
1—exp[ [—im § 1 (JT |r ) dy] 


aigXo 
je exp [in| (JL fr?) du| ax 


x 


X= 


Remarquons alors que si on pose: 
z, 
0 (x.4) = \ (J TI) da 


e 
a 


on a: 
6 (x + x0, 4) = 9 (x, 4) + % 

donc: 
8(4,4) = p(x54) + 90 x! xo 

avec: 


P(X + 104) =P (x4) 
ce qui permet d’écrire: 
—exp [imp(x,4) + im9o x/ xo) 


— exp(— io) 


xX ——— 
(16) 


| E(x) e—?mp(x,a) eo tMO,t/ x dx 
Comme ¢ (x) e—”? est une fonction périodique en x 
(période yo) on peut développer cette expression comme 
suit: 


€(x) e-imnx,a) — 2 im erinkx/ x, (17) 
Ge Gal reportant, on trouve: 
= es Qink 
toe Ho oe Tae ai aa! ize) 


k 


Remarquons que le changement du point de référence a 
se traduirait par un simple changement des a,,,. Pour 
la suite, nous écrirons p (x) . 

En analogie avec l’hypothése de Suydam [1], nous 
étudierons les déplacements (&, #) tels que 27k — my 
est voisin de 0. 

Soit une surface y = constante telle que 6, /2M ne 
soit pas rationnel. On peut toujours choisir les entiers 
ket m de telle sorte que 27k — m0, soit aussi voisin de 0 
que l’on veut. 

On peut supposer que ces déplacements sont parmi 
les plus dangereux. Ils coincident d’ailleurs avec les 
déplacements de Suydam si lon se place dans le cas 
particulier de la symétrie cylindrique. 
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Aussi allons-nous étudier les déplacements définis par 


gimp) e2inky]| Zo 


E = Aym (19) 
2 Xo 1 
xX —p7, —~— em p(y) etn hy | 9 
Prem i(2k — 9p) 
On constate que: 
= —1Y =f 
ASS eS 
2ak— my y 
AVEC: 
20k — mBq 
y —= ~ — 
Xo 


Posons alors: 
== OG &— VX 17 0e* 
Donec: 


e= ely) 
C= mp7) = 2mke A+ C(y)=m0(x,4) + vx + Cy) 


C (yp) étant une constante arbitraire. 
Introduisons la fonction / telle que: 7 = )X 


LL arity ee 
1 02 L blivX + -—% x| 
=x (2 pbs) 
O€ fs oR é ee ; og] 
ae i id. 0 —~* a y| 
D’aprés la relation (14): 
Dieta O€ 01) _ Ov . a op ah 
CED a * apt ie oy ox 
ivh sl 
NE va as it Ov Vv on 1v OM 
€- 72 A eX 2 ee, Wy, Sar 
70h Y ak 
a Ox A re q 
te 
+ Re (b+ 7p 
D’ou 
sca fi. 1 ov 7v OM 
We \ BNO ae 
7 0A vy aire yo’ 2 ¥ j 
Se m =) b+ = LIB b+ a) —2yK} 
(20) 


Les calculs qui ont permis de passer de l’expression (9) 
de OW ,, 4 expression (20) sont valables en dehors des 
surfaces y = const. sur lesquelles les lignes de force 
sont fermées. Pour un # donné ces surfaces correspon- 
dent a des valeurs discrétes de y. L’équivalence entre (9) 
et (20) est donc établie si ?on exclut de Vintégration 
sur y un petit intervalle autour de chaque valeur singu- 
liere de y. Mais comme les quantités a intégrer dans (9) 
et (20) sont continues en y, les expressions (9) et (20) 
sont aussi égales sans restriction sur lintégration en y. 
L’elimination du terme provenant de y p (div &)? neces- 
site une étude plus détaillée an voisinage des surfaces 
singuliéres de y. 


UN CRITERE DE STABILITE HYDROMAGNETIQUE 


Nous allons chercher un critére local, cest-a-dire que 
nous nous limiterons aux déplacements localisés au 
voisinage des surfaces y = y, et tels que y soit voisin 
de 0. 
ov 
oY 

Ce choix de déplacement est analogue au choix fait 
par Suydam. En effet, le critére de Suydam est en fait 
un critére local au voisinage des surfaces r —r, od y 
cylindrique = & + mB,/rB, (voir introduction). 

Ce y cylindtique correspond exactement au y» défini 
ci-dessus en symétrie de révolution. Les cylindres r — r, 
définis par Suydam se généralisent selon les surfaces 
y = Yo 0u les lignes de champ magnétique sont fermées. 

Ainsi dans 6W; nous remplacerons pat Cul We) ce 
et les quantités telles que r, B etc. ... relatives 4 la solution 
d’équilibre seront prises constantes et égales a la valeur 
qu’elles prennent pour y = yy. Elles dépendront encore 
en général de y. 


"(y)=(9—vo) (50) =(v—vo)o’s (21) 


idh fo) +78 [6+ 4a] —2e4 (22) 
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Nous limiterons encore les déplacements en supposant 
la fonction / réelle. Dans ce cas le seul terme contenant 
\2 
04/07 dans lintégrant est de la forme ae ("| et il est 
m*\ox!) 
alors évident que lexpression sera minimum pour 
m = co; OW ,,. ne contiendra que la fonction / ce qui 
permet de minimiser par rapport 4 elle sans difficulté. 
On trouve: 


Np Be paa lea) 2} 
Wr 3 (eda laa rs eee 
: pe 
Pee Nea pl yl ey 
se P= Yo) 2 0 F r B2 | r 3 
YQ" ere 
+ (Y — Yo) - 3 + er] 2jKe (23) 


En intégrant par partie le terme o 0’ par rapport a y, 
on peut écrire dW... sous la forme Wp... = C?R 


R= J dy [(p— yo)? 0? + Ao? (24) 
avec 
[v'o\2 er? Bt 
C= : dy 
Al 2 , Bps 
ay aa 
A= sel(r'+ Za) 4 + 35 ayK\ , i 
JOURS (TIP) oe 


Or, on a démontré (voir annexe) que pour og (1) = 
o (—1) = 0 et 9 9’ continues 


+1 ; +1 
\ (e* = Ae*)ay> (4+; | eax 
==] asl 


et qu’il existait des déplacements qui rendent la valeur 
du premiére membre aussi voisine que l’on veut de alle 
du second. 


On peut appliquer ce résultat expression OW p.., et on 
obtient le critére suivant: 


» [ae ral( 4 LT \Z 4 top 
aed | ies ial i |e 


En tenant compte de la formule (4), Q.., s’écrit encore: 


r2 B2 J? B Py is y\9 
a ie gyi ee a eae 
Ou=$4x (042) (a ayl 7a2) + 3m) 2K 
(25) 
Dans cette expression 
1 
Oe iwi 07 Hose ae Fr 
mi dy ee os ates 
Te 
en posant pe 
Ee = 
! ow re) se odyx 


et en tenant compte des équations d’équilibre du 
systéme, on peut réécrire O., sous une autre forme: 


ows me Jax 
a Patt Daa RB 
ee Jax 
+70 Va at sO 


| Pe S44, 
PED apa + (PAB) 


Dans ces expressions, il reste encore une fonction 
arbitraire, c’est le paramétrage des courbes y = cons- 
tante, c’est-a-dire que lon peut choisir arbitrairement 
pour x une fonction quelconque de yp. 


4 =F (xo) 
Jd x%=Jy4 Ko (27) 


Il nous reste 4 minimiser sur cette fonction f soumis 
a la restriction: 
bate io) ale 1 (28) 
Dans l’expression (26) seuls les deux derniers termes 
varient dans la transformation y =f (7). La minimi- 
sation est algébrique et la valeur minimum pour ces deux 
termes s’écrit: 


(26) 


4, étant tel que: 


1 
4§ [1+ (22/7? B?)] (J/r?) dz - 
r , 
ls ue dy—2p'T eel 


Jax 
"PP sat PR] 

L’expression finale minimum des 75s ectit: 
oO iy Oneal ts 
oo T2 Jax oy Pe) 

+G(1 eagle 
ee pA Is 
+0 § so dx—p9Gge (29) 
Notons que J dy = ds/B, s designant l’abscisse curviligne 
des courbes a pression constante. 


2 


; Jax 
dy—2p'T WD “aps 


Sli 
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Le critére obtenu s’exprimera donc par: O. > 0. 
Si nous nous placons en symétrie cylindrique on a / 
-1/B,, T=rBo, dx devient dz, 0/dp devient 0/07 
et 6 — JT /r? = B,|(rB,). Toutes les fonctions ne dépen- 
dent que der. On obtient dans ce cas le critére de Suydam 
en utilisant ’équation d’équilibre 
4/2 2p Le 
qt 
it dit 
ea 
° dr 


05 (30) 


‘ef ww signifiant alot 
f dr 


NOTE 

Ce travail terminé, nous avons regu un article de 
B. B. Kadomtsev [3] dans lequel Dauteur généralise 
aussi le critére de Suydam pour une géomeétrie torique 
quelconque mais en imposant a Tavance une forme 
particuliére de déplacement localisé. Le résultat obtenu 
revient dans notre méthode a supposer 4 = constante. 
Dans ce cas nous retrouvons en géométrie de révolution 
son résultat. L’expression (29) est différente et plus 
sévére cat nous avons supposé seulement que / est 
réelle et minimisé dans ce cas. 

Le critére de B.B. Kadomtsev 
notations: 


s’éctit avec nos 


1 i 
w &/ JE% (4 4 [b+ aS: JT 2|r2 dy, 
Oe 6 /B dy | 
: ey Rage | 
<6 ‘| e dy— 2pT? e J i 
is § JB dy 


($ Jaz) 


/ 20 / 2) 
TED apt —P" gia 


Annexe 
Il s’agit de démontrer que la quantité 


est toujours supérieure a 1/4 pour toutes fonctions y 
s’'annulant 4 -— 1 continues ainsi que leurs dérivées, la 
valeur 1/4 pouvant étre approchée d’aussi prés que l’on 
veut par une fonction de cette classe (C). 

Nous indiquerons les grandes lignes d’une démonstra- 
tion basée sur le développement de y en polynomes de 
Legendre, orthogonaux sur (— 1, + 1). Prenons pour 
une combinaison linéaire des NV premiers polynomes 


om (x) 


N 
b= wanes) (A1) 
il 
N 
y(—1) = > (—1)" a,=—= 0 
1 
Ae 
I-41) => 2, = 9 
il 
D’ou 
INV N 
rye Sa, 0 
pe Ba, n= 0,2, 4 
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En utilisant les propriétés intégrales des polynomes 
de Legendre on peut écrire: 
IN (n—m = 0, 2,4...) 


Css - =e aN om >) Ay? sees ay (72 = 1) 
32a) + 1) 
tes, 
aes te) 


R R R {Ry, ayant que des » pairs 
y= Byy T Bx (Ry; ayant que des » impairs. 
La suite traite R,,,, le calcul pour R,, est identique. 


En tenant compte de Dy 4, = 0 On peut poser: 


gn 
dg = X4— Xe 
an Xn +2 Xn 
Eh = OS 35 
Ry s’éctit alors: 
N N—2 
a (20 1) x2, + 2 (24+ fy) Seeoee aks 
n= n=4 


L 


forme quadratique que l’on peut décomposer ainsi: 


2n+1 . 
> pile = : 2 
Ry = Xn E n 20 Xn +2 + Oy Xy 
Ae 
(= 7 
AVEC, —= DEL 
(24 — 3) 
Ci = 2n =o 
ZEN 
ce qui fournit 
2n+1 il 
on = Bi == aa. 


Xx, est donc toujours positif ce qui prouve que Ry Vest 
toujours aussi. 

Ainsi: S> 4 pour tout ydelaforme (A1). S’ilexistait une 
fonction de i classe (C’) telle que S (y,) < 4 on pourrait 
trouver un_y de la forme (A 1) tel que [S (y) — S (,)] soit 
arbitrairement petit, donc tel que S (y) < 4 ce qui est 
impossible. 

Ilsen) résulte: que af 9) 23 
classe (C). 

Pour montrer que l’on peut approcher d’aussi prés 
que lon veut la borne 4, choisissons pour les x, la 
forme suivante: 


iad 1 
xy = (—1)”/2 ce —77] tn 


pour tous les y de la 


avec P > 0 et 
more 2n+l 
om 2a, - avec f, a ()) 
Alors 
N-2 1 ° N—2 
il x 
= rh |e es —. 
se : oe (a? ad nies 2PHi 


COL Onsen 


D’auttre part: 


Z, 
D,= igece Le ts. 
> 2n+-1 


N—2 
=a 
—= Gee: = 2 
OES A rae 
i 
OCA ie ire Sc 1 \72 
tae oe feeaee | 
N—2 
BiNae eft La 
Dy > 3 -s le al cat 2u+1< 3n 
N—2 
D 2 
»>3|> an aH a 
N—2 5 
Zz ze De A 
Or ye > ah Dae (P41) 
2 1 


€ — fonction de Rieman. 
Or on sait que 


E(P +1) = -- fonction entiére 
D’ou 


N—2 


UN CRITERE DE STABILITE HYDROMAGN ETIQ UE 


Ainsi pour P >0,& (NV) tend vers zéro lorsque NV + co 


] 
N—2 
1— K (N) > Nae aa 


n=2 


Ry <12/P2 
De 


Done Ry/Dy pour N assez grand sera inférieur a 
disons 24 P? et ’on peut donc dctermine: des fonctions y 
telles que Ry/Dy soit aussi voisin de 0 que lon veut 
en choisissant P > 0 assez petit. 
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RAPID COMPRESSION HEATING OF A PLASMA IN THE LINEAR Z-PINCH 


H. A.B. Bopin, A. A. NEwTron AND N. J. PEACOCK 


A.W.R.E., ALDERMASTON, BERKSHIRE, ENGLAND 


The electrical characteristics of discharges with high and low initial rates of current rise are briefly described and 
the electrical measurements ate used to determine the work done on the plasma by rapid compression. The particle 
and the magnetic energy in the discharge are discussed for both high and low initial rates of current rise (dl di ©) and 
for discharges with and without an axial magnetic field (H7,). It is shown that the gas is heated much more efficiently 
both with and without H, at high d//d¢ (0). For dZ/dt (0) ~ 1.5 - 10" amps/sec (H7, = 0) 15% of the initial bank 
energy goes into heating the gas and 10% into the magnetic field at the first maximum contraction ; for al/dt ~ 10" 
amps/sec (H, = 0) 11% goes into the field and only 6.5% into the gas at this time. The energy input to the gas is 
found to become limited well before peak current, and the initial implosion is not followed by much of the 
subsequent isentropic compression heating which might be expected, so that the temperature never exceeds about 
300 eV; it follows that most of the observed neutron yield (~ 108) is not thermonuclear. Data obtained from the 
electrical waveforms and from high-speed photographs are presented and analysed, and it is shown that large currents 
build up outside the central plasma column after the first maximum contraction. Further evidence suggests that 
these currents flow in wall material boiled off by radiation from the central column. These wall effects account for the 
absence of further isentropic compression heating and also contaminate the central plasma column; they are such 
as to prevent any significant gain in heating by increasing the power input to the discharge. It is concluded that 
similar effects probably account for the failure to achieve a thermonuclear temperature in many high power pinch 


experiments. 


1. Introduction 


The rapid heating of a plasma is desirable for a number 
of reasons. For example, if a cold gas is heated up it 
passes through a stage when it is only partially ionized 
during which time it loses energy very rapidly by transi- 
tions involving bound states*. It is desirable to minimise 
the duration of this phase by making the power input large 
compared with the radiated power. Also, the con- 
tainment problem is simpler at high temperatures. 
Heating by rapid magnetic compression (see, for 
example, Andrianoy et al. [1], Bickerton [2], Fitch [3]) 
is one of the most promising methods of obtaining a hot 
plasma. The simplest geometry for this approach is 
the linear z-pinch. 

In earlier experiments on fast linear x-pinches carried 
out at A.W.R.E. the neutron emission was studied at 
high intial rates of current rise (d//d¢(0)); this work 
showed that even from high-power discharges the 
neutron emission was not thermonuclear (Bodin et 
al. [4]). There followed a more detailed investigation 
of such discharges, with both high and low values of 
dl/dt(O), (hereafter described as “fast” and “slow”), 
both with and without an axial magnetic field (H,). 
The main qualitative features of fast and slow discharges 
and a few of the measurements on rapid compression 
heating, were reported by Bodin et al. [5]. 

In this paper, which is divided into two main parts, 
the experiments on rapid compression heating are 
described in detail. The first part discusses the energy 
put into the gas and the magnetic field up to the time 
of the first maximum contraction. It is then shown that 
in high power discharges the energy input to the gas by 


* This has also been pointed out recently by J. L. Tuck in Phys. 
Rey. Letters, 3, 313 (1959). 
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compression never significantly exceeds its value at that 
time. The second part of the paper presents and analyses 
results which show how the heating is limited by wall 
effects. The consequences of these effects are discussed. 


2. Apparatus and method of measurement 


The experiments were carried out using the Maggi 
condenser bank, described by Fitch and McCormick [6]. 
The capacity used was 100 wf and most of the measure- 
ments were made at 20 kV. The total inductance, L,, 
external to the gas discharge was 5 x 10-° henries for 
fast discharges and 200 x 10°® henries for slow dis- 
charges, giving d//dt(O) in the gas of 1.5 x 1012 and 
1011 amps/sec respectively. A clear quartz discharge 
tube, 34 cms long, 15 cms bore and 8 mm wall thickness 
was used, with plane copper electrodes. Weak preioniza- 
tion was always used in the experiments described. 
Streak photographs were taken and the voltage (I/), 
current (/) and d//dt were recorded. 


3. The electrical characteristics of discharges with 
high and low initial rates of current rise 


Examples of the voltage and current waveforms are 
shown in figure 1 for fast (A) and slow (B) discharges. 
A (above) shows a typical current oscillogram for a 
discharge in which the bank inductance is comparable 
to the initial tube inductance, but much less than the 
final tube inductance. In this case the current rises very 
tapidly to a high value (360 kA in 0.6 / sec), and then falls 
to a lower value (208 kA) at the first maximum contrac- 
tion (first pinch) of the plasma. This form of the variation 
of current with time is desirable because (see Fitch [3]) 
at the first maximum contraction more energy can go 
into the gas than into the magnetic field associated with 
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Fig. 1 
A. dlI/dt (0) ~1.5- 10! amp/sec, Vz = 20 EV. fg = 320 Ds 


the discharge. (See section 4.) The voltage is not 
sensitive to changes in the tube inductance when the 
bank inductance is small, and the oscillogram in this 
case (A) dces not show any of the rapid fluctuations at 
the time of the pinch which are characteristic of a 
high-inductance bank (see B). The high frequency 
oscillations seen on the voltage oscillogram, and to a 
lesser extent on the current, are due to cable reflections 
in the bank (Fitch and McCormick [6]). The current 
in the slow discharge (B) is predominantly determined 
by the large series inductance and shows only slight 
evidence for pinching. In this case the current at the 
first maximum contraction is 190 kA. 

The application of H, results in several apparently 
stable oscillations of the plasma, as found by Bezba- 
tchenko et al. [7] The oscillations were seen on streak 
photographs and also on the current waveform for a low 
bank inductance (Fig. 2). The peak current is larger 
for H, =0, contrary to predictions based on simple 
inductance considerations. (See also section 6.) 


596 kA 


S7SkA 


Hsec 


Fig. 2 Tracings of typical current waveforms with and without 
H,. Po = 350 uw Dg, Vz = 20 kV, dl /dt (0) ~ 1.5 - 10)? amps/sec. 
iN, fT, = 0 B, E77 = 1480 oe 

4. The particle and magnetic energy in discharges 


with high and low initial rates of current rise 


THE METHOD OF CALCULATING THE DISCHARGE ENERGY 
The total amount of electrical energy put into the 
discharge tube up to time, 7, is 
t 


) Vidt (1) 
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Tracings of typical current and voltage waveforms at high and low d//dt (0), H, = 0. 


B. dl/dt (0) ~ 101" amps/sec, V, = 21 kV, py = 350 wu Dy. 


where |’ is the voltage across the tube and / the current. 
Using equation (1) it can be shown* that at time ¢ 


\VIdt=4L i244) 2 (dL]di di (2) 
0 


0 


where L = tube inductance. 

By definition } LJ? is the energy in the azimuthal 
magnetic field, and therefore the last term must represent 
the total mechanical work done on the system as a result 
of the motion (due to the magnetic forces) of current- 
carrying conducting material, i.e. as a result of changes of 
inductance. The total mechanical work done on the 
discharge is 


W—=+\ 12 (dL|dt)dt (3) 


Equation (3) is independent of the geometry since W 
includes the work done due to the motion of all the 
currents in the tube. Where all the currents flow in 
a central plasma column all of the mechanical work is 


done on this column. 
For discharges with no axial magnetic field W goes 
into particle energy W, and 


Vee (4) 
For discharges with H/, 
V,=W-—W’ (5) 


where W’ is the energy required to compress the 
trapped field. 

The total magnetic energy, W,,, in the discharge is, 
for unstabilized discharges, 


W,=Wyy=iLlP (6) 


The total magnetic energy in the discharge for dis- 
charges with H, is 


Wy ~ W no a Aen (7) 


where IW,,, is the energy in the trapped axial field. 


* Resistive heating, which can be shown to be small in times of 
interest, is neglected. 
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AND MAGNETIC ENERGY IN THE DISCHARGE 
FIRST MAXIMUM CONTRACTION 


‘THE PARTICLE 
UP TO THE TIME OF THE 
Discharges without axial magnetic field 


Figure 3 shows the particle and magnetic energies 
in the discharge as functions of radius up to the first 
maximum contraction for the fast and the slow discharges. 
In the fast discl harges it is found that whereas IV, rises 
steadily to a maximum at the first maximum contraction, 
W,,, levels off so that W/W, = 1.5 at this time, i.e. more 
energy is put into the gas than into the magnetic field. At 
the first maximum contraction about 15°% of the total 
energy initially in the bank (20 kJ) has gone into the gas, 
and about 10° into the magnetic field. 

In the slow discharge IV,, is always 
and both increase steadily as the compression proceeds. 
At the first maximum contraction 6.5% of the energy 
in the bank (20 kJ) has gone into the gas and 12% into 
the magnetic field. In the slow discharge most of the 
work is done at small radii, in contrast to the fast 
discharge in which a substantial amount of energy is 
put into the gas during the early part of the compres- 
sion, when the current (Fig. 1) is high. 
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Fig. 3 Total particle energy (IV) and magnetic energy (W,,) in 
the discharge as functions of radius at high and low dl lat (0) 
without axial magnetic field. 


Thete is a gain by a factor 7 in mean power input and 
by a factor 2.3 in W, for the fast discharge. An increase 
in W, by a factor of only 2.3 for a factor 15 increase in 
the inital d/ /d¢ is not unexpected when d//dt (0) is changed 
by varying the bank inductance, using a relatively low 
voltage. It arises because the very high initial d//d¢ in the 
low inductance case is only sustained for a small fraction 
of the compression time near the start while the tube 
inductance remains comparable to the bank inductance. 

At the first maximum contraction the plasma radius 
in the fast discharge is larger by a factor 2.1 than that 
for the slow discharge; this is a consequence of the fact 
that more energy is put into the gas in the former case 
during the earlier part of the compression, and the 
plasma pressure limits the contraction at a larger radius. 
Another example of this is seen inthe streak photographs in 
fig. 9, for Hl, = 0, which show that the plasma radius at the 
first maximum contraction increases with the voltage. 
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the discharge as functions of radius at low and high d//dt (0) with 
axial magnetic field, H7,. 


Discharge with axial magnetic field 


Moderate axial magnetic fields were applied to reduce 
the fastest growing modes and improved stability was 
found (Bodin et all. [5]) for fields ~ 1000 oe. Fig. 4 
shows the particle and magnetic energy (including that 
in the trapped field) in the discharge as functions of 
radius up to the first maximum contraction for the fast 
and the slow discharge with H, = 
that with an axial magnetic field 7, is always less than 
W,, 11 both cases, although only by 13% for the fast 
discharge. Comparing figures 3 and 4 it is seen that the 
application of 7, does not have a marked effect on W, 
and most of the magnetic energy is still in the H, Geld) 


Ion temperatures 


The values of W,, and ion temperature, 7, at the first 
maximum contraction are shown as functions of initial 
pressure for the fast discharge in fig. 5. For H, = 0 the 
concept of temperature is fully valid in most cases and 
is always approximately valid, as ion thermalization will 
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Fig. 5 Total particle energy (IV,) in the discharge and ion tem- 


perature (7’;) as functions of initial pressure at high dl/dt (0) with 
and without axial magnetic field, Tol. 


cur in ~ 0.2 usec in the longest case given (7, ~ 300eV 
ut intitial pressure ~ 25 uw). For H, — 1480 oe the 
tapped field will delay the thermalization unless, for 
-xample, there is a shock wave, and values of 7, are not 
siven for this case. 7, is calculated from W, assuming 
10 losses and that all fe energy is in the i fon i.e. 


T~,.=- 


2 


W, |.N/ 


Oo! bo 


where V/ = number of deuterons in the tube. 

It is seen that the highest temperature which could 
have been achieved at the first maximum contraction is 
300 eV at 25 uu, the lowest pressure at which a discharge 
can be initiated with a bank voltage, Vz, of 20 kV. 


THE WORK DONE ON THE GAS UP TO PEAK 
CURRENT-MAXIMUM TEMPERATURES ACHIEVED 


Figs. 6 and 7 show W, and W,, as functions of time 
for peas: == (0) eine! 80) we respectively *. In the fast 
discharge W, tises to a maximum at the time of the 
first maximum contraction, and then falls as the plasma 
expands against the field; as the plasma contracts for a 
second time W,, begins to rise again but it is found that, 
for AZ, — 0, it never exceeds the aere reached at the 
first maximum contration. For H, = 
it only exceeds the early value by 20°, during the ee 
g2ces. Thus, for 7, — 0 there is no further - compres- 
sion heating after the first maximum contraction, while 
for HH, = 1480 oe there is only a little subsequent 
heating. The observation that the plasma temperature 
never significantly exceeds that achieved at the first 
maximum contraction is general in the present experi- 
ments and is discussed in the next section. 
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Fig.6 Total particle energy (IW,) and total magnetic energy (WV) 
n the discharge as functions of time at low and high dl/dt (0), 


H, = 0. 


‘ W, after the first maximum contraction includes a small 
Daecibution due to currents which are found (see section 5) to 
stow outside the central plasma column after this time. At earlier 
imes all W, goes into the central plasma column. 


RAPID COMPRESSION HEATING OF PLASMA 


The thermonuclear yield expected from the calculated 
ion temperature in these experiments (see preceding 
section) is only 108?—10*, which confirms the cee 
reached in the study of the neutron emission (Bodin et 
al.[4]) that most of the observed yield (~ 105) was not 
thermonuclear. 
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Fig. 7 Total particle energy (W,,) and total magnetic energy (IV) 
in the discharge as functions of time at low and high d//d¢ (0) 
with axial magnetic field. 


5. The heating limitation in high-power Z-pinches 


The processes limiting the heating were revealed 
when data from the streak photographs and the elec- 
trical waveforms were compared. 


THE METHOD OF ANALYSIS OF THE STREAK 
PHOTOGRAPHS AND THE ELECTRICAL MEASUREMENTS 
The inductance L () of the discharge can be obtained 

from Il’ and J by the equation 


V =(digs) (£7) (8) 
which can be written 
yea 9) 
L(?) 


Resistive terms have again been neglected because, 
except in the early sheath formation stages and in the 
later stages when impurities from the walls enter the 
discharge, the current and voltage waveforms are con- 
sistent witha purely inductive impedance of the discharge. 
Furthermore, if resistance is included the main conclusion 
(that /, can substantially exceed 7, — see below) would 
be strengthened. 

Hine plasma radius r() can be found from the induct- 
ance using the relation 

L (t) = 2/log [R/r (4)] * 10-* henries (10) 
where R = radius of outer coaxial current return and 
7 —tube length, on the assumption that the current 
flows ina thin sheath. (The inductance for a thick sheath or 
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for a uniform current is larger.) The radius of the current 
sheath found in this way from the electrical measurements 
is denoted by r,. The radius of the outer boundary of the 
central luminous column (r,) is determined directly 
from the stteak photographs. Values of the current 
corresponding to r, and r, ate, respectively, /, (measured 
directly by the monitor) and /, (determined from the 
streak measurements, using equation (10) to find a value 
of L(A) from, and inserting this value and the observed 
voltage into equation (9)). 7, only equals /, if all the 
current flows on the outside of the luminous central 
column; in this case r, and r, are also equal. 


Results obtained from analysis of the streak and 
electrical measurements 


Fig. 8 (A, above) shows J/,, 7, and also the current 
given by theory as functions of time for three different 
discharges; fig. 8 (B, below) shows the corresponding 
values of r,, 7, and theory. The initial gas pressure is 
arranged so that in all three cases the mass of gas, M, 
per cm length of the tube is approximately the same, the 
voltage is 20 kV in each case so the pinch-time scaling 
factor (1”2/M)t is the same for the three discharges to 
a good approximation. The theoretical curves were 
computed for the present experiments by Taylor and 
Nisbet [8] from the Rosenbluth [9] snowplough equa- 
tions, taking into account the effect of condenser-bank 
inductance and H7,. 

It is found that /,, 7, and theory are in fairly good 
agreement in the early part of the discharge in all cases 


NITROGEN DISCHARGE H, =0 


DEUTERIUM DISCHARGE H, = 0 


(see fig. 8, A). The biggest discrepancy during this time, 
for deuterium with H,=0, can be explained by 
uncertainties in determining r, during the very early 
stages when the whole curve is very sensitive to the 
values of r, just after the plasma leaves the walls. 

The agreement between /, and /, in the early stages 
is sustained for more than two successive contractions 
in the case of deuterium with H, = 1480 oe. For 
deuterium with H, = 0, 7, begins to exceed /, soon 
after the first maximum contraction; /, falls as the 
central column contracts for a second time (see observed 
radius time curve) but /, continues to rise. For nitrogen, 
I, becomes much greater than /, before the first maximum 
contraction and the first dip in the observed current 
occurs well before the first maximum contraction of the 
central column. When the discrepancy appears the tube 
is presenting a lower inductance to the bank than it 
would if all the current flowed in the central region. 
This lower inductance represents a larger “‘effective”’ 
radius for the current channel and the radius/time curves 
bear this out. For nitrogen 7, is about 5 cm while r, is 
less than 0.5 cm at the first maximum contraction, and 
most of the current must flow in the outer regions of the 
tube. The observed peak current for a nitrogen discharge 
is 800 kA, while for deuterium it is 600 kA, although the 
radius of the luminous central column is smaller for 
nitrogen. 

The evidence from the current measurements could be 
consistent either with all the current flowing at radii in 
the neighbourhood of 7,, or with a current /, flowing on 
the outside of the central column together with other 


DEUTERIUM DISCHARGE H, = 1480 ce 


< 5 a < < 
"a i. aN Vig ae pene sea = po = 350, 
i = ie , 
600 6 ? ee 600 Vg 20kY 600F Vg 20kV 
(0)~1-5 110) F 
“at 110° amp/ sec of 4, ()~ 1-52 10"? amp / sec 47/4 @)~ 152107 amp/ sec 
soot Z 500 # 
iS : A yee $00 H 1480 o€ in 
7 STREAK (1 a 

_ 400h 7 CURRENT MONITOR ( 1, ) 400 F (1) ye royal STREAK (1, ) 
Fed fe 3 / & 
¥ 
S s CURR EK ey = CURRENT 
3 3 300 Sepals 


3S00F 


100 


MONITOR (Ie) / 
oe 


MONITOR (le) 


LO 


4 CM @ o 


RADIUS 
- 
RADIUS 


ELECTRICAL MEASUREMENTS (re) 


STREAK (r,) 


TIME 


Fig. 8 A (above) Comparison between the current/time variation measured directly (/,), and calculated from the streak measurements Ce) 
and according to theory for three different discharges. ‘7 
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surrents at radii >r,. However, the agreement between /, 
ind theory and r, and theory during the whole of the im- 
slosion phase in all three cases shows that the dynamical 
yehaviour of the discharge is consistent with a current /, 
lowing in the central column, even though, for the 
litrogen discharge, the total current becomes much 
rreater than /, before the first maximum contraction. 
After the first maximum contraction of the radius/time 
yehaviour of the central column is still largely consistent 
with a current /, flowing in it, and the most likely ex- 
ylanation of the discrepancy is the growth of additional 
surrent outside the luminous central column. 

The observation that /, can exceed /, is a general one, 
of which the cases shown are typical. The current dis- 
tepancy is observed to depend on the energy input to 
he gas, the compression ratio and the presence of high-Z 
mpurities in the gas. The growth of current outside the 
entral region has also been observed by Andrianov 


mal. [1]. 


Qualitative observations from the streak photographs 


At about the time current begins to grow outside the 
entral volume an increase in the background luminosity 
s observed on streak photographs, consistent with a 
natked increase in luminosity near the discharge tube 
valls. Similar effects have been observed by Burganov 
t al. [10] and Komelkov [11]. The onset of this light 


is measured from the outer edge of the luminous sheath. 


varies with energy input, compression ratio and the 
presence of impurities in the gas in the same way as the 
current discrepancy. In nitrogen discharges faint back- 
ground light appears during the initial implosion and 
becomes very intense at the first maximum contraction. 
In deuterium discharges at pressures > 100 w this light 
appears at or soon after the first maximum Contraction 
and is very much brighter than the central colum, which 
becomes completely indistinguishable. Figure 9 shows 
streak photographs in deuterium at different bank volt- 
ages (Vg) both with and without an axial magnetic field. 
The middle two pictures correspond to the curves for 
deuterium in fig. 8 and it is seen that with H, = 0 the 
onset of background light and the current discrepancy 
both begin at 1.5 usec; with HZ, = 1480 oe the light 
comes later, and /, and /, are in agreement for much 
longer. 

At initial deuterium pressure < 100 @ the central 
column is invisible, but background light still appears 
vety brightly at or soon after the first maximum con- 
traction (determined electrically). Even in discharges 
with 99.99 % pure deuterium the same general phenomena 
ate observed on the streak photographs. The addition 
of 21/,% of nitrogen to a deuterium discharge radically 
alters its character, and it behaves more like a pure nitro- 
gen discharge than a deuterium discharge in the respects 
described above. 


39 


H. A.BODIN, A. A.NEWTON &N. J. PEACOCK 


Interpretation of the results obtained from the streak 
photographs and electrical measurements 
The observations decribed above can largely be ex- 
plained by the following qualitative theory. At some 
time during the discharge the surface of the tube wall is 


vapourised by radiation from the central plasma, and 


current (in the same direction as the central current) 
builds up in the gas produced in this way, which becomes 
highly luminous due to the presence of oxygen and 
silicon. This model is consistent with the observation 
that wall effects are more pronounced and occur sooner 
if the energy input is raised, or if the content of impurities 
is increased or if the compression ratio becomes larger. 

A discharge in wall material, once initiated, will grow 
rapidly and the radiation need only provide a trigger. 
As the temperature of the wall surface depends on the 
radiation flux, the total intensity of radiation required 
can be quite small. It is easy to account for large radiation 
fluxes in heavy gases, but in deuterium the bremsstrah- 
lung can be shown to be negligible and the radiation 
could either come from impurities drawn in from the 
walls in the initial stages or from partially ionized 
deuterium, or possibly from some kind of turbulence. 
It has been shown (see, for example, McWhirter [12] and 
Kogan [13]) that the radiation from a deuterium plasma 
can be greatly increased by the addition of very small 
amounts of impurity. The wavelength of the radiation 
must be short, as wall effects have been observed when no 
visible light is detected from the central region; further- 
more, quartz is transparent to wavelengths > 2000 1S 

Wall effects can also be initiated by plasma which 
escapes from the central region and hits the walls. In 
some slower unstabilized discharges in light gases-there 
is evidence for plasma jetting out and causing a highly 
localised increase in luminosity at the walls (see, for 
example, fg. 9 fot _ = 15 kV, 42, =.0). For dli/dt = 
1.5 x 10!" amps per sec the onset of wall light is usually 
symmetrical (see fig.9 for Vg =20 and 25 kV and 
HT, = 0 and 1480 oe), and occurs in times rather short 
to be explained by plasma hitting the walls; in these 
conditions the radiation mechanism is much more likely. 

Further evidence for a discharge in wall material is 
the appearance of a second, highly luminous, sheath 
which contracts slowly in some conditions at low power; 
this sheath begins to form when the background of light 
from the vicinity of the walls appears. In other conditions 
the luminous region appearing near the walls is seen to 
move in slowly without forming a distinct sheath. The 
electron temperature in the discharge in vapourised wall 
material will be low and resistive effects important, lead- 
ing to a thick current sheath. 

These wall effects can explain many anomalies; for 
example, the appearance of the first “dip” in the current 
waveform well before the first maximum contraction in 
nitrogen discharges; the observation that the peak cur- 
rent for discharges in nitrogen is 30°/, greater than for 
those in deuterium; the decrease in peak current for 
deuterium discharges when H, is applied. 


The influence of wall effects on the compression heating 


The production of conducting material at the walls 
and the subsequent growth of currents in this material 
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screens off the inductive electric field from the central 
column and the rate of current rise in this column falls. 
The current in the central plasma is therefore limited and 
there is little further compression heating. At later times 
the central plasma becomes heavily contaminated, as 
observed by Andrianov et al. [1], Komelkov [11] and 
others. 

If attempts are made to increase the energy input to 
the gas by, for example, increasing the voltage, then 
wall effects appear earlier, and the heating becomes 
limited sooner, with negligible gain. Fig.9 shows 
evidence for this from streak photographs at Vg = 15, 
20 and 25 kV; it is seen that the background light comes 
sooner at higher voltages. Increasing the axial magnetic 
field or reducing the voltage both reduce the energy 
density in the plasma at the time of the first maximum 
contraction and are observed to delay the onset of wall 
effects. In fact it is found that over a wide range of con- 
ditions there is a limit to maximum energy input to the 
gas. There is evidence for this in fig. 6 and 7 which show 
that the maximum value of IW”, is almost the same with 
ot without H,. With H, — 0 more energy is put into 
the gas at the first maximum contraction, and there is 
no further heating. With AZ,, the energy input to the 
gas is less at the first maximum contraction, wall effects 
are delayed, and there is some further heating —suf- 
ficient to increase 7, to. about the value achieved with 
H,—0. In general the results suggest that, in the 
present experiments, as soon as the energy density in 
the central plasma exceeds a certain value, wall effects 
appear and prevent further energy input to the central 
column. In some experiments using very high voltages 
(120 kV) Andrianov et al. [14] found that large cur- 
rents flowed at the walls right from the start of the dis- 
charge. 

6. Conclusion 


In low-power x-pinches insufficient energy is put into 
the discharge to achieve a thermonuclear temperature, 
even if the initial implosion is followed by further com- 
pression heating; in this case such compression takes 
so long that the heating is limited by instabilities, im- 
purities from the walls and other cooling processes. In 
very fast pinches wall effects limit the heating at a much 
earlier stage and so in no case is a thermonuclear plasma 
produced. However, the ratio of particle to magnetic 
energy in the high-power case, before the heating be- 
comes limited, is very favourable for further compres- 
sion if wall effects can be eliminated. It seems probable 
that the type of effects observed in the present experiment 
are responsible for the failure to achieve a thermonuclear 
temperature in many of the high-power, fast-pinch ex- 
periments so far carried out. 
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ENERGY DISTRIBUTION OF NEUTRONS PRODUCED 
BY A THERMONUCLEAR REACTION* 


W.R. Faust AND E. G. Harris ** 


U.S. NavaL RESEARCH LABORATORY, WASHINGTON D. C., U.S.A 


The energy distribution of neutrons in a thermonuclear reaction is calculated. It is assumed that the reacting 
mixture of gases is in thermodynamic equilibrium and that the differential cross section for neutron production is 


isotropic in the center of mass system. 


The energy distribution of neutrons produced by a 
thermonuclear reaction has recently become of interest 
because of the possibility of producing such reactions [1] 
in the laboratory. Some calculations [2] have been made 
previously but were completely different in approach 
from the methods presented herein. 

In order to calculate the neutron energy distribution, 
it is assumed that the reacting mixture of gases is in 
thermodynamic equilibrium at a temperature 7 and that 
the equilibrium is not disturbed appreciably by the 
nuclear reactions which are taking place. It is also as- 
sumed that the differential cross section for neutron pro- 
duction is isotropic in the center of mass coordinate 
system of two interacting particles. 

It follows [3] from the assumptions made above that 
the number of pairs of ions with velocities in the range 
de,, dc, is 


Vi V5, (Wy Mio) 's 


yy (20k1)3 ~ 


xp zepinct + meat} dey dey (1) 


& = Boltzmann’s constant 

i = temperature 

C,, C, tepresents the velocities of ions of mass 
m, and m, respectively 

¥,, Yo = numbers of ions per unit volume of 
mass #7,, #, respectively 


where 


lif m, + mM, 
‘5 = 2 
Vie \2 if My = Me 


The differential rate for neutron emission into the 
solid angle d!2 —2 a sin6 d@ is obtained from (1) by 
multiplying by go(g) where g = |c, —c,| and a(g) is 
the appropriate reaction cross section, namely 
VV (My Mg) '!2 


il dQ 
€Xp \-serlmne? + mee) | dey dey 4 (2) 


The energy distribution of neutrons observed in a co- 
ordinate system at rest relative to the gas may be obtained 
from (2) by transformation of coordinates and subse- 


* This work was jointly sponsored by the Office of Naval Research 
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quent integration over variables extraneous to the prob- 
lem. The transformation involves a change of variables 
from c,, c, and cos@ to a new set g, G and & the 
neutron energy. It is convenient to perform the trans- 
formation in two separate steps; the first [4] 


My 


Cc = G’ A 
- a lily = 
My 
c, = G’ — — 9’ 3 
z my ® (3) 


cos.0 == cos 0 
where 7% = m, + mz 


The second transformation is 


g=g 
G2aG 
0) 


SE —(e}—# a) 


cos 0’ =m 
2Gu 


In the last equation of (4) & is the neutron kinetic 
energy measured in the laboratory coordinate system 
while # is the neutron velocity in the center of mass co- 
ordinate system and ~ is the neutron mass. The Jacobian 
transforming the original set of differentials to the new 
set da, dG, d£ 1s 


0(c,,C2,cos8) 9 (cy, cos) 9(g’, G’, cos 0’) 
SiG.) ole. Gyend) 9 ole Gs (5) 


Chapman and Cowling [5] have evaluated the first 
determinant on the right hand side of (5) and show it to 
bequal to unity, while the second determinant is 

0 (g’, G’,cos 0’) il 
BetG ie ame (6) 
wherenG ==s|.Gi. 

Limits of domaine for g, G and E are obtained from 
the original limits of the old variables subject to the 


restrictions implied by the transformations, thus it is 
found that 
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Since the limits of G are determined by E and u, the 
> integration must be performed first. By utilizing the 


Poobians found above, the new differential diceibadion 
an be written, 
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Integration of (8) over all directions of G and g yields 
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Next integration with respect to G and g in that order 
yields the desired energy distribution of neutrons, namely 
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For convenience the integration variables will be 
changed to EF, = (m,m,/2 m,) g?, and with the substitu- 
tion 


m : 
Te E a(QO+E, ) where « = (1 = i. equation (10) 
becomes, 
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EAR REACTIONS 


An excellent approximation to the actual form of the 
neutron energy distribution function (11) may be ob- 
tained by noting that FE ~ w IOrVwhile Bea? is-very 
much less than E. Under these condition equations (11) 
can be approximated as 
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It follows from equation (12) that the width of the 
neutron energy distribution at e } of its maximum value 
is (to terms of first order) 
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DESAY C Ts IN ENGLISH 


Plasma oscillations, (1), I. B. BERNSTEIN (Project Matterhorn, Princeton University, N.J.,U.S.A.) aNd S.K. TREHAN 
(New Delhi, India) Nucleat Fusion 1, 3 — 41 (1960) 


his article is the first part of a review of the theory of wave phenomena in plasmas. The basic kinetic theory is developed, taking 
Into account coulomb collisions. From this, magnetohydrodynamic equations are derived. It is demonstrated that it is legitimate to 
deal with a closed subset of these equations either in the limit of collision-dominated phenomena, ot in the limit where the effective 
phase velocities of the phenomena of interest are much greater than thermal speeds. In case of collision-dominated plasmas the 
theory of transport coefficients is discussed. 

These equations are then applied to an extensive treatment of small-amplitude wave phenomena in plasmas. A discussion of the 
dissipative effects on hydromagnetic waves is given. Hydromagnetic waves are also considered from the Chew, Goldberger and 
Low theory. Longitudinal and transverse oscillations in current-carrying plasmas are also discussed. 

Oscillations of a cylindrical plasma are considered and the phenomenon of ion cyclotron resonance is discussed. The possibility of 
radiation by plasma oscillations by a uniform sphere is exhibited. Some general results on the stability of longitudinal electron 
oscillations in non-uniform plasmas are given. A brief treatment of large amplitude electron oscillations is given and the breaking 
of these oscillations as a dissipative mechanism for the organized plasma motion is discussed. 

Part II of this paper, to appear in a future issue of this journal, will be devoted to the discussion of plasma oscillations directly 
from the kinetic theory. 


Rayleigh-Taylor instabilities of magnetically accelerated plasma, T. S. GREEN AND G. B. F. Nisert ( Atomic 
Weapons Research Establishment, Aldermaston, Berks., England) Nuclear Fusion 1, 42 — 46 (1960) 


High speed photography has been used to study the compression of a deuterium plasma by an axial magnetic field generated by a 
wide single-turn coil. In this experiment the inductance of the coil is small compared with that of the external condenser bank 
so that the plasma implodes at a constant acceleration of 5 - 1012 oms/sec?. Flutes which develop on the outer surface of the plasma 
are interpreted as magnetohydrodynamic Rayleigh-Taylor instabilities produced by the inward accelaration of the interface sepa- 
rating magnetic field and plasma. The observed growth rates are in agreement with simple theory. It is shown that the tatio of 
the diffusivities of magnetic field and momentum determines whether the plasma resistivity or its viscosity is the dominant mecha- 
nism for damping the short wavelength instabilities. At the temperatures and densities of this experiment the resistivity is primarily 
responsible for the damping and the wavelength of “maximum instability” calculated on this basis agrees with the observed value. 


A necessary condition for hydromagnetic stability of plasma with axial symmetry, C. MercriER (Group de 
recherches de 1’ Association Eurotom-C.E.A., Fontenay-aux-Roses, Seine, France) Nuclear Fusion 1,47 — 53 (1960) 


In this work a necessary condition is derived for the stability of a plasma with axial symmetry. This condition corresponds to 
displacements which are localized in the neighborhood of those surfaces of constant pressure on which the lines of force ate closed. 
This criterion reduces to Suydam’s criterion in the case of cylindrical symmetry. 


Rapid compression heating of a plasma in the linear, z-pinch H. A.B. Bopin, A. A. NEwTON AND N. J. PEACOCK 
(Atomic Weapons Research Establishment, Aldermaston, Berks., England) Nuclear Fusion 1, 54 — 61 (1960) 


The electrical characteristics of discharges with high and low initial rates of current rise are briefly described and the electrica 

measurements are used to determine the work done on the plasma by rapid compression. The particle and the magnetic energy in 
the discharge are discussed for both high and low initial rates of current rise (d//d7(O)) and for discharges with and without an 
axial magnetic field (H/,). It is shown that the gas is heated much mote efficiently both with and without 17, at high d//d7(0). 
For d//d¢(0) ~1.5 - 10! amps/sec (Al, = 0) 15° of the initial condenser bank energy goes into heating the gas and 10% into the 
magnetic field at the first maximum contraction; for d//d/ ~ 10" amps/sec (H], = 0) 11° goes into field and only 6.5% into the 
gas at this time. The energy input to the gas is found to become limited well before peak current, and the initial implosion is not 
followed by much of the subsequent isentropic compression heating which might be expected, so that the temperature never 
exceeds about 300 eV; it follows that most of the observed neutron yield (~10*) is not thermonuclear. Data obtained from the 
electrical waveforms and from high-speed photographs are presented and analysed, and it is shown that large currents build up 
outside the central plasma column after the first maximum contraction, Further evidence suggests that these currents flow in wall 
material boiled off by radiation from the central column. These wall effects account for the absence of further isentropic compression 
heating and also contaminate the central plasma column; they are such as to prevent any significant gain in heating by increasing 
the power input to the discharge. It is concluded that similar effects probably account for the failure to achieve a thermonuclear 


temperature in many high power pinch experiments. 


Energy distribution of neutrons produced by a thermonuclear reaction, W. R. Faust AND EG. EARRIsS 
(U.S. Naval Research Laboratory, Washington, De Te Ses) Nuclear Fusion 1, 62 — 63 (1960) 


The energy distribution of neutrons in a thermonuclear reaction is calculated. It is assumed that the reacting mixture of gases is in 
thermodynamic equilibrium and that the differential cross-section for neutron production is isotropic in the center of mass system. 
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RESUMES EN FRANGAIS 


Oscillations du plasma (1), |. B. BERNSTEIN ET S. K. TrEHAN (Project Matterhorn, Princeton University, Princeton, 


N. J., Etats-Unis d’ Amerique) Fusion Nucléaire 1, 3 — 41 (1960) 


Cet article est la premiére partic d’une étude d’ensemble portant sur la théorie des phénomenes ondulatoites dans les plasma. Les 
auteurs y développent la théorie cinétique fondamentale, en tenant compte des chocs coulombiens et parviennent ainsi a établir 
des équations magnétohydrodynamiques. Ils démontrent que l'emploi dun sous-ensemble complet de ces equations constitue une 
approximation légitime pour l’étude de phénoménes dans lesquels soit les chocs prédominent, soit les vitesses effectives de 
phases sont tres supérieures aux vitesses d’agitation thermique. Ils examinent la théorie des coefficients de transport dans le cas 
de plasmas ot les chocs prédominent. 

Ces équations sont ensuite appliquées 4 une étude approfondie des phénoménes ondulatoires de faible amplitude dans les plasmas. 
Les effects de dissipation sur les ondes hydromagnétiques font objet d’un examen critique. Les ondes hydromagnetiques sont 
également étudiées du point de vue de la théorie de Chew, Goldberger et Low. On examine aussi les oscillations longitudinales 
et transversales dans les plasma traversés par un courant. 

Les auteurs étudient les oscillations d’un plasma cylindrique ainsi que le phénoméne de résonance cyclotronique des ions. Ils 
montrent que les oscillations d’une sphére uniforme du plasma peuvent produire un rayonnement. Ils donnent quelques résultats 
généraux concernant la stabilité des oscillations longitudinales des electrons dans des plasmas non uniformes. 

L’atticle renferme une bréve étude des oscillations de grande amplitude des electrons et de l’interruption de ces oscillations en 
tant que mécanisme de dissipation pour le mouvement du plasma organise. 

La deuxiéme partie de article, qui paraitra ulterieurement dans un numéro de cette revue, sera consacrée a l’étude des oscillations 


du plasma, directement fondée sur la théorie cinétique. 


Instabilités de Rayleigh-Taylor dans un plasma soumis 4 une acceleration magnétique, T.S. GREEN ET 


G.B.E. Nisterr (Avomic Weapons Research Establishment, Aldermaston, Berkshire, Royaume-Uni) 


Fusion Nucléaire 1, 42 — 46 (1960) 


Les auteurs ont eu recours a la photographie ultra-rapide pour étudier la compression d’un plasma de deutérium par un champ 
magnétique axial créé par un entoulement large 4 spire unique. Dans cette expérience, l’inductance de la bobine est faible, comparée a 
celle des condensateurs extérieurs, de telle sorte que le plasma «implose» 4 une accélération constante de 5 >< 10! cm/sec?. De l’avis 
des auteuts, les stries qui apparaissent 4 la surface extérieure du plasma correspondent aux instabilités magnétohydrodynamiques de 
Rayleigh-Taylor, produites par l’accélération centripéte de Vinterface qui sépare le champ magnétique et le plasma. La vitesse de 
croissance obsetvée concorde avec la théorie. Les auteurs démontrent que, le rapport entre la diffusivité du champ magnétique et celle 
du moment determine si c’est la tésistivité du plasma ou sa viscosité qui contribuent surtout a amortir les instabilités de courte 
longeur d’onde. Dans les conditions de température et de densité de l’expérience, l’amortissement est dd surtout a la résistivité 
et la longeur d’onde «d’instabilité maximum», calculée a pattit de cette hypothése concorde avec la valeur effectivement obsetvée. 


Un critére nécessaire de stabilité hydromagnetique pour un plasma en symeétrie de révolution, CLAUDE 


MERCIER (Groupe de recherches de association E:uratom-C..E...A., Fontenay-aux-Roses, Seine, France) 


Fusion Nucléaire 1, 47 — 53 (1960) 


Pour ce travail, on établit une condition nécessaire de stabilité pour un plasma en symétrie de révolution. Elle correspond a des 
déplacements localisés au voisinage des surfaces 4 pression constante sur lesquelles les lignes de force magnétiques sont fermées. 
Ce critére se réduit a celui de Suydam en syméttie cylindrique. 


Chauffage du plasma par compression rapide dans la striction linéaire avec champ magnétique axial, 


H. A. B. Bopin, A. A. NEwron Er N. J. Peacock (Atomic Weapons Research Establishment, Aldermaston, Berk- 
shire, Royaume-Uni ) Fusion Nucléaire 1, 54 — 61 (1960) 


Les auteurs décrivent bri¢vement les caractéristiques électriques des décharges avec accroissement initial du courant rapide ou lent, 
les mesures électriques suivant 4 determiner l’effet de la compression rapide sur le plasma. Ils étudient l’énergie corpusculaire et 
lénergie magnétique inhérentes 4 la décharge pout des acctoissements de courant initiaux (d//d/(0)) rapides et lents, et pour des 
décharges avec et sans champ magnétique axial (H/,). Ils montrent que le chauffage du gaz est beaucoup plus efficace lorsque la 
valeur de d//dt(0) est élevée, qu’il y ait ou non un champ magnétique axial. Pour MLO Ss) SC UO IN IS WEE, == (0), IA ele 
Pénergie initiale fournie par les condensateurs servent 4 chauffer le gaz et 10°% sont communiqués au champ magnétique au moment 
de la premiére contraction maximum; pour d//f ~ 10" A/s (HI, = 0), 11% sont absorbés pat le champ magnétique et 6,5°% seule- 
ment par le gaz, au moment considéré. On constate que la quantité d’énergie fournie au gaz se trouve limitée bien avant que le 
courant n’ait atteint sa valeur maximum, et que l’échauffement par compression isentropique qui suit implosion initiale est beaucoup 
plus faible qu’on pourrait le prévoir, de sorte que la température n’est jamais supérieure 4 enviton 300 eV; il s’ensuit que la plus 
tare partie nombre observé de neutrons émis (~ 10) n’est pas d’origine thermonucléaire. Les auteurs exposent et analysent 
es données obtenues par l’étude des formes d’ondes électriques et des images photographiques prises A grande vi Ral 

que des courants importants se forment a l’extérieur de la Sina centrale ee ee eee parieioh LE ee 
données semblent indiquer que ces courants passent dans la matiére constitutive des parois chassée sous l’action du tayonnement 
provenant de la colonne centrale. Ces effets de paroi empéchent le chauffage par compression isentropique de se poursuivre et 
contaiminent également la colonne centrale de plasma; leur nature est telle qu’ils rendent impossible tout accroissement significatif 
de la température par une augmentation de la quantité d’énergie communiquée a la décharge. 

Les auteuts concluent que c’est ptobablement a cause d’effets analogues que l’on n’a pas pu obtenir de températures thermonuclé- 
aites, au couts de nombreuses expériences de striction du plasma mettant en jeu de grandes énergies., 


Distribution énergétique des neutrons produits au cours d’une reaction thermonucléaire, W. R. Faust Er 
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BE. C. Harris (U.S. Naval Research Laboratory, Washington, D.C., Etats-Unis dL Amérique ) 


Fusion Nucléaire 1, 62 — 63 (1960) 
Les auteurs de la monographie font le calcul de la distribution énergétique des neutrons au couts d’une teaction thetmonucléaire 


Ils partent de ’hypothése que le mélange gazeux réagissant se trouve en état déquilibre thermodynamique et que la section efficace 
différenticlle relative 4 la production de neutrons est isotrope dans le systeme du centre de gravité 


AHHOTAIMNM IO-PyCCKU 


OcunnnaAnMK na3zmEr (I), UW. B. Bepuurein u C. K. Tpexan (MarrepxopucKkuit npoext, IpuncronucKkui yuu- 
sepcurer, IIpuncron, How Jacepcu, CILIA) AmepHpim cuutTe3 1, 3-41 (1960) 


Wanna cTatbA ABIIAeTCAH TMepBol uwacTbro O6380pa TeOpuyt BOMHOBbIX ABIeHMi B MmmagMaXx. PaspaooTaHa OCHOBHaA 
KuneTM4ecKand TeOpMA C yUeTOM KYJIOHOBCKMX CTONKHOBEHMM. M3 9TO TeOpuM BLIBORATCA MarHMUTOrMAposAMHaMMueckMe 
yPpaBHeHuaA. TlOKa3aHoO, YTO MOXKHO OOpalsaTben C B3AMKHYTOIt NOACMCTeMOM 9TMx ypaBHeHMult unm B Mpenenax ABIIEHMIM 
AOMMHMPYHOUIMX CTOIKHOBeHMM UNM B cIyuYaAX, KOTZa scbcdbeKTMBHbHIe CGa3so0BEIe CKOpPOCTM MHTepecyroulMx ABIEHMIt 
SHAUYMTeIbHO OOJbUIe TeNIOBBIX CKOpOcTes. pm pacecmMoTpeHum OMMHMPYyIOUMxX MWas3MeCHHbIX CTOJIKHOBeHMi1 O6Cyx- 
AaeTCA TCOPMA TPAHCIOPTHbIX KOS MUMeHTOB. 

SaTeM 9TM ypaBHe€HMA MpPMMeHAIOTCA kK MUHTCHCMBHOMY MCCJICNOBAHMIO BOJIHOBbIX ABNeCHMM Cc Manow aMMIMTy_ON B 
ma3sMax. IIpmMBOAMTCA paccMoTpeHMe Bompoca O BIMAHMM paccesHNiA Ha TupOMarHMTHEIe BONHEI. TugpomMarHuTHpre 
BOJIHBI PaCCMaTPMUBAWTCA Tax7Ke C TOUKM 3PeHMA TeopuM Ubro, Tompaoepra u Jloy. Odcy2xmarwrTrca Tak2Ke MpOOMbHHIE 
M MOMepeyHbie OCUMIINIAWUMM B TOKOHeECYUIMX TMa3sMax. 

PaccMaTpMBAIOTCA OCHMINAUMN WUMIMHapMUecKOIt mma3sMBI UM o0cyxKyaercH ABMEHMe MOHHO-IMKJIOTPOHHOTO pPe3s0HaHCa. 
Tloka3aHa BO3MO?KHOCTh NOMYYeCHMA MSNYYeHMNA B pPesybTaTe OCUMANAUMM WasMbI B OAHOPOAHOIt eqdepe. I[pusonzatca 
HEKOTOpble OOMUMe PeSYIbTATbI CTAOMMBHOCTM MPOAONLHHIX 9NECKTPOHHbIX OCUMMMAUM B HEOMHOPOAHbIX Mmla3MaXx. 

Kpatko o0cy2xaeTcCA TpaKTOBKa SMECKTPOHHBIX OCUMINAWUMM BbICOKOM AMIIIMTYbI a Takx#e MepepbiB 9TUMX OCIMI- 
MAWMMm KaK Me€XaHM3M paccesHMA ANA OpraHM30BaHHOTO LBYM2KeHMA WN1agMbI. 

Mactp Il stom crarbu, HaMeyeHHaH K OMyONMKOBAHMIO B MOCMeAyIOULeM HOMepe AAHHOTO 2WypHaNa, OyeT MOCBAUIeEHA 
paccMOTpeHMIO OCUMIIAUMM WiasMbIl HemocpesCTBeHHO C TOUKM 3PeHMA KMHETMMeCKOM TeOopMn. 


Heycromunsoctn Peinu-Telinopa B nia3mMe, yckopsemoit B MarHuTHoM none, T. C. [pun u Ix. B. ©. Huobnerr 
UccnredosarenrvcKxut wexTp no aTomHOMY Opyseu10, OnrdepmacTon, Bepxrwup, Axn2nuUn) 
AmepHpm cuntes 1, 42-46 (1960) 


TIpm M3sy4eHuuM cmaTMA WeuTepMeBonw MNa3MbI C MOMOWbIO aKCMaNIbHOrO MAarHMTHOLO MONA, OOpasoBaHHOTO WIMpoKoM 
KaTyUuIKOM C OAHMAM BMTKOM, ObIJIa MCMONb30BaHa BbICOKOCKOpOcTHaA chboTorpadua. B faHHOM 9KCTepuMeHTe MHAYK- 
TUBHOCTB KaTYUIKM HaCTOJIbKO He3HauMTesbHa MO CpaBHeHMIO C MHAYKTMBHOCTbIO BHeLIHeM GaTapew KOHAeCHCaTOpPOB, 
4TO Tiwa3Ma TeyeT C MOCTOHHHBIM ycCKOpeHMem B 5X 10!" cm/ceK®. 

SaBuxpeHMA, oOpasyeMble Ha BHeUIHe/ MOBeEPXHOCTM MasMbl, paccCMaTpMBaIOTCH KaK MarHuTorMp~pomMHamuueckaA 
HeycTomumBoctp Pelinu-Teitnopa, BbIgsBaHHaH BHYTPCHHMM yCKOpeHMeM IOBePXHOCTM, pa3semArourleM MarHMuTHOe Tose 
M UsasmMy. Hadmtogaemoe BO3pacTaHMe cKOpocTeit cooTBeTcTBYyeT Mpocrom Teopuu. TloKka3biBaeTcA, UTO OTHOWIeHMe AMCb- 
dy3suu MarHuTHOTO MONA M WMITyIbCa OMpenenseT, ABIHeTCA 2M COMpOoTMBNeHMe MWNa3MbI MUIM ee BABSKOCTb JOMMHM- 
PYOUIMM MeXaHU3MOM NIA TaulieHu4A KOPOTKOBONHOBOM HeycTolMunusoctu. Iipu tTemmeparypaX M MNOTHOCTAX 9TOTO 
SKCNepuMeHTa COMpOoTMBNeHMe ABNIHeTCH OCHOBHOM MpM4MHOM raulieHMA M DNMHA BONHbI «MAaKCMMANbHOM HeYyCTOM- 
4MBOCTM», MOACUMTaAHHaAAH Ha YSTOM OCHOBe, COOTBeETCTBYeT HaOmIOMmaeMOM BeIMUNHe. 


Heo0xocquMEM KpPMTepuM TMAPOMAaTHHTHONM CcTaOUTbHOCTH AIA WN1a3MbI C OCeBOM CuMMeTpHen, Kron Mepcpe (Mcciie- 


OosaTenvocKana zepynna Accoyuauuu Esparom — Komuccapuar no aromnot snepeuu — MDonrene-o-Pos, 
PDpanwugn). AmepuHpm cuute3 1, 47-53 (1960) 


B stom padote yeTaHaBlMBaeTcaA HeEOOXOAMMOe yYCNOBMe cTaOMMIbHOCTU AIA WWasMbI c OCeBOM cuMMeTpMenm. OHO CooT- 
BeTCTByeT TICEPeCMeCLICHUAM, NOKAJIM3SOBaHHbIM IO coceyxCTBy Cc MOBePXHOCTAHMM MOCTOAHHOTO aBJIEHUA, Ha KOTOPBIX 3a- 
MBIKaIOTCA CMMIOBbI€ MarHMTHbie JIMHMM. OTOT KPMTepMmt cBoqMTCA K KPMTepMio Camema ANA cnyyaA UMNMAApMAecKOU 
CMuMMeTpun. 


BsicTtpoe cmaTNe HarpeBaeMOM 1a3Mbl B IMHeEMHOM % TMHYe, X. A. B. Bogun, A. A. Hproron u H. XK. MuKox 
(AccredosaTerocKut wexTp NO aTOMHOMY OpYyocuI0, OndepmacTon, Bepxwup, Anenrun) 
AgepHpm cuuytTes 1, 54-61 (1960) 


B w#oKsage WaeTcH KpaTKOe ONMCcaHMe SJIEKTPMAeCKUX XapaKTePUCTUK Pa3sPAAOB C BbICOKOM MU HU3KOIt MepBOHAdAIb- 
HOM MHTCHCMBHOCTbIO POCTa TOKa M 9JIEKTPMUYeCCKUX USMEPeHMM, MCNONb3SYeMbIX AIA OMpeeueHMA BBEINOMHEHHOM pa- 
6OTbI MO ObICTPOMy CaTKIO wasMbIl. PaccMaTpMBaeTcA BotIpoc O06 SHeEprMM wacTMIL mM MarHMTHOTO MONA B paspaye 
ANA cnyuaeB C BbICOKOM U HU3KOM MepBOHAYAaIbHOM MHTCHCMBHOCTS1O pocta TOKa ( d//d+(0)) u ANA paspaAnOB UpM Hamu- 
4MM MU B OTCYTCTBMM aKCMaIbHOTO MarHMTHOrO Tosa (H,). TloKa3biBpaetca, TO ropa3syqzo0 Gomee sddbexTuBHOe HarpeBa- 
Hue Ta3a MpoMcxoquUT Kak pM HanMuuu, Tak MU B OTCyTCTBUe (H,) pm SonbuIOM BenmuMHe (d//dt (0). a dl/dt (0) ~ 1,5x 10" 
ammep/cex. (H,=0)15% nmepBpoHayanbHOM 9HeEpruu KOHAeEHCaTOpOB MpeT Ha HarpeBaHMe rasa u 10% wa mMarHuTHoe moe 
mpu MepBOM MaKCMMaJIbHOM coKpaljeHun; asa d//dt~ 10" amnep/cex (H,= 0.) 3a Tor me nmepuon BpemeHu 11% wager na 
nome uu wMulb 65% Ha ra3. Tlepeqaua 9HeEprum ra3a OKAa35IBAaeTCAH OTpaHMYeHHOM B3aqonro AO AocTuxKeHMUA MAKCMMamb- 
HOTO 3HAYeCHMUA TOKA, UM NepBoHadaNbHOe pe3sKoOe veRaTMUe HE COMpOBOwAaNOCcb MOcNeAyIOULMM HarpeBOM M309HTpOMMU- 
uecKOTO C#KaTMA, KOTOPOrO MOXXHO ObINO ObI OMMNLATb, Tak 4YTO TeMMepaTypa He MpeBHIuaNa MOpAgAKA 300 9B. OrTcroma 
emeaqyeT, 4YTO OONbUIMHCTBO HaOIOfaeMOTO BbIxXOfa HeEUTpoHOB (~ 10°) He TepMOoARepHOrTO MmpomcxompeHna, TIpusBo- 
ATCA U aHanMusupyioTcaA WaHHble, MOMyYeHHbie 10 GopMaM 9s7eKTpM4eCKUX BOMH M MO BbICOKOCKOPOCTHbIM GboTorpaduamM; 
M0Ka3biBaeTcA, UYTO Mocne NepBoro MAaKCMMAaIbHOTO COKPallljeHMA BO3SHUKAIOT CMUNbHbIe TOKM BHe IeHTpasIbHOTO UIHypa 
mma3Mpl. JlanbHelMume BaHHbie MOATBep*RAalwT, ATO BTM TOKM poTeKaroOT B MaTepMase, BbIMapeHHOM M3 CTeCHOK M3Iy- 
YeHUeM IECHTpaJIbHOrTO WiHypa MitasMsIl. BTO BIMAHMe CTEHOK OOBACHHET OTCYTCTBUe [aNbHeMulero UZO9HTPOMMUYeCKOLO 
HarpeBa c*#KaTMeM M MPMBOAMT Take K 3arpPASHEHMIO WeHTpaNbHOTO MwasMeHHOTO UWiHypa; OHO TaKOBO, 4YTO MpesoTBpa- 
ujaeT m1o60e 3aMeTHOe yBesIMy4eHNUe HarpeBa pM poctTe MorsioulaemMom 9HeEprMu B pa3spAne. B 3aKN1O4eHMe TOBOPUTCA, 
UTO MOZOOHbIe BScbcsbeKTbI, BEPOATHO, ABJAIOTCA MpMuMHONW Heyywauy MOCTM*eHUA TeMIMepaTypbl TePMOAMePHbIX peak 
BO MHOTMX S9KCMepuMeHTax TIO MMHYY BbICOKOM MOUIHOCTH. : 


Pacnpezenenue 9HeEprMu HeMTPGHOB, OOpaso0BaHHbIx B pe3yibTaTe TepMOATepHon peakuuu, Y. P. PaycT u E. 


Jix. Xappuc (Mopcexaa uccnedosarenvcKkan naboparopusx CLA, BawuneTrou, Oxpye Konymbus, CLITA) 
Anepupm cunte3 1, 62-63 (1960) 


TloxcuntTpipaeTcaA pacupemeneHue 9HEPrMM HesITPOHOB pu TepMoAepHoOM peakuMu. Wpennomaraetca, UTO pearupyio- 
ulaHd cMecb ra30B HaXOZMTCA B TePMOZMHAMMYeCKOM PaBHOBecuM uM uTO AMcbcepeHI{MalbHOe NMoMepedHoe ceueHue AIA 
reHepallMm HeMTPOHOB ABIIAETCH UZOTPOMHbIM B CMCTeMe WeHTpa Macc. 
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Oscilaciones del plasma (1), I. B. BERNSTEIN Y 
Princeton, New Jersey, los Estados ¢ nidos de Norte America) 


RESUMENES EN ESPANOL 


S. K. Trenan (Proyecto Matterhorn, Universidad de Princeton, 
Fusién Nuclear 1,3 —41 (1960) 


Este articulo constituye la primera parte de un estudio sistematico de la teoria de los fenédmenos de tipo ondulatorio en los plasmas. 
Los autores desarrollan la teoria cinética en que se basan, teniendo en cuenta las colisiones coulémbicas, y derivan de ella ecuaciones 
magnetohidrodinamicas. Demuestran que cabe emplear un sistema cerrado de ecuaciones, tanto hasta el limite en que predominan 
las colisiones como hasta el limite en que las velocidades de fase efectivas de los fendmenos de interés son muy superiores a las 
velocidades térmicas, Comentiin la teorfa de los coeficientes de transporte pata el caso de plasmasen los que predominan las colisiones. 
A continuacién aplican estas ecuaciones a un tratamiento matemético muy extenso de los fenomenos ondulatorios de pequces 
amplitud que tienen lugar en un plasma. Analizan el efecto de disipacién en las ondas hidromagneticas, y estudian éstas también 
desde el punto de vista de la teoria de Chew, Goldberger y Low. Discuten las oscilaciones transversales y longitudinales en plasmas 
portadores de corriente. ~y 
Seguidamente estudian las oscilaciones de un plasma cilindrico y discuten el fendmeno de resonancia ciclotronica de los i 
Muesttan la posibilidad de que una esfera uniforme del plasma emita radiaciones por oscilacion. Dan algunos resultados de caracter 
general sobre la estabilidad de las oscilaciones longitudinales de electrones en plasma no uniformes. 
También comentan un tratamiento abreviado de las oscilaciones de gran amplitud de los electrones y la interrupcion de estas oscila- 
ciones como mecanisiro de disipacién del movimiento de un plasma organizado. ' 

La segunda parte del trabajo apatecera en otro numero de esta revista y tratara del estudio de las oscilaciones del plasma aplicando 


directamente la teoria cinética, 


Inestabilidades de Rayleigh-Taylor en un plasma acelerado magnéticamente, T. 5. GREEN y G. B, F. NiBLETT 


(Atomic Weapons Research Establishment. Aldermaston, Berkshire, Reino Unido) Fusion Nuclear 1,42 —46 (1960) 


Los autores han estudiado la compresién de un plasma de deuterio por un campo magnético axial producido por un solenoide 
de una sola espita de didmetro grande. En este experimento la inductancia del solenoide es pequena compatada con la de bateria 
de condensadores extetiores, de manera que la implosién del plasma se produce con una aceleracién constante de5  10!2 cm/s?. 
En la superficie exterior del plasma se producen irregularidades que se consideran como inestabilidades magnetohidrodinamicas 
de Rayleigh-Taylor, producidas por la aceleracién centripeta de la superficie de separacion entre el campo magnético y el plasma. 
Se ha observado que la velocidad del proceso responde a suposiciones tedricas sencillas. En el articulo se demuestra que la relacion 
entre la difucibilidad del campo magnético y la cantidad de movimiento determina si la resistividad del plasma o su viscosidad 
es el mecanismo dominante pata amortiguar las inestabilidades de corta longitud de onda. En las condiciones de temperatura y 
densidad en que se efectué el experimento, Ja resistividad constituye el principal factor de amortiguamiento y la longitud de onda 
de la “inestabilidad maxima” calculada sobre esta base coincide con el valor obsetvado. 


Un criterio necesario de estabilidad hidromagnética para un plasma de simetria axial, CLAUDE MERCIER 


(Groupe de Recherches de ? Association Euraton — C.F. A. , Fontenay-aux-Roses, Seine, Francia) 


Fusié6n Nuclear 1, 47 — 53 (1960) 


En la memoria, el autor establece una condicién necesatia fara la estabilidad de un plasma de simettia axial. Esta condicion cortes- 
ponde a desplazamientos que quedan localizados en las proximidades de aquellas superficies a presi6n constante en las que las 
lineas de fuerza son cerradas. En el caso de una simetria cilindrica, este criterio queda reducido al de Suydam. 


Calentamiento de un plasma por compresién rapida en la constriccién lineal con campo magnético 


axial, H. A. B. Bopin, A. A. Newron yN. J. Peacock (Afomic Weapons Research Establishment, Aldermaston, 
Berkshire, Reino Unido) Fusion Nuclear 1, 54 — 61 (1960) 


Los autotes describen sucintamente las caracteristicas eléctricas de descargas efectuadas con indices altos y bajos de crecimiento 
inicial de la cortiente, y examinan las mediciones que se utilizan para determinar los efectos que tiene sobre el plasma una com- 
presion rapida. Se estudia la energia de las particulas y la energia magnética resultantes de la descarga con respecto a indices altos 
y bajos de crecimiento inicial de la corriente (d//dt(0)) y con respecto a descargas con campo magnético axial (/7,) y sin él. Para 
valores altos de d?/d7(0) el calentamiento del gas es mucho mayor, tanto con un H, como sin él. Para un d//di(0) de aproximadamente 
1,5 X 10% amp/s (campo HT, = 0), el 15 por ciento de la energfa inicial de la bateria de condensadores pasa a calentar el gas y un 
10 por ciento al campo magnético al producitse la primera contraccién maxima; para un d//d¢ ~ 101! amp/s (H/, = 0), el 11 por 
ciento pasa al campo magnéetico y sdlo un 6,5 por ciento al gas, La energia absorbida por el gas queda limitada bastante antes de que 
la corriente aleance su maximo, y la implosién inicial no va seguida de un calentamiento por compresién isentropica de la magnitud 
que podia esperarse, de manera que la temperatura no rebasa nunca los 300 eV aproximadamente; de ello-se deduce que la mayor 
parte de Ja cantidad de neutrones emitidos (~ 10°) no es de origen termonuclear, Los autores exponen y analizan datos obtenidos 
de las formas de onda y de fotografias ultrarapidas, indicando cémo después de la primera contraccion maxima se forman grandes 
corrientes fuera de la columna plasmatica central. Otras pruebas sugieren que esas cortientes pasan a material de las paredes des- 
ptendido pot la accién de las radiaciones de la columna central. A este efecto de las paredes se debe el que no se produzca un 
ulterior calentamicnto por compresién isentrépica, asi como la contaminacion de la columna plasmatica central; su naturaleza es 
tal que impide todo aumento importante del calentamiento cuando se aplica mayor energia a la descarga. Los autores llegan a la 
conclusién de que es a efectos andlogos a los que probablemente se debe la imposibilidad de alcanzar una temperatura propia 
de las reacciones de fusi6n en muchos experimentos realizados sobre la constriccién aplicando clevada energia. 


Distribucién energética de los neutrones producidos en une teaccién termonuclear, W. R. Faust y 
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E. G. Harris (U.S. Naval Research Laboratory, Washington, D. C., Los Estados Unidos de Norte America 5 


Fusion Nuclear 1, 62 — 63 (1960) 


Los autores calculan la distribucién energética de los neutrones producidos en una reaccién termonuclear, Al efectuar el calculo 

: ays Gj anit NEG On « 2 
se supone que la mezcla reactiva de gases se encuentra en equilibrio termodinamico y que la seccion eficaz diferencial correspondiente 
a la produccién de neutrones es isotrépica en el sistema del centro de gravedad. 
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HHP®OPMANNA AH ABTOPOB 


3a cmpaBkKamnu cmenyer oOpaljaTbca mo ampecy: «Pe- 
AaKTOpy WypHata «AWepHEt cures», MexxazyuHapoy- 
Hoe areHTCTBO TIO ATOMHOM sHeprum, Bena I, ABcrpua.» 


Hopayok odopmuerua pyKonnceli 


AxBTopaM lpeymaraerca CHEAYIOUIMM MOpAyOK ocbopm- 
JIGHUA PykKoMcen: 


a) Dopma pyxonucn. Pykommcb MoxKeT ObITB mpen- 
CTaBweHa Ha aHTIUMCKOM, pyccKom, cbpaHiy3cKom 
UNM UCMaHCKOM A3bIKe. TekcT pykKonmcu ouKeH 
ObITh OTIIeuaTAaH HA WMLUryuyet mMammaKe uepe3 Oba 
UNU TPU UATepBaNd, UMETh WuUpOKUe NONA U GOBLITS 
mipeyqcTaByleH Ha OymMare xopoulero KauecTBa B Tpex 
IK3SEMNAADAX; 


b) Ypasnenua ONRKH’I SbITH MocmeqOBaTeNBDHO TIpo- 
HYMepoBaHbI apabcKumu ywuq@panu B@ cKobKax «( )», 
KOTOPBI€ JOJXKHBI OBIT pacMOOReHBI cupaBa OT 
ypaBHennua. Ilo BosMoxHOcTM ypaBHeHua 6yzyT Ja- 
BaTbCA B OAMHAaKOBOM cbopme B nMepeBowax Ha BCe 
ASBIKU. IlosTOMy 2KeAaTeMbHO UCHONb30BaTL B ypaB- 
He€HMAX TOJIbKO Mex] yHapOyHbIe 3HAKM M M30eraTb 
TaKRe COKpaljeHuu COB, KOTOPbIe MOHATHLI TOIbKO 
B AaHOM A3bIKe (T. e. Enap, Xmaxe MU T. Z.). Cue- 
AyeT UCHONIb3OBATb MpocTHle 3HaKM, WaB MUX MOACHe- 
HMe B TCKCTe; 


c) Pucyunu u cxemot. 


1) Kaxyem pucyHok m cxema 2OmxKHEI OBIT TaHbI 
Ha OTZe€bHOM cTpaHuye M MpoHyMepoBaHBI m0- 
CHEFOBAaTeIbHO apabcKumu wWuqdpamu, Hanp., 
«Pue. 1»; 


2) Hagmmeu KO BceM PMCyHKaM M CXeMaM OIKHDI 
ObITh ZAHbI Ha OTOeNDUOM NUCTe YM COOTBeTCTBeH- 
HO IIPOHYMepoBaHB!; 

3) HanpaBiaTb doTocHumMKM cmeqyeT TONbKO B cIy- 
wae KpalinHei HeoOxoqumMocTM B HMX; 


4) 2KematenmbHO UpescTraBnaTb uepTexKu, BBINON- 
HeHHbIe YEPHOM TyWIbIO Ha OeOM BOCKOBKe; 
5) Ha Bcex pucyHKax 2KemraTenbHO MCHOUb30BAaTb 


Mex] yHapOwHbie ycJIOBHbIe 3HaKN, OnpenzemeHua 
STUX SHAKOB AOJZKHbI BKIOUATbCA B DOApUcy- 


HOYHbIe HaanMcu. OTO oberuuT mepesox mare- 
pMasIoB; 

d) CHocku foIKHbI OBITS MPOHyMepoBaHBI apabcKumu 
UUPpamu MOCcUeqOBAaTeNbHO U yKaB3bIBATbCA B TeK- 
cTe. SaTeM BCe CHOCKM JOJIKHbI ObITD TAaHbI Ha OT- 
QeIbHOM mMCcTEe; 


e) HammenoBaHua ciupaBounomt mMTrepatTypbl osKHbI 
ObITb MpOHYMepoBaHbI apabcKumu WUGpanu B KBaT- 
PpaTHbIxX cKOOKax. ITlomHbIe HaMMeHOBAaHMA BCex 
VUCHONIb30BAHHbIX MUSTAHMt TOIKHbI ObITh mpMBere- 
HbI B OYeEPEFHOM MOPAKe Ha OTAeIbHOM crpanmne 
mo cmeqyrouemy oOpasiy: 

«{1] UEXOB, K., Aromuan snepeun 14, 137 (1959)» 
umu [2] ABAHOB, ®., busuxa naazmvor (AH CCCP, 
1956)»; 


f) TaOnuubl OKI ObITS MpOHyMepoBaHbI mocme- 
WOBAaTeIbHO PUMCKUMU YUQOpamu. Hagmmen K Tab- 
IMUaM JOJ2KHbI ObITh MeCPCUMCMeHbI HA OTTEIbHOM 
auctre. Ilo Bo3mMoxxHOCTM B TabmMIAax cmenyer uc- 
NOJIb30BaTb Mew TyHApOAHbIe YCMIOBHbIe 3HAKU, 3Ha- 
YeHUA KOTOPBIX OJIZKHbI ObITb FAHHbI B TeKCTe 
Hagnueemw K TaOsuijam. PasmMepbl JOU2KHBI TAKHKe 
WaBaTbesA B TeKcTe HagnMcem K TadmMIaM; 


g) K pykonmcu KaxxOM cTaTbu OxKHO ObITB mpu- 
JIO#KeHO KpaTKoe pesiome (IpMMepHO 300 cm0B). Pe- 
310Me OyayT NepeBoxMTHcaA Ha padoune ASHIKM ATéHT- 
CTBAa U BKJIOUATbCA B KypPHa. 


Iipumeuanne: 


Bce pucynkKu, cxemMbl, TAOHNUYSI, AUCToL CO CHOCKAMU 
K TeKCTY U ChpaBOUHOuU NUTepaTYypOoK wHNeobxodU.mo 
NpPeOCTABAATS B Tpex 2K3eMnNNApAX. 


TIpm mpeycTraBmeHum pykKommen aBTop yommKeH yKa- 
3aTb, KOMy M NO KakOMYy agpecy OKI OITH Ha- 
WpaBJIeHbI TpaHkKM cTaTbu. BuaHKM Ha 3aKa3 ONON- 
HMTCIIbHOTO KOJIM4CCTBA SKSCMMIAPOB OyayT HampaB- 
JIATbCA OJHOBDCMCHHO c TrpaHKaMmi. 


B cayuae NPMHATMA pykKOoTMcu Kam TOMy aBropy Oy- 
net OecnmaTuHO HalpaBNaTbeaA 50 9K3eEMNONIAPOB TMMO- 
rpadckoro u3qaHMAa ctatew. JOMNONMHUTeIbHbIe 9K3eM- 
IIAPbI credyeT 3AKaA3vLIBATb TPM BOsBpaljeHum B 
ATeHTCTBO TpaHOK CTAaTbM. 
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La correspondencia debe dirigirse al «Redactor de la revista 
FUSION NUCLEAR, Organismo Internacional de Energia 
Atémica, Viena I, Austria». 


Presentation de los Originales 


Al preparar los originales deberdn tenerse en cuenta las siguien- 
tes normas: 


a) Texto de los articulos: Los articulos podran redactarse en 
espafiol, francés, inglés o ruso. El texto, con dos copias, 
se enviarad mecanografiado a doble o triple espacio, con 
amplio margen y en papel de buena calidad. 

b) Todas las ecuaciones se numeraran con nimeros arabigos 
colocados entre paréntesis a la derecha de cada ecuacién. 
En la medida de lo posible, las ecuaciones se reproduciran 
exactamente en todas las traducciones; por ello, conviene 
emplear los simbolos reconocidos internacionalmente y 
evitar el uso de abreviaturas que sdlo tengan sentido en un 
determinado idioma (por ejemplo, Vinegs tmax» etc.).. De- 
berdn utilizarse simbolos simples cuyo significado se defi- 
nira en el texto. 


c) Figuras: 

1) Las figuras se presentardn en hojas aparte y se numera- 
ran con nimeros arabigos; por ejemplo: «Fig. 1»; 

2) Los textos que acompafien a las figuras se agruparan en 
una hoja aparte con la numeracién correspandiente; 

3) Sélo se enviaran fotografias cuando sea indispensable; 

4) Se recomiende que las figuras se hagan en dibujo lineal, 
utilizando en lo posible tinta negra y papel de dibujo 
blanco; : 


5) Para facilitar el trabajo de traduccidn conviene que en 
las figuras sdlo se empleen simbolos reconocidos inter- 


nacionalmente; estos simbolos se definiran en el texto 
que vaya al pie de la figura. 


d) Las notas de pie de pagina se indicaran en el texto con 
numeros arabigos y se agruparan en una hoja aparte. 


e) Las fuentes bibliograficas se indicaran en el texto con nume- 
ros arabigos colocados entre corchetes, y se agruparan en 
una hoja aparte en orden numérico. Las menciones se haran 
de conformidad con los siguientes ejemplos: «[1] GRA- 
NADOS, A. J., Rev. esp. fis. 17, 483 1958)» o «[2] JONES, 
L. M., Plasma Physics. (XYZ Book Co., New York, 1957) 
59». 


f) Los cuadros se numeran con nimeros romanos y los textos 
explicativos se agruparan en una hoja aparte. Siempre que 
sea posible se utilizaran en las columnas simbolos aceptados 
internacionalmente, dando su significado en los textos expli- 
cativos. 


g) Junto con el original se enviara un resumen (de unas 300 
palabras). El resumen se publicara en la revista en los cuatro 
idiomas de trabajo del Organismo. 


Nota: 


Deberan enviarse un original y dos copias de las figuras y de 
los cuadros, ast como de las hojas separadas que contengan las 
notas de pie de pagina y las fuentes bibliograficas. 

Al enviar el original se indicara el nombre y la direccién 
exacta de la persona a la que deberdn remitirse las pruebas de 
imprenta. Con las pruebas se remitir4 un formulario para en- 
cargar las separatas. 

Los autores de los articulos publicados en la revista recibiran 
gratuitamente 50 separatas. Si desean recibir una cantidad 
mayor habrdn de indicarlo al devolver las pruebas de imprenta. 
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